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Knickung gerader Stabe durch Druck und konservative Torsion *. 
Von M. Beck. 


1, Einleitung. L. Euler hat im Jahre 1757 als Erster die axiale Druckkraft berechnet 1; 
unter der ein schlanker Stab ausknickt. Ein seitliches Ausbiegen kann auch bei einem Stab 
auftreten, der ein Drehmoment zu ithertragen hat. Die kritische Belastung eines gedriickten 
und tordierten Stabes hat fiir kreisférmigen Querschnitt zuerst A. G. Greenhill 2? untersucht, 
wahrend H. Grammel * fiir die Lagerung V von Abb. 4 eine Lisung fiir beliebige Querschnitte 
perechnet hat. Diese Arbeiten beruhen auf der Annahme, da8 der Momentvektor beim Aus- 
knicken seine urspriingliche axiale Richtung beibehalte. Das merkwiirdige Ergebnis, daB bei 
den Lagerungen III und IV in Abb. 4 ein beliebig kleines Moment den Stab zu knicken vermag, 
‘egte H. Ziegler * zu einer genaueren Uberpriifung des belastenden Momentes an. Dabei zeigte 
sich zunachst, daB ein axiales Moment nichtkonservativ ist. Nichtkonservative Stabilitats- 
probleme miissen aber mit dem kinetischen Stabilitatskriterium (Untersuchung der kleinen 
Schwingungen um die Gleichgewichtslage) behandelt werden, wahrend das von Greenhill und 
Grammel verwendete einfachere statische (Bestimmung kritischer Belastungen durch Er- 
mittlung der nichttrivialen Gleichgewichtslagen) und das energetische Kriterium (fiir die 
kritische Belastung ist die potentielle Energie erstmals nicht mehr positiv definit) nur bei 
konservativen Systemen verwendet werden diirfen®. Die von A. Trésch® durchgefihrten 
Schwingungsrechnungen ergeben, daB selbst bei Einfihrung von Dampfungskraften die er- 
wahnten unbefriedigenden Resultate nur zum Teil korrigiert werden. 

Indem H. Ziegler * auf die den Momentvektor erzeugenden Einzelkrafte zuriickging, stellte 
er fest, daB sich fir die praktisch wichtigsten Belastungsarten konservative Momente ergeben 
(vgl. auch Abschnitt 2 b), wahrend Kriafteverteilungen, die axial bleibende Momente liefern, 
wohl ziemlich selten vorkommen. Mit den von ihm herausgegriffenen Belastungen definierte 
er das semi-, das quasi- und das pseudotangentiale Moment und zeigte, da mit diesen Mo- 
menten beim Stab mit zwei gleichen Biegesteifigkeiten fiir alle fiinf Arten der Lagerung in 
Abb. 4 plausible Knicklasten gefunden werden * 9, 

Die folgende Arbeit stellt eine Erweiterung der Rechnungen von H. Ziegler dar. Unter 
allgemeinen Voraussetzungen werden im Abschnitt 2 die Ausgangsgleichungen abgeleitet. 
Eine Lésung ohne Einschrankungen fiihrt indessen auf sehr uniibersichtliche Knickgleichungen. 
Wir lassen daher im Abschnitt 3 einen Stab mit zwei verschiedenen Biegesteifigkeiten zu, der 
zwar sogar eine spannungsfreie Verwindung aufweisen kann, jedoch vor dem Knicken nur auf 
Torsion beansprucht sein soll. Im Abschnitt 4 wird dann bei einem Stab mit zwei gleichen 
Biegesteifigkeiten auch eine Druckkraft mitberiicksichtigt. Insbesondere die Knickgleichungen 
dieses Abschnittes werden eingehend diskutiert und die kritischen Dynamen in Diagrammen 
wiedergegeben. 


* Von der ETH in Ziirich genehmigte Promotionsarbeit. 

1]7,, Euler, Histoire de \’Académie, Bd. 13, Berlin 1757. 

2 A. G. Greenhill, Proc. Inst. Mech. Eng. 1883, S. 182. 

3 R. Grammel, Z. angew. Math. Mech. 3, (1923), S. 262; siehe auch C. B. Biezeno u. R. Grammel, 
Technische Dynamik, Bd. 1, 5. 607, 2. Aufl. 1953. 

4 H. Ziegler, Z. angew. Math, Physik 2, (1951) S, 265, H. Ziegler stellte ferner fest, daf fiir ver- 
schiedene Arten der Lagerung jede Drehzalil kritisch ist, sobald die Welle ein beliebig kleines axiales 
Moment zu iibertragen hat, 

5 H. Ziegler, Ing.-Arch. 20 (1952), S. 49 und H. Ziegler, Z. angew. Math. Physik 4 (1953), S. 90 und 
168. Die erste Arbeit demonstriert das Versagen des statischen und des energetischen Stabilitats- 
kriteriums bei einem einfachen nichtkonservativen Stabilitatsproblem, wahrend die zweite eine 
grunds&tzliche Analyse der linearen Stabilitatsprobleme und der Methoden zu ihrer Losung darstellt. 

6 A, Trésch, Ing.-Arch. 20 (1952), 5. 258. 

7 H. Ziegler, Z. angew. Math. Physik 3 (1952), S. 96. 

8 Auch die Berechnung des Einflusses der Torsion auf die kritischen Drehzablen ergibt fiir diese 
drei Momente keine Uberraschungen mehr. H. Ziegler, Ing.-Arch. 20 (1952), 5. 377. ‘: 

9 Einen weiteren, mit dem pseudotangentialen Moment verwandten Fall, die Ubertragung des 
Momentes durch ein Kreuzgelenk, behandelt W. T. Koiter in einer noch nicht veroffentlichten Arbeit. 
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2. Grundlagen. a) Die Differentialgleichungen. Die kritische Belastung 
eines schlanken, geraden Stabes, der durch eine axiale Druckkraft P und durch ein vor dem | 


Knicken axiales, konservatives Moment W (vgl. Abschnitt 2b) beansprucht ist, wird nach 
dem statischen Stabilitatskriterium dann erreicht, wenn zum erstenmal nicht nur fur die | 
urspringliche gerade, sondern auch fiir eine ausgebogene elastische Linie Gleichgewicht még- 
lich ist. 
Fir die Ableitung der Differentialgleichungen der elastischen Linie fihren wir zwei Ko- | 
ordinatensysteme cin (Abb. 1 und 2): das jedem Stabquerschnitt zugeordnete Hauptachsen- | 
system (&,7,¢) und ein ortsfestes System (x, y,2), dessen | 

ee z-Achse parallel zur urspriinglichen Stabrichtung sein soll, 
wahrend sein Ursprung sowie die x- und die y-Richtung in 
jedem Knickfall geeignet gewahlt werden. Die Stellung der 
beiden Koordinatensyste- | 
me zueinander beschreiben 
wir mit den drei Eulerschen | 
Winkeln 0, pg, p (Abb. 3). 
Der Stab besitze die 
Lange J; sein Gewicht | 
werde vernachlassigt. Die 
beiden Biegesteifigkeiten « 
und fs beziiglich der -, baw. 
der 7-Achse werden langs 
des ganzen Stabes konstant | 
vorausgesetzt, ebenso die | 
Torsionssteifigkeit y und 
der spezifische Verwin- 
dungswinkel mw) des un- 
belasteten Stabes. Die |} 
GréBen x2, %,, %¢ beschrei- 
ben die elastische Defor- | 
mation: x, ist die Krim- 
mung in der Hauptebene 
(yn, 6), %, diejenige in der 
(€,¢)-Ebene, und x, stellt | 
den spezifischen Torsions- 
winkel dar. Die einzelnen 
Stabquerschnitte charak- 
Abb. 2. Schnitt durch den knickenden Stab. terisieren wir durch die 
vom Ursprung des (x, y, z)- 

Systems aus gemessene Bogenlange t der elastischen Linie; 
x’, y’ usw. sind Ableitungen nach t. 


Im Stabende B (t = b) greifen ein Moment W (W,., Wy, W.) 


Abb. 1. Der unter einer Druckkraft und 5 % ees | 
einem Torsionsmoment ausknickende Stab. und. eine Kraft K Cie. Kos oa P) an (Abb. Ly Dabei sind 
K, und Ky, Reaktionen; fir gewisse Lagerungen sind es auch 


W, und Wy. Die Dyname (K, W) ergibt im Querschnitt F (Abb. 2) eine Beanspruchung, die | 
durch die Kraftschraube (R, M) wiedergegeben wird. Bei einem schlanken Stab haben wir 


nur das Moment M zu beriicksichtigen, das im (x, y, z)-System die Komponenten (vgl. Abb. 1) 
M, = Wz — Ky [2(b) — 2(t)] — P [y(b) — y(t] , | 
M, = W, + Kz [z(b) — 2(t)] + P[x(b) — x(t)], (2.1), 
M, = W, — Kz [y(b) —y()] + Ky[+(0) — x(0)] | 


besitazt. Fir die Umrechnung ins Hauptachsensystem gelten die mit Hilfe von Abb. 3 ableit- 
baren bekannten Transformationen 


a 
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M; = M;g (cos p cos y — sin y sin y cos 3) + My, (sin y cos y + cos p sin ¢ cos #) 
+ M, sin sin 0, 


M, = M,(— cosy sing — siny cos 9 cos 0) + My (— sin yp sing + cos p cos ¢ cos #) 


-+ M, cos y sin #0, (2) 
M, = M, siny sin 0 — My, cos y sin # 
+ M,cos?. 
Der Zusammenhang 
Mz =oxg, M,=Bu,, Mp =y xe (2.3) 


zwischen der Beanspruchung M und der Deformation liefert uns die Differentialgleichungen 
fiix die elastische Linie und die Torsion im (&, 7, €)-System. 


. Wie sich durch Energiebetrachtungen begriinden lat, diirfen wir uns auf eine Theorie 
erster Ordnung beschranken. Kleine GréBen sind x,y <1, x’, 7',9 <1, He, Hy K jl, ferner 


eM Wn, Wy < Lia und K,, Ky < P, W/l, so daB sich, wenn wir mit y =q@ + yp den 
Winkel zwischen der x- und der &-Achse bezeichnen, die Gleichungen (2.1) und (2.2) auf 


DS (bY 8) — F f(b) — y(t) 


2a 
M, = W, + K, (6 —t) + P[x(b) —x«@)], Mt = W,= oT ( ) 
M; = M, cosy + M,siny + W@sing, ne 

M, = — M,siny + M,cosy + W#écosq, ate Cy 


reduzieren. 
Um (2.3) auf geeignete Variabeln transformieren zu kénnen, driicken wir x;, %,,%; durch 
die Winkel 3, g, y und ihre Ableitungen aus. Durchlaufen wir gleichférmig die elastische Linie, 
so bedeuten x:, x, %- + Wy die Kom- 
ponenten der Winkelgeschwindig- 
keit des Hauptachsensystems. Die 
gesuchten Beziehungen sind also 
die Transformationsgleichungen 
x =y' singsind + 0 cos, 
%,=y' cospsin} — # sing, 
He + Wy =y cosh +9’ 
(sin ® ~ 3, cos # ~ 1) 


der Winkelgeschwindigkeit der 
Kreiseltheorie. Die dritte Gleichung 
von (2.3) lautet nun W =y(x' — wp) 
und fihrt mit w* = W/y auf 
x = (w* + @))t+6=wt+ 6, (2.5) 
wobei 6 eine Integrationskonstante 
ist. 

Die einfachste Darstellung der 
folgenden Rechnungen erhalten wir 
mit den Variabeln 


(2.4) 


f(t) a a cos p ? g(t) Fe Os sin P ” Abb. 3. Die Koordinatensysteme (, y, 2), (€ 7, ¢) und (u, v, 2). 

Wt) = Pp +o. (2.6) e 
Geometrisch bedeuten f, g die negativen Komponenten des Tangenteneinheitsvektors v der 
elastischen Linie im aufgerichteten Koordinatensystem (u, v, 2) von Abb. 3 af = — Vy, = — Vy 


Mit den Komponenten x’, y’ von vim (x, y, z)-System sind f, g somit durch die Beziehungen 


x =fsin(wt + 6) —geos(wt+ 6), y' = — feos(wt + 6) —gsin(wt +4), (2.7) 
f = x' sin (wt + 6) —y' cos(mt + 6), g = — x’ cos (wt + 6) — y’ sin(wt + 4) (2.8) 
verbunden. 


Los 
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Aus (2.4) gewinnen wir mit den neuen Variabeln die Ausdriicke x; =f’ +08, %, = —& 

+ wf, die wir zusammen mit M;, M, aus (2.2a) in (2.3) einsetzen. Wir erhalten die Differen- 
tialgleichungen 

af’ + («0 — W)g = M, cos (wt + 6) + My sin (wt + 0), (2.9) 

Bg’ —(Bw — W)f = M, sin (wt + 6) — My, cos (wt + 9), : 
welche sich in dieser Form fiir die Behandlung des nur durch ein Torsionsmoment hbelasteten 
Stabes cignen. Die Druckkraft fiihrt indessen nach (2.la) auch x, y und damit die Integrale 
von f, g in (2.9) ein, so da® diese Beziehungen [unter Beriicksichtigung von (2.7)] nach ¢ abge- 
leitet werden miissen und wir beim gedriickten und tordierten Steb von den Differential- 
gleichungen ! 
of’ +[(a +8)o—W]e’ +(Ww— fo? + P)f=K, cos(wt-+ 6) — K,sin(wt + 9), 
Be” —[(a+f)wo— Wf’ +(Wo—aw?+ P)g = K,sin (wt + 6) + K, cos (wt + 6) 
auszugehen haben. Die weiterhin giiltigen Gleichungen (2.9) ergeben die Randbedingungen in 
den nicht eingespannten Stabenden. 


(2.10) 


Die kritische Belastung ist die klein- 
QW, WG ste, fiir welche die Gleichungen (2.9) 
Zz Li ee KL bzw. (2.10) eine die Randbedingungen 
erfiillende nichttriviale Lésung  be- 
sitzen. Wir untersuchen die folgenden 
finf Arten der Lagerung (vgl. Abb. 4): 
I. den beidseitig eingespannten 
Stab, 


II. den beidseitig eingespannten 
Stab mit gegeneinander verschieb- 


baren und um die urspriingliche Stab- 


achse verdrehbaren Einspannungen, 
III. den einseitig eingespannten 
Stab, 
IV. den auf der einen Seite ein- 


V/A gespannten, auf der andern gelenkig 
7 ae ie a 
V. den beidseitig gelenkig gelager- 
I I Il WV W 


ten Stab. 


Abb. 4. Die untersuchten Lagerungen des Stabes. b) KonservativeTorsi ons- 
momente. Das belastende Moment 

(Abb. 1) ist dann konservativ, wenn seine virtuelle Arbeit als Variation eines Potentials dar- 
gestellt werden kann. Fiir die virtuelle Arbeit 6Ajy der Momente in den beiden Stabenden 


B und A, Wea W,,,W. = + W) und W 4 (Wax, Way, Woz = — W), erhalten wir 
b 
bAy = W{ du, dt — (Wyx dys — Wey 5x5) — (Wax Oy, — Way dxi) (2.11) 


(dx, ist die z-Komponente des vektoriellen Winkels, um den der Stab pro Langeneinheit 
virtuell gedreht wird). Gleichung (2.11) 1a8t sich auf? 


b 
1 2 Mt Ul aa / 7 
b4y = sa” [= yl! —y! x +20) a + (Ws — 2M. ys] 8x4 — (Waa Wo 35) oy] 


1 fA / il / / 
ae ( V5 5 We, ya) ono té —+W,, ve) bye (2.12) 


umformen. Auch den Randtermen muB ein Potential zugeordnet werden kénnen. Damit das 


* Die Gleichungen (2.10) entsprechen den Beziehungen (6.1) bei H. Ziegler, Z. angew. Math. — 


Physik 2 (1951), S. 274. 
2 H. Ziegler, Z. angew. Math. Physik 3 (1952), S. 105. 
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méglich wird, ist an beiden Stabenden die Bedingung 
fa) 1 , — 7) 1 / 
ee ay a2 —W,) (2.13) 


zu erfillen. Es sind W,, W, kleine GréBen; um sie als Funktionen von x’, y’ darzustellen, 
geniigt deshalb ein Ansatz 


WoW, (Ata -- B*y')5 WW, (Gx = D*y'). 
Die Bedingung (2.13) verlangt, daB 
A* + D* = 1 (2.14) 


semitangentiales Moment: W- w(42;1) W- wo, fy;t) 
WW (Zend ys1) 


Abb. 5a. Das semitangentiale Moment. 


guasitangentiales Moment: W=W (2; 01) W-W(0,0,1) 
, > W_W(2}0,1) 


Abb. 5b. Das quasitangentiale Moment. 


sei. Bei einem konservativen Moment treten also beim Ausknicken die Komponenten 


W,,=W,(A* x + Bry’), Wy = W,[C*x’ + (1 — A*) y'] (2.15) 
1 
auf. Die einfachsten konservativen Momente erhalt man fir A* = D* = aD Beene 0) 


W,,= +, bo V7 = = W,y' (semitangentiales Moment) 
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und fur A™ = 1 Bee Da 
Wo=W x’, Wy=0 (quasitangentiales Moment). 


Diese beiden Momente lassen sich, wie aus Abb. 5 hervorgeht, durch einfache Kraftever- 
teilungen verwirklichen 1. Das semitangentiale Moment gehért zu einer gleichmaBigen Ver- 
teilung von drei oder mehr Kraften von gleichem, konstantem Betrag und konstanter Rich- 
tung, die an einer zum Stab normalen, starr mit ihm verbundenen Scheibe angreifen. Ein 
quasitangentiales Moment entsteht, wenn nur zwei solche Krafte auf die Scheibe wirken. 
Greift dieses Kraftepaar an einem rechtwinklig und starr mit dem Stab verbundenen Quer- 
balken an, liegt ein pseudotangentiales Moment vor mit A*= 1, B* = tg r und C*= D* = 0: 


W,,=W, (x’ +tgt-y')=Wo, cost (x +tgt-y’); Wp 0 (Wop 1a 


Dabei ist 7 der Drehwinkel des Querstabes. Beim pseudotangentialen Moment sind zwei Falle 
zu unterscheiden: ist der Stab im zweiten Ende starr eingespannt (a), so ist t= (W/y) l= w* 1: 


pseuodotangentiales Moment: W=2Pa(z'. cost, 0,c0st) = Pa ( 'y'sint, 0,008 T) 


W=2Pa(zx'cos t+y’'sinz, 0,008 T) 
Abb. 5c. Das pseudotangentiale Moment. 


kann er sich in jenem Ende axial drehen (b), so verteilt sich die durch W bewirkte Torsion: 
t = W1/2 y =w* 1/2. Dat hichstens den Wert 7/2 erreichen kann, ist auch W beschrankt: 
W<ay/2lim Fall (a), W Sa y/lim Fall (b). 

Weder das axiale Moment (A* = B* = C* = D* = 0) noch das tangentiale (A* = D* — 1, 
B* = C* = 0) erfillen die Bedingung (2.14); beide sind deshalb nichtkonservativ. 

Bei den Lagerungen I und II muf die Stabtangente in beiden Enden ihre Richtung bei- 
behalten: dx, = dy, = 6x, = Oy, =90, und die Randterme in (2.12) fallen weg. Das 


belastende Moment ist somit immer konservativ. 


3. Knickung gerader Stibe durch konservative Torsion bei zwei verschiedenen Biege- 
steifigkeiten und spannungsloser Verwindung. a) Allgemeine Lésung der Diffe- 
rentialgleichungen. Die Bestimmungsgleichung fiir das kritische Torsionsmoment 
geht beim Fehlen einer Druckkraft aus den Differentialgleichungen (2.9) hervor. Diese lauten, 
wenn wir (2.la) mit P = 0 einsetzen: 


af’ + (aw — W)g =(W, — Ky b) cos (wt + 6) + (Wy + K, b) sin (wt 4 6) | 
+ K,tcos (wt + 6) — K,tsin(wt + 0), | 
Bg — (Bw — W)f = (W. — K,b) sin (wt + 6) — (Wy + Kyb) cos (wt + d) | “_ 
+ Kytsin (wt + 6) + K,tcos (wt + 0). 
Die homogenen Differentialgleichungen 


of'+(ao—W)g=0, pe’ -(Bw—W)f=0 


1 H. Ziegler, Z. angew. Math. Physik 3 (1952), S. 99. 
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lésen wir mit dem Ansatz 
fn = AycosAt + BysinAt, g,=Cycosdt + Dy sindt. 


Dieser fiihrt auf die charakteristische Gleichung mit der Wurzel 


4 =| o—W) (Bo —W) 


xB (3.2) 
und auf die Beziehung zwischen den Koeffizienten 
ar 
© =—qB,, Dy = q Ay [ee aA 


[A ist immer positiv oder positiv imaginar; um das Vorzeichen von q klar festzulegen, soll der 
Ausdruck dafiir nicht durch Einsetzen von (3.2) vereinfacht werden.] 
Die partikulare Lisung ist die Summe von zwei partikularen Integralen 


fi = A, cos (wt + 6) + By sin(wt + 6), g, = C,cos (wt + 6) + D, sin(wt + 4), 
fe = A, cos (wt + 6) + B, sin (wt + 6) + Etcos(wt + 6) + Ftsin(wt 4 6), 
8 = C, cos (wt + 6) + D, sin (wt +- 6) + Gtcos (wt + 6) + Htsin (wt + 4), 


von denen das erste die Stérglieder ohne den Faktor t, das zweite diejenigen mit diesem Faktor 
bericksichtigt. Fiir die einzelnen Koeffizienten ergibt die Rechnung: 


Wy + b Kx 
Ww 9 
B (20 @—W) Ky 
[(a + 8)» —W]W?’ 
a(2Bo—W) Kz 


CG =latAo WW?’ Nes 


pee fe ee aE 


B(200—W) Kx 
PCr ace 

a(2po—W) Ky CS 
[Ca Po 


An = 


— __ Ky 
EH, Ee 


Beim Zusammenfiigen der einzelnen Bestandteile f, + f, + fo, gn + 8: + 8 der allgemeinen 
Lésung von (3.1) fiihren wir noch die Abkiirzungen 


a — fp 
Ae A, By by te GEap or 


ein, beriicksichtigen ferner die Beziehungen 
C,+C,= —B,—B,—nE, D+D,=A4,+4,—nF 
und erhalten damit 
f(t) = A, cosAt + By sind t + Aj, cos (wt + 6) + By sin (wt + 6) + Et cos (wt + 8) ] 
+. Ftsin(wt-+ 6d), | 
g(t)= — q By cosAt +q Ay sinAt — (By +E) cos (wt + 6) + (Ay —nF) sin (wt + 0) 
— Ftcos(wt + 6) + Etsin(wt + 6). 


Wir suchen nun die zu den Lagerungen I bis V gehérenden spezicllen Lisungen, beschranken 
uns jedoch in den Fallen IV und V auf das semitangentiale Moment. 


b) Die Knickfalle II und II. Da sich das Lager B gegentiber A in der x- und 


(3.4) 


y-Richtung frei verschieben kann, tritt keine Lagerkraft K auf: K, = K,=0. Die Aus- 
driicke (3.4) fiir f und g vereinfachen sich deshalb wesentlich: 


f(t) = Ay cosAt + Bysindt aes: (wt + 6) + 8 sin (wt +6), 
(3.5) 


; Wx Wy . 
g(t) = —q Bycosdt + q A, sind t — = cos (wt + 9) — yy sin (wt + 6). 
In beiden Fallen legen wir den Ursprung des (x, y, z)-Systems in das eingespannte Stabende Fs 
Der Neigungswinkel # der Stabtangente in A (Abb. 3) hat dauernd den Wert Null; fiirt= a= 0 
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sind somit [vgl. (2.6)] die Randbedingungen 


f(0) = Ay — 2 cos 6 + Tsin 5 =0, | 


We Wy ee 
g(0) = — qB, — 7 008 0 — Fein b= 0] 


zu erfillen. 


In den nicht eingespannten Stabenden ergeben sich die Randbedingungen daraus, da die | 


GroBen W,, W,, die zwar nach (3.3), baw. (3.5) den Verlauf der elastischen Linie bestimmen, 


ihrerseits nach Abschnitt 2b durch die Tangentenrichtung in diesen Stabenden festgelegt sind. _ 


a) Fall Ill, semitangentiales Moment. Das Moment in B schreibt keine | 
x-Richtung vor; wir kénnen diese mit der -Richtung im Einspannquerschnitt A zusammen- | 
fallen lassen und erreichen damit das Verschwinden der Integrationskonstanten 6. Aus (3.5) | 


und (3.6) gehen dann die Gleichungen 


SUZ Wena alie } 
flO) = a (cosAt — cos wt) ee oe: + sin wt) a | 
3.7) 
i aed eee Pa 
g(t) => (qsin At — sin wt) + > (cost — cos w t) | 
heryor. Wir formulieren sie speziell fiir t =1 und stellen diese Ausdriicke fiir f(l), g(l) den- 
jenigen gegeniiber, die bei einem Moment mit W,,= oe «x (l), Wy= = y' (I) aus (2.8) gewonnen | 
werden: 
pe cae Wy 
ie 2-7 sin wl — 2 =; cos wl, 
i‘: We We 
a(l) = —2 ep cos wl — 257 sin wl. 


Auf diese Weise entstehen die beiden homogenen, linearen Gleichungen 


Ws 
W 


(— sina — sin ol) + 72 (cool +coswl) =0, 
7 (cos Al + cos wl) +42 (qsinal +sinwl) =0 
fir W,, W, (und damit fiir By, A,). Die Bedingung fiir nichttriviale Lésungen 
242 coshleosml +(q + 7)sindlsin wl =0 (3.8) 


(III, semitangentiales Moment) 


gibt implizit den Wert des kritischen Moments an. Die Knickgleichung ]a8t sich in ein Produkt 
aufspalten: 


; l l l l 
(tg 5 +ateo-5)(qergd 5 + tg >) <0. (3.8a) 


Mit dem Ursprung des (x, y, z)-Systems in der Stabmitte fiihren die Rechnungen direkt auf | 


diese Form der Knickgleichung. Die Knicklast wird durch die scharfere der beiden in (3.8a) 
enthaltenen Forderungen bestimmt. 


6) Fall Ill, quasitangentiales Moment. Wenn wir das Moment in B dar- 
stellen als W = W (x’ (1), 0, 1), ist nach Abb. 5b die x-Achse festgelegt. Der Winkel 6 ist nicht 


mehr Null, hingegen vereinfacht sich (3.6) durch Wy = 0. Statt (3.7) erhalten wir 
Wx |. 1 
= Sai, sin 0 cos At + ee Osindt — sin (wt + . 


, 
Wx G2) 
g(t) = =F [cos 6 cosAt — qsin 6 sind t — cos (wt + 6)]. | 


(3.6) | 
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Die mit (2.7) gebildete Beziehung 


= = x'(l) = f(l) sin (wl + 6) — g(l) cos (wl + 8) 


wird zur Knickgleichung, sobald darin fir f(I), g(l) die Ausdriicke (3.9) eingesetzt sind: 


cosAlcoswl-+sinAl =; £08 6 sin (wl + 6) — qsin 6 cos (w1+ 6)| = 0. (3.10) 
(III, quasitangentiales Moment) 


y) Fall OI, pseudotangentiales Moment. Wahlen wir die x-Richtung 60, 
wie es Abb. 5c verlangt, so gelten auch hier die Gleichungen (3.9). Der Momentvektor hat im 
Fall (a), d.h. beim in A starr eingespannten Stab (vgl. Abschnitt 2b), die Komponenten 
W, = W(x'(l) + y'() tg w* l), Wy = 0; aus (2.7), (2.5) folgt 


TF cos w* 1 = f(l) sin (ivy 1 + 8) — g(t) cos (v1 + 6), 
und nach der Substitution von (3.9) fiir f(I), g(l) entsteht daraus die Knickgleichung 
cosAl cos wl + sindl te cos 6 sin (wy! + 6) — qsin 6 cos (a) 1 + 0) =O u( anv ka) 
[1II, pseudotangentiales Moment (a)] 


Das gleiche Vorgehen ergibt im Fall (b) (Stab in A drehbar) 


cos Al cos [iS + @) + sind ae cos 6 sin ee + a)! — 6| 


— qsin 6 cos le + 0) (l + | =i 0. (3.11b) 
[1II, pseudotangentiales Moment (b)] 


6) Fall II. Die Komponenten W,, W, des Moments in B sind Reaktionen. Wie im 
Fall III (semitangential) kénnen wir die x-Richtung so wahlen, da die Beziehungen (3.7) 
gelten. Im eingespannten Stabende B miissen die Forderungen f(l) = (C, g(l) = 0 erfillt sein. 
Sie verwandeln (3.7) in zwei homogene Gleichungen fiir W,, W,, die auf die Knickbedingung 


2 —2eosAl cos ol — (4 +7) sindlsin ol = 0 (II) (3.12) 


oder 


(4 tga “5 — tg05)(tea-z —qtg 05) ='() 
fiihren. 
Hier wie auch im Knickfall I, ist es fiir die GréBe des kritischen Moments unwichtig, mit 
welcher Krafteanordnung es erzeugt wird. 
c)Die Knickfalle I, IV und V. Bei den Lagerungen J, IV und V sorgt die Lager- 


kraft (Abb. 1) dafiir, daB das Stabende B sich nicht senkrecht zur urspriinglichen Stabachse 
verschiebt. A und B miissen also dieselben x- und y-Koordinaten aufweisen: 


; +1/2 , +12 
ol CE a 
, —1/2 —1/2 * 


(Die Stabmitte sei der Ursprung des ortsfesten Koordinatensystems.) Die Integranden sind 
mit Hilfe von (2.7) auf die Variabeln f, g zu transformieren. Falls bei TV und V nur semitangen- 
tiale Momente beriicksichtigt werden, darf in den Ausdricken (3.4) fur f, g die Integrations- 
konstante 6 weggelassen werden. Bei der Integration von (3.13) fahren wir noch die Ab- 
kiirzungen 
fee a (2B e—W) me B(2a@—W) 
[(a + B)o—W]WIl ” [(« + B)o—W]Wl 
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ein und erhalten 


sin @— cos @ = 
abe re! 2 | es 
By — 2rsinA— cos w-— + 2sq cosd > sin w > + By+En oe =o 


wl 


j ’ is sin —- COB @— 
A Nee er i baa, re a hol ae 
(3.13a) 


a) Fall I. Die Lagerungen in A und B sind Einspannungen, in denen f, g den Wert Null 
beibehalten. Wir kénnen somit die Beziehungen 


Hs) -3)=0 af) a) 6) A) 208 af) a4) 


aufstellen, die zusammen mit (3.13a) je drei homogene, lineare Gleichungen fiir Aj, A,,, F und 
By, By, E bilden. Eine der beiden Koeffizientendeterminanten mu8 beim Ausknicken Null | 
sein, entweder 


‘ , 1 sin ®— cos ers 
— 2sqsini— cos w—> + 21 cosd—sinw > 1 Leos an ol | 
2! l eee =! | 
cos =e cos W > 5} 5) 
ind y i l s n sin W : 
q sin sin W > gq C08 W > o> 
oder (3.14) 


: l 
Sates! I eon eur eee eee 
— 2rsind = cosw> + 2s qcosd —sin w > 1 n\—+ 


2 2; wl 
; l : l l 1 =—(0. 
sin A — sin 0) —> yz 008 W-> 
ae Ba ae a Bye Wie 
end 2 2 2 


Welche der beiden Forderungen das kleinere Knickmoment liefert, hangt vom Verhaltnis von 
a, Pp und y, sowie von Wy ab. 


B) Fall IV, semitangentiales Moment. Die Randbedingungen in A stim- 
men mit denen des Knickfalls I tiberein: f(— 1/2) = 0, g(— 1/2) = 0. Im Gelenklager Blassen | 
sich W,/W und W,/W auf zwei Arten durch die zu bestimmenden Koeffizienten der allge- | 
meinen Lésung (3.4) ausdriicken: erstens indem wir (3.3) mit b = 1/2 und unter Verwendung 
der eingefiihrten Abkiirzungen nach jenen GréSen auflésen, und zweitens indem wir beriick- 
sichtigen, da} W, und W, die Komponenten eines semitangentialen Momentes sind 


Ws l at ei) 

= po Bat Fa tEsl=>e(5), 
W. l Le fi Gy) 
i ee = aa tig ae PTH ee 
dey eth Ep Bal yi 


[rechts setze man zunachst (2.7) dann (3.4) ein]. Nach Erginzung der vier fiir die Stabenden 
A und B aufgestellten Forderungen durch (3.13a) wird die Knickgleichung durch Nullsetzen 
der Determinante dieses Gleichungssystems gebildet. Diese lat sich zwar auf (3.16) ver- 
einfachen aber nicht in ein Produkt aufspalten, da die Randbedingungen nicht sym- 
metrisch sind [die linke Seite von (3.16) ist als 6-reihige Determinante zu lesen]: 
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9 ond l l 1 1 1 BUR COS, 
— Ss sin A — = SS ot 
q 5 COR W + Ar cosA-— sin w 1 Se <j 
roves 0 Seer 
2 2 
agian ae 0 sr lieta oe 
2 9 9 
0 0 0 
0. cos ee John ! ! l 
5 5 w= — sleosw— 
A l l l l 
0 = es ee ee i ae 
sin Ww > 7 C08 W _ slsin w 5 
0 0 0 
0 — sin uf ea ae l I si : 
9 5) OF Ss ee cy 
0 on oe eee 
az ees Pe onan a 
9 5 Os sloosw 
1 1 1 (: sin ®— cos “7 =0. 
= YiPs =e sili aa eZ 1 
rsin 4 cos w-> + 2sq cosh sin w = n soe a 
ile 0 isla 
DP 2) 
= nee 0 ay po aines 
q 5) 2 
(IV, semitangentiales Moment) (3516) 


y) Fall V, zwei semitangentiale Momente. Die Beziehungen (3.13a) 
und (3.15) kénnen als Randbedingungen sofort ibernommen werden. Die (3.15) entsprechen- 
den Gleichungen fiir das Gelenklager A, in dem das Moment ebenfalls semitangential einge- 
fihrt wird, lauten 


u{—) Wa Shy l 1 I 
ae a ay ee Est=tv(—4), 
(3.17) 
m,(—3) en, Hee 1 1 I 
W : a 7 ri sare mie cee oe Ay 
Diese Randbedingungen fihren auf die beiden Knickgleichungen 
3 l 
sin ®— Cos @W— 
: l l ope: l 1 2 2 
sin 1 = cos Ww + 2r cos sinw > 1 _.(4 7 
l eee l l == ()), 
cosh — ©08 W > — 7 sinw= +2sleosw > : 
l ; oe! l . l 
q sind > — sin W—> cos w= + (281 — n) sin wo 
(3.18) 
sin ® — cos w@— 
Pl eas fale le ky ne ) | 
nA cos WT qs cosA- sinw > a PY 
sin A —sinw 5 — cos w+ — 2sIsin wo = 0. 
yee Posie Ge (251 =n) cos aye 
q 2 TaD 2 2 


(V, zwei semitangentiale Momente) 
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d) Die Auswertung der Knickgleichungen, Um fir eine vorgegebene | 
Welle das Knickmoment zu ermitteln, hat man in die zugehérige Knickgleichung Q (W, a/s 
Bly, Wp», 6) zunachst die Zahlenwerte fiir «/y, B/y, Wo, 6 einzusetzen und darauf Q mit Ver- | 
suchswerten fiir W solange abzutasten, bis man der Nullstelle mit dem kleinsten W genigend | 
nahe ist. Obschon einzelne Knickgleichungen nicht leicht zu tberblicken sind, gelangt man 
auch dort mit einer Handrechenmaschine relativ rasch zum Ziel. Die drei folgenden Hinweise 
werden das Aufsuchen der Knickmomente in vielen Fallen erleichtern. 


a) Pseudotangentiale Momente. Fir den durch ein pseudotangentiales Mo- | 
ment belasteten Stab besteht im Fall (a) keine Knickgefahr. Mit w) = 0, 6 = 0, 2/2 1a8t sich | 
die Knickgleichung (3.1la) der Lagerung III 


cos = jot) =a () 
ap 


unter den allgemein giiltigen Bedingungen («4 + 8) >y,20 >y, 26 > y mit |w*| U < 7/2 nie 
lésen. Auch fiir die ibrigen Werte von 6 und fiir w) ~ 0 wird, wie eine ausfihrliche Unter- | 
suchung gezeigt hat, mit dem gréBtméiglichen Moment die Knicklast nicht erreicht. Da die 
kritische Belastung eines konservativen Systems beim Hinzufiigen von Bindungen, die keine 
Arbeit leisten, auf keinen Fall abnimmt 1, gilt dieses Resultat auch in den Knickfallen I, II, 1V, | 
wahrend bei der Lagerung V der Fall (a) kaum in Frage kommt. Hingegen darf im Fall (b) 
nach (3.11b) die Knickgefahr nur bei nichtschlanken Querschnitten mit zwei annahernd 
gleichen Biegesteifigkeiten unberiicksichtigt bleiben. Zwar laBt sich auch ein Stab von schlan- 
kem Querschnitt (« >y >f) mit w) =0 nicht knicken, wenn die das pseudotangentiale 
Moment erzeugenden Krafte (Abb. 5c) in der urspriinglichen Langsrichtung des Querschnittes 
zeigen (0 = 0). Fir 6 = 7/2 jedoch, mit 


2) © Bie eee 
4, =) ene |w*| und qi ACen | 
lautet die Knickgleichung (3.11b) 
Coj A; 1 cos _ = qi Gin Al sin 2 —0 


und weist eine Lésung |w*| 1 < a auf. 


6) Nichtschlanke Querschnitte. Die von R.Grammel bei des Lagerung V 
fiir zwei axiale Momente und fir nichtschlanke Querschnitte (« >/ > y) berechnete Kurven- 
schar ? ]aBt sich, wie Tabelle 1 zeigt, teilweise auch bei unseren Knickgleichungen verwenden. 


Tabelle 1. | 
Knickfall | I | II III, s III, q IV, s V,s—s by aoe 
es aa 35,97) | 25,13 | 12,57 | °6,285 | 27,25 |§19,65 25.13 | 
un 22,76 | 16,02 7,59 | 3,80 17,20 | 12,35 | 16,2 
&O=2B = 4y 
— 1,58 1,57 1.66 fitael 65 1,59 1,59 1,55 


In der ersten Zeile sind die Knickmomente fiir zwei gleiche Biegesteifigkeiten eingetragen, 
die aus der Tabelle 2 (Abschnitt 4) auf die dimensionslose GréBe w, = Wl/y umgerechnet 
wurden. Die Knickmomente der zweiten Zeile wurden aus den Knickgleichungen dieses 
Abschnitts ermittelt und beziehen sich auf einen Stab, dessen beide Biegesteifigkeiten noch 
wesentlich gréRer sind als y. Da das Verhaltnis der beiden kritischen Momente in allen unter- 
suchten Fallen ungefahr denselben Wert aufweist, lat sich bei solchen Staben die gefahrliche 
Last mit Grammels Kurven oder mit der einfachsten unserer Knickformeln gut approximieren. 


* H. Ziegler, Z. angew. Math, Physik 3 (1952), S. 114. 
2 R. Grammel, Z. angew. Math. Mech. 3 (1923), S. 268 oder C. B. Biezeno u. R. Grammel, Technische 
Dynamik, Bd. 1, S. 605, 2. Aufl. 
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y) Verwundene Stabe. Abb.6 demonstriert den Einflu8 der spannungsfreien 
Verwindung w, auf das Knickmoment!. Dieser du®ert sich besonders deutlich bei einem sehr 
schlanken Querschnitt « S y >, fiir den die Knickgleichung im Fall III bei einem semi- 
tangentiaten Moment [vegl. (3.8)] 


2+ 2 cos joo cos 0 22 [ 
\ 
BES gin gosing = 0 20 ag 
eo \ 
78 
lautet. Dabei bedeuten 
6 
Wi 
pe OO 
o= wl + Shae, . 12 
In Abb. 6a finden sich die zur Be- ” 
stimmung der Knicklast maBgeben- —, 
den Kurvendste dieser Funktion 
wiedergegeben. Um daraus die 6 
eigentlichen Knickkurven zu_be- , 
rechnen, hat man sich auf eine 


spezielle Querschnittsform undein 2 
spezielles Material festzulegen: fir 
einen Eisenstab (Querzahl m=10/3) 7 i 
mit sehr schmalem Rechteckquer- _, 


schnitt? (6 = - y) gelten die Kur- 


ven von Abb. 6b, welche, wie die 
Kurven in Abb. 6a, symmetrisch 
sind beziiglich des Koordinaten- -g 
ursprungs; das Vorzeichen von Wp 
bezieht sich auf den Drehsinn des . 

. ‘ : Abb. 6a. Hilfskurven zur Berechnung des Knickmoments des verwundenen 
belastenden Moments. Die Knick- Stabes mit sehr schlankem Querschnitt. 


| 
| 
= Bl ees 


My, Opb+Tdy| 


-16 -4 =12, -10 =6. -6 -4 ~ 0 Z 4 6 8 


Abb. 6b. Das kritische Moment fiir einen verwundenen Hisenstab mit sehr schmalem Rechteckquerschnitt, 


1 Wir setzen hier voraus, daB die Torsionssteifigkeit von w, unabhingig sei. Diese Annahme ist nur 
fiir kleine Werte von w, bewiesen: P. M. Riz, C. R. del’ Académie de VURSS 23 (1939), S.17 u. 141 und 
A. Lourie, G. Janelidze, C. R. de l’Académie de ?URSS 24 (1939), Ss Wh ti, BRA. 

: 2 Vel. z. B. Timoshenko u. Goodier, Theory of Elasticity, 5. 277, New York 1951. 


244 Beck: Knickung gerader Stibe durch Druck und konservative Torsion. Ingenieur-Archi v 


kurve in Abb.6b weist bei w)/ = — 6,8 einen Sprung auf, mit dem sie itber die Asymptote 
Wy! + w, =0 hinitberwechselt. Mit einer ahnlichen, jedoch kleineren Unstetigkeit wird 
1 + w, = —2 itbersprungen. Allgemein besitzt die zur Erzeugung einer ausgebogenen 
elastischen Linie erforderliche Belastung Maxima bei w)l + w, =n27,n = 0, —1,+2,... 


(Ausnahme: n= +1). Die Kurve, die in Abb. 6b die Verwindung vor dem Knicken 
wiedergibt, weist bei diesen Funktionswerten Spitzen auf. Sie korrespondieren mit den | 
Schnittpunkten der je paarweise gegeneinander oszillierenden méglichen Knickkurven, von 
denen jeweilen die der Abszissenachse nichstgelegene das Knickmoment bestimmt. Die 
Ordinate der Knoten nahert sich asymptotisch dem Wert w, = 22 f/y = 4,08. Die kleinste 
Knicklast erhalt man nicht wie beim gedriickten Stab! fiir wy) = 0, sondern fur, ! = 07 | 
Daf die spannungsfreie Verwindung und die Torsion gleichartige GréBen sind, bewirkt fiir das 
kritische Moment eine ausgepragtere Abhangigkeit von Wp als fiir die knickende Druckkraft. 


4. Knickung gerader Stabe mit zwei gleichen Biegesteifigkeiten durch Druck und Torsion. 
a) Allgemeine Lésung der Differentialgleichungen. Die Aus- 
gangsgleichungen (2.10) lassen sich, sobald die beiden Biegesteifigkeiten einander gleich- 
gesetzt werden (x =f), mit K = Ky +i K, zu einer einzigen Differentialgleichung fir die 
Variable F = f + 1 g zusammenzufassen: 


oF’ —i(2a0m— W)F +(Ww —o1w? + P)F=Ke”’. (4.1) 
Wir fiihren gleichzeitig die komplexen GréBen M = M, —i My,W = W, —i Wy,z=x +iy | 


ein. 

Dain jedem Querschnitt alle Geraden durch den Schwerpunkt Hauptachsen sind, brauchen 
wir eine spannungsfreie Verwindung nicht mehr zu bericksichtigen. Bei geeigneter Wahl der 
x-Richtung verschwindet die Integrationskonstante 6 in (2.5). 

Im Falle von zwei gleichen Biegesteifigkeiten kénnen wir ferner zu dimensionslosen GréBen 
ibergehen. Statt mit W und P rechnen wir mit 

_ Wi _ PP 
Cr res etn 


und. bilden entsprechend 
Wel Wyl 
Wy = 9 [Agi = 


_ KP 
20 - ; 


(os 


Ki 


M1 
20 


Im folgenden sind x, y, z,Z,t mit einem dem Stab die Lange 1 zuordnenden MaBstab gemessen, 
w bedeutet nun die ganze Verwindung, und Ableitungen sind nach der dimensionslosen Bogen- 
lange t auszufiihren. 
Die Differentialgleichung (4.1) wird damit auf 
F’ —i2(w —w)F + (2ww —w? + p)F=ke”’ (4.2) 


transformiert. Wir fiigen zur Lisung der zugehérigen homogenen Differentialgleichung 
Fe=Ae+Be# (AA=4,+i14;, B=B, +i B,) 
mit den Wurzeln 
Aye =o — w+ Jw? + p (4.3) 
der charakteristischen Gleichung das partikulare Integral 


F es G BO k APN: 
he eases 

P 
hinzu und erfassen mit 


P=Aeo 2B etl Gee (4.4) 
alle Lésungen von (4.2). 
b) Die Knickfalle IT und II. Die kinematischen Randbedingungen kénnen aus 


Abschnitt 3 direkt tbernommen werden. Bei den Lagerungen III und II tritt keine Lagerkraft 
k auf: k =C = 0. Das eingespannte Stabende A [in beiden Fallen der Ursprung des (x, y, z)- 


1 H. Ziegler, Schweiz. Bauztg. 66 (1948), S. 463 und E, Hui, Knickung verwundener Stabe, Diss. 
ETH, in Vorbereitung. 
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Systems] stellt die Forderung F(0) = 0, die (4.4) auf 
Fe) =A (e'* _ eat) (4.5) 


vereinfacht. Beim Knickfall ITI greife in B ein allgemeines konservatives Moment [vgl. (2.15)] 
an. 


a) Fall Ill, w, = A*wsx', w,=(1—A*)wy’. Die mit Hilfe von (2.la) und fiir 
t = b = 1 formulierten Gleichgewichtsbedingungen (2.9) 


f'(l) + (@ — 2 w) g(1) = 2 w, cos w + 2 wy sin w, | 
g'(1) — (w — 2 w) f(1)= 2 w, sinw — 2 w, cos w | 


lassen sich bei einem allgemeinen Moment nicht komplex zusammenziehen, weshalb wir F 
zerlegen miissen: 


f= —2sin Jw? + pt[A, sin (w — w) t + A; cos (w — w) t], 
g = 2sin Jw? + pt[A, cos (w — w) t — A; sin (w — w) ft]. 


(4.6) 


Indem wir diese Ausdriicke differentiieren, gewinnen wir die Beziehungen 

f' =Vu? + petg Jw? + ptf —(w —w)g, 

g' = (w — w) f + fw + p etg fu + ptg. 
Wir setzen sie in (4.6) ein, ebenso w,, wy mit Hilfe von (2.7) und erhalten 
f(1) [V6 + pctg Jw? + p+ 2w (1—2 A*) sin wcos «| + g(1) w (1 — 2 A*) (sin? w— cos? w) = 0, 
f (1) w (1 — 2 A*) (sin? w — cos? w) 4+ g(1) [ Vw? + pctg Vw? + p— 2w (1 —2 A*) sinw cos co| == O)e 
Die Bedingung fiir die Existenz nichttrivialer Lésungen, 

[Yu + p cos fw? + p— w (1 — 2 A*) sin fw? + p] [w® + p cos Ju? + p 
+ w (1 — 2 A*) sin Jw? + p] = 0 (4.7a) 


liefert zwei mégliche Knickgleichungen. Wir werden im Abschnitt 4 f zeigen, da die kleinere 
Knicklast fiir A* > 1/2 dem zweiten, fiir A* < 1/2 dem ersten Faktor entspringt; d. h. es gilt 


Jw? + pcos fw? + p—w|1 —2 A*|sinfw? + p =0. (4.7) 
DUTT 40, =tA* winx’, 2, = (1 = A*) a0] 
In (4.7) ist mit A* = 1/2 die Knickgleichung fiir ein semitangentiales Moment 
a 2 
cosJw? + p=0, 7 = w? + p (4.8) 


(III, semitangentiales Moment), 


mit A* = 1] die Knickbedingung fiir ein quasitangentiales Moment 


Vw? + pcos Vw? + p —w sin yw? + p=0 (4.9) 
(III, quasitangentiales Moment) 


je als Spezialfall enthalten. , 
6) Fall III, w, =w (*’ + B* y’), wy = 0. Mit den gleichen Uberlegungen wie im Unter- 
abschnitt «) kénnen wir hier die Knickgleichung 
- + ( 5 cos 2¥w? -+ p+ B* w/w? ap sin 2 Jw +p=0 (4.10) 
[III, w, = w(x’ + B*y’), wy =0] 


ableiten. Setzen wir zunachst B* = 0 (quasitangentiales Moment), so erhalten wir eine mit 
(4.9) aquivalente Gleichung, wihrend mit B* = tg aus (4.10) die Bestimmungsgleichung 


z+ (« | cos 2 Vw? + p+tgt:w)w? + psin2/w? + p= 9 (4.11) 
(III, pseudotangentiales Moment) 


fiir die Knickdyname bei einem pseudotangentialen Moment hervorgeht. 
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Beim in A starr eingespannten Stab ist tT = 2 (a); gestattet hingegen das Lager A 
eine Drehung des Stabes, so setze man T = (b). Ferner beachte man die Schranken 


Wimax = = . = im Fall (a) und wax = + . + im Fall (b). 


y) Fall Il. Mit der Bedingung F(1) = Ale © — ae) sorgen wir dafur, daB auch im 

Lager B die Stabtangente parallel zur z-Achse bleibt. Sie fihrt auf 
+(4,—A) = 2)w?+p=—22, =w’ +p (4.12) 
(II). 

c) Der Knickfall I. Wir legen fiir die Lagerungen I und V den Ursprung des (x, y, z)- 
Systems in die Stabmitte, damit sich die Zwischenrechnungen durch die Symmetrie der Rand- 
bedingungen vereinfachen. In allen drei Fallen darf sich A gegeniiber B nicht verschieben, 
d. h. die elastische Linie mufi die Bedingung . 


+1/2 
[za | -iy) a= 0 


1/2 


erfiillen, die wir mit Hilfe von (2.7) sofort integrieren kénnen: 


+1/2 +1/2 exp | le | ] 
dz = iFe('d=A 
A, —©@ 

1/2 ie (4.13) 

exp | ie aera | | 

+B A, — © qn Gnne Wh | 
Zwei weitere Bestimmungsgleichungen fiir A, B und C gewinnen wir aus den Forderungen 
F == ==i0, a = 0, welche den knickenden Stab den Einspannungen A und B an- 

passen. Aus der Liésbarkeitsbedingung dieses Gleichungssystems geht die Knickgleichung 

psinjw? + p +2 |w? + p (cos jw + p — cos w) ==.) (4.14) 


(1) 
hervor, die mit der von A. Trésch1 fiir die Lagerung I abgeleiteten Knickbedingung gleich- 
wertig ist. 

d) Der Knickfall IV. Im Knickfall IV erweist es sich als giinstig, die Bogenlange t 
der elastischen Linie vom Stabende A aus zu messen. An die Stelle von (4.13) tritt dann die 
weiter unten formulierte Beziehung (4.17b), wahrend die Einspannbedingung auf (4.17a) fihrt. 
Die fiir das Stabende B aus (2.9) abzuleitende dynamische Randbedingung wird, im Gegensatz 
zu den beiden ersten, von der Art des belastenden Moments abhangig sein. 

o) Semitangentiales Moment. Wir fassen (2.9) komplex zusammen: 


F(t) —i(w —2w) F(t) =2me, (4.15) 
haben im Stabende B speziellt = 1 und 
m(l) =w = ¥(«(1) —iy(l)) = i= F(1) e” 
zu setzen, und erhalten 
F’(1) —i(w — u) F(1) =0, (4.16) 


woraus durch Substitution von (4.4) die Beziehung (4.17c) gefunden wird. Die Koeffizienten- 
determinante des Gleichungssystems 


A +B +-C =), 

exp t (A, — w) — exp i (A, —w) —1 “ 
ar aye ee sop Pee moter Mio, (4.172, b, ©) 
A jw? + pe” —Byw? + pe” +Cwe” =0 


1 A. Trésch, Ing.-Arch, 20 (1952), S. 265. 
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hat fir die kritische Belastung den Wert Null: 
P Vw? + pcos Jw? + p — (2 w? + p) sin Jw? +p+2 w yw? + p sinw = 0. —(4.18) 
(IV, semitangentiales Moment) 


B) Quasitangentiales Moment. Fir w, = w x’'(1) ,wy=0 kann aus (2.9) nicht 
mehr eine (4.16) entsprechende komplexe Randbedingung gebildet werden, sondern es entstehen 
die beiden reellen Forderungen 

f'(1)— 2 w sin w cos w f(1)+ (w — 2 w sin? w) g(1)= 0, 
8 (1)— (@ — 2 w cos? w) f(1)+ 2 w sin w cos w g(1)= 0. 
Vor der Zerlegung von F(t) in Real- und Imaginarteil kann C mit Hilfe von (4.17a) eliminiert 


werden. Indem wir nun versuchen 4,, 4;, B,, B; aus (4.19a, b) und der zerlegten Gleichung 
(4.17b) zu ermitteln, stoBen wir auf die Knickgleichung 


(4.19a, b) 


Dy 2 ®, iY 
@, sin M, P, cos P, ®, sin ®, — P, cos @, 
— ®,(1—cos ®,) —@, sin ®, + p ®,(1—cos ®,) —@,sin®,+p | ~ pee) 
— M, sin ®, + p ®,(1— cos ®,) —@,sin®,+ p — @, (1— cos ®,) 


(IV, quasitangentiales Moment) 
Es bedeuten 


0, = Jw? + p —w, @,=|wtptw. (4.21) 
Die Zwischenrechnungen vereinfachen sich, wenn man (4.19a, b) durch [(4.19a) sin w — (4.19b) 
cos w| und [(4.19a) cos w + (4.19b) sin w] ersetzt. Ferner beachte man (A, — w) (A, — w) = — p 


und (A) (2) — w) + 2 w = — (dea) — @) . 
y) Pseudotangentiales Moment. Als Randbedingungen in B erhalten wir hier 


aus (2.9) die Beziehungen 


f' (A) —2 w (sin mw — tg Tt cos ~) cos wm f(1)+ [w — 2 w (sin w — tg t cos w) sin w] g(1)= 0, 
} g (1) — [w— 2 w (cos w + tg 7 sin w) cos w] f(1)+ 2 w (cos wm + tg tT sin w) sin w g(1)= 0. G2) 
Im ibrigen ist das Vorgehen bei der Ableitung der Knickgleichung dasselbe wie beim quasi- 
_ tangentialen Moment und fiihrt auf 


—@D 

Say ee ee P, cos D, tattle P, 

— P, (1— cos @,) — @, sin @, + p 

— P, sin @, + p @, (1 — cos ®,) 
| P, P 

@, sin D, + 2wtgt (1—cos ®,) — P, cos P,— 2 w tgt sin ®, 
P, (1 — cos ®,) — , sin ®,+ p sa ee 
— @,sin ®, + p — ®, (1— cos @,) 


(IV, pseudotangentiales Moment) 


e) Der KnickfallV. Bei der beidseitig durch Gelenke oder kurze Lager gehaltenen 
Welle kann die Tangente sowohl im Stabende B als auch in A beim Ausknicken die Richtung 
. andern; wir miissen deshalb fiir beide Enden wissen, wie das belastende Moment eingefiihrt wird. 
Aus der groBen Zahl der zu unterscheidenden Fille greifen wir diejenigen heraus, bei denen die 
_ beiden Momente gleichartig sind; auch verlangen wir, daB, wenn quasi- oder pseudotangentiale 
_ Kraftepaare wirken, die erzeugenden Kriafte in A und B gleichgerichtet seien (Abb. 5b, c). 

Die Randbedingungen in B sind die gleichen wie im Fall IV [vg]. (4.16) beim semitangentialen, 
| (4.19) beim quasitangentialen und (4.22) beim pseudotangentialen Moment], doch ist diesem Stab- 
'ende jetzt die Bogenlinge t = 1/2 zugeordnet. Durchwegs t= 1 durch t= 1/2 ersetzend 
/kénnen wir ihre dortige Formulierung iibernehmen, und mit t = — 1/2 entstehen aus jenen 
} Beziehungen auch die Gleichungen der Randbedingungen in A, da diese ja mit denen in B iiber- 
j einstimmen. Ferner ist (4.13) zu beriicksichtigen. 

| 17 
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Aus diesen Randbedingungen geht fiir semitangentiale Momente in A und B die Knick- 
gleichung 


2 w (cos jw? + p— cos w) +p jw? + p sin |w? +p=0 (4.24) 
(V, zwei semitangentiale Momente) 


hervor. Im Falle von zwei quasitangentialen Momenten 1aBt sich die Knickdeterminante in ein 
Produkt aufspalten. In der Bedingung 


(2 w? + p) sin Jw? + p — 2 w jw? + p sae) — 0 (4.25) 
(V, zwei quasitangentiale Momente) 
stecken die kleineren Knicklasten, so daB die zweite Forderung 
p (2 w? + p) sin Jw? + p + 2p ww +p sin w + 8 w? yw? + p (cos yw? + p— cos w) = 0 


nicht beriicksichtigt werden mu8. Wirken in beiden Stabenden pseudotangentiale Momente, so 
erhalten wir die maBgebende Knickgleichung 


p (2 w? + p) sin /w? + p-+ 2 pw |w? + p sin w+ 8 w? |/w? +p (cos Jw? + p— cos w) 
+4w ig le w? + p) sin Jw? + p—2 wy w? + p sin w Sa +p (cos Jw? + p+ cos w) —0 
(4.26) 


(V, zwei pseudotangentiale Momente), 
wahrend 


2 w jw? + p sin w— (2 w? + p) sin Jw? + p+ 2w lw? + p tg > (cos /w® + p — cos w) ==) 


auBer Betracht fallt. 


Fiihren wir die entsprechende Rechnung mit zwei axial bleibenden Momentvektoren durch, 

so findet die von R. Grammel! abgeleitete Knickformel 
w? + p= 7? 

ihre Bestatigung. Es zeigt sich ferner, daB diese Bedingung auch gilt, wenn die Momente in A 
und B stets die Richtung der Stabtangente aufweisen. Diese beiden nichtkonservativen Knick- 
falle muBten eigentlich mit dem kinetischen Stabilitatskriterium untersucht werden. Die 
Schwingungsrechnung wurde fiir den nur durch axiale Torsionsmomente belasteten Stab von 
A. Trésch? durchgefiihrt und bestatigte die statisch berechnete Knicklast. Es ist infolgedessen 
zu vermuten, da auch beim zusatzlich auf Druck beanspruchten Stab und bei tangentialen Tor- 
sionsmomenten die statisch ermittelte Knickbelastung mit der kinetischen iibereinstimmt. 

f)Diskussion der Knickformeln. In den meisten Knickfallen wird die kritische 
Dyname durch eine transzendente Beziehung zwischen w und p festgelegt. Nur bei den Lage- 
rungen II und III (semitangentiales Moment), bei denen die Knickkurven Parabeln sind, laBt 
sich zu jedem w sofort das zugehérige p angeben. Um die Ermittlung der Knicklast auch in den 
ubrigen Fallen zu erleichtern, wird zu jeder Knickgleichung die zugehérige Kurve im (w, p)-Dia- 
gramm wiedergegeben (Abb. 7 und 8) und die Knickfunktion 2 (w, p) = 0 durch ein einfaches 
Polynom modglichst gut angenihert (Tabelle 3) °. 

Bei der Belastung durch ein pseudotangentiales Moment tritt die Torsionssteifigkeit y 


P= 2m + 


in den Knickformeln auf. Es bezeichnen: E den Elastizitatsmodul, G den Schubmodul, m die 
Querzahl, J, und J, = 2 J, das axiale und das polare Trigheitsmoment, wahrend k* ein von 
der Querschnittsform abhangiger Zahlenfaktor ist. Die kritische Belastung wird damit von der 


1 R. Grammel, Z. angew. Math. Mech. 3 (1923), S. 206 und C. B. Biezeno und R, Grammel, Technische 
Dynamik, Bd. 1 S. 608, 2. Aufl. 

2 A. Troésch, Ing.-Arch. 20 (1952), S. 258. 

* Fir die Tabellierung der verhaltnismaBig komplizierten Knickbedingungen der Lagerung IV stand 
mir die programmgesteuerte Rechenmaschine des Instituts fiir angewandte Mathematik der ETH zur 
Verfiigung, wofiir ich dem Schweizerischen Schulrat und Herrn Prof. Dr. E. Stiefel zu Dank ver- 
pflichtet bin. 
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Form des Querschnitts und vom Material abhingig. Wir beriicksichtigen hier nur Kreisquer- 
schnitte (k* — 1) und Stabe aus Flufstahl (m = 10/3): 
. _ mE Ja 10 wi 


ee ie OU, IS 


Beim pseudotangentialen Moment gelten dann die numerischen Schranken w < 0,605 (a) und 
w < 1,21 (b). 

In Tabelle 2 sind zunachst die kritischen Lasten fiir w = 0 und p = 0 zusammengestellt. Als 
Schnittpunkte der Knickkurven mit der p-Achse erhalten wir die Eulerschen Knickkriifte, 
wahrend die Knickmomente w, die von H. Ziegler! gefundenen Werte bestatigen. Der Torsions- 
winkel 2,6 w, gibt einen Anhaltspunkt fiir die GréBenordnung der Knickmomente; es zeigt sich, 
daf} dem Knicken durch Torsion allein im ungiinstigsten Fall IIT, g eine Verdrehung der Welle um 
117° vorausgeht. Tabelle 2 enthalt ferner in einigen Fallen die Reihenentwicklung der Knick- 
funktion fiir kleine Werte von w: sie zeigt, da nicht alle Knickkurven unter einem rechten 
_ Winkel gegen die p-Achse laufen; bei symmetrischer Lagerung ist die Anderung der knickenden 


Tabelle 2. 
Knickfall Gleichung | P(w = 0) = Po| py bestimmt durch W(p = 0) = % | m bestimmt durch | pay 
I (4.14) 39,478 gin! Po 0 4,494 12 Wy = We Pe= P= 3a 
II (4.12) 9,8696 = 3,142 — 
III, s (4.8) 2,4674 — 1,571 — — 
III, q (4.9) 2.4674 cos |p, = 0 0,7854 cos2w,=0 | pp—p=2|u 
III, p (4.11) 2.4674 cos | py= 0 — — Po —P = 2 |u| 
IV, s (4.18) 20,192 te Vp, =Vpo| 3.406 ct gw) = — —wW)| Po — Pp = 2,11 wv? 
LV, ¢ (4.20) | 20,192 tgVp, =Vp.| 2.475 : = 
IV, p (4.23) 20,192 tg Vp, =V Po = == aa 
ies 6 |< (4,24) 9,8696 | sin|/p, =0 2,456 |tgw)+—3?=0| pp—p =1,81 w* 
V.q—4 (4.25) 98696 sin |p, = 0 1571 = Pp —P = 5 w? 
V, p—p (4.26) 9,8696 sin |p, = 0 = — Po —P = 0,2 w? 


Druckkraft durch kleine Torsionsmomente von zweiter Ordnung klein, quasi- und pseudotangen- 
tiale Momente jedoch machen sich bei unsymmetrischer Lagerung schon in erster Naherung be- 
merkbar. 

Die in Tabelle 3 wiedergegebene genaherte Darstellung der Knickgleichungen besteht aus 
einem zweigliedrigen Ansatz, dessen Koeffizienten mit Hilfe der angegebenen kritischen 
| Dynamen ermittelt werden (d. h. die exakte Knickkurve schneidet die Naherung in den Be- 
stimmungspunkten). Da man fir praktische Anwendungen keinen Wert auf die genaue Kennt- 
nis der Bestimmungspunkte legen wird, darf die letzte Stelle der Koeffizienten weggerundet 
werden. Mit der zweiten, meist dreigliedrigen Naherungsgleichung erfassen wir auch den auf 
) Zug und Torsion beanspruchten Stab. Die letzte Kolonne miBt die Gite der Approximation, 


, indem sie den auf © bezogenen maximalen Fehler angibt, der bei einer Berechnung von p aus w 


mit der Naherungsgleichung im Gebiet zwischen den Bestimmungspunkten entsteht. 

Die Knickkurven von Abb. 7 bestatigen den Satz, daB eine Vermehrung der Bindungen 
eine Erhéhung der kritischen Belastung bewirkt: entsprechende Kurven der Lagerungen III, 
II, I und III, IV, I schneiden sich nicht; verschiedenartige Bindungen hingegen lassen sich nicht 
» vergleichen, wie die Schnittpunkte der Kurven III, s und V, s—s, IV und II usw. zeigen. 
| Fiir die kritischen Dynamen eines Stabes mit zweiannadihernd gleichen Biege- 
| steifigkeiten kann man aus den Kurven von Abb. 7 eine gute Naherung gewinnen. In einem 
) Koordinatensystem (w,, pz) mit 

a Re 


UE ae, oes 
a y yo Re B 


(6 Sa) 


| w 


| 1 H. Ziegler, Z. angew. Math. Physik 3 (1952), 5. 112. 
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Abb. 7. Die kritischen Lasten fiir a = B. 
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sind von der gesuchten Kurve zwei Punkte bekannt: die kritische Druckkraft Pp. und das mit 
Hilfe von Abschnitt 3 zu berechnende Moment fiir p = 0. Man zeichnet nun die Knickkurve 


fir « = f, die ebenfalls durch Pp. geht, ein und erhalt daraus die Naherungskurve durch eine 
affine Transformation. 


Wie aus Abb. 8 hervorgeht, ist Knicken durch pseudotangentiale Momente erst méoglich 
bei einer Druckkraft p > pmin (gilt auch fir I, p und II, Pp), Wobei in allen Fallen auBer III, p (b) 


die minimale Knickkraft gréBer ist als 0 Die Knickfalle ITI, p (a) und V, p—p bilden 


_merkwirdig scheinende Ausnahmen von der Regel, da® ein gréGeres zulissiges Torsions- 
moment nur auf Kosten einer kleineren kritischen Druckkraft gewonnen werden kann. Das 
_Wiederansteigen dieser Knickkurven mit zunehmendem |w| mu damit erklart werden, daB, 
wie sich im Verlauf der folgenden Diskussion des Knickfalls III zeigen wird, der Charakter des 
| belastenden Moments sich bei der Verdrehung wandelt. 


| Knickfall III, w,= A* w x’, wy = (1 — A*) wy’: In der Beziehung (4.7a) sind zwei 
mégliche Knickkurven enthalten. Der Schnittpunkt der zum ersten Faktor gehérenden Kurve 

mit der w-Achse ist bestimmt durch tg wy = 1/(1 — 2 A*), derjenige der zweiten durch 

tg wy) = 1/(2 A* — 1). Fir A* = 1/2 liefert somit die erste, fir A* > 1/2 die zweite das 

kleinere Moment. Zum gleichen Ergebnis kommt man bei jedem p < 22/4, so daB wir auf die 

Knickformel (4.7) gefithrt werden. Neben 

der Beziehung tg w) = 1/|1 — 2 A*| zeigt Sir) 

auch die Entwicklung 

qr? 


= p=—Ap=2|1—244||ul 7 


fiir kleine w, daB von den hier erfaBten 
konservativen Momenten das semitangen- 
tiale (A*= 1/2) die gréBten Knicklasten 
ergibt, wahrend das Knickmoment spe- 
ziell im quasitangentialen Fall (A* = 1, 
A* = 0) kleiner ist und allgemein fiir 4*—> 
+ co, A*—> — oo gegen Null abnimmt. 

Knickfall IIT, w, = w(x’ + B* y’), 
wy = 0 [Knickgleichung (4.10)]: Die Ent- 

wicklung fiir kleine 


Ap = 2 |w| (B* (-L) yl + B*?) 4 
und die Schnittpunkte mit der w-Achse 
tg 2 w) = — 1/B* lassen erkennen, dafB 
fir B* > 0 die Kurven auBerhalb der- 
jenigen von I, gq (B* = 0) verlaufen und 
| fur B* > + co gegen die Kurve I, sstreben; 
' fiir B* > — co werden die Knickmomente 
beliebig klein. Beim pseudotangentialen 
/Moment ist B* eine Funktion von w. 7 
' Fassen wir auf einer Parallelen zur w-Achse 
‘im Abstand p jeden Punkt als Schnitt- - é 
| punkt einer der hier diskutierten Knick- Abb. 9, Zur Diskussion des Knickfalls III. 

-kurven auf und bestimmen wir das so zu 

»jedem w gehérige B*(w), so erhalten wir fiir p = 2,090 die in Abb. 9 wiedergegebene Kurve. 
Sie besitzt zwei Schnittpunkte mit B* = tg 2,6 w, eine Tatsache, die sich in dem itber- 
raschenden Wiederanstieg der Kurve III, p auBert. 


Fron =2,090) 


cS! 


G3 OY OS 06 


g) Die elastische Linie des ausgeknickten Stabes. Vom stati- 
»schen Stabilitatskriterium ausgehend versuchten wir die elastische Linie des ausgeknickten 
'Stabes zu ermitteln; beim Anpassen der allgemeinen Lisung an die Randbedingungen 
/ergaben sich die Knickgleichungen, und da diese ja das eigentliche Ziel dieser Arbeit sind, 
| wurde die Berechnung der elastischen Linie an jener Stelle abgebrochen. Wir wollen sie 


nun wenigstens fiir einige einfache Falle zu Ende fiihren. 


| 
| 
| 
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Bei den Lagerungen III und II haben wir von (4.5) auszugehen. Die Projektion der ela- 
stischen Linie auf die (x, y)-Ebene la8t sich [vgl. Ableitung von (4.13)] komplex berechnen: 


t 


Zoot [ite dA 


exp i (A, —@)t—1 exp i (A, —w)t—1 
(A, — @) (A, — @) 


Dabei miissen w und p die zum untersuchten Knickfall gehérige Knickbedingung erfillen. 
Fir III, s und Ilist A= A, + i A; frei wahlbar, da durch die Randbedingungen keine Richtung 
ausgezeichnet ist. Mit A = -- 1 hat ein allgemeiner Punkt P(x(t), y(t), = t) der elastischen 
Linie die Koordinaten 


= mao! — cos @, t) + @, (1 — cos ®,1)], | 
n (4.27) 
y =< (@, sin @, 1 — ®, sin 1) | 
[D,, BD, vgl. (4.21)], woraus fiir p = 0 die Ausdriicke 
1 —cos 2wt sin 2 wt 
oe yt (4.28) 


entstehen. Einige dieser Kurven sind in Abb. 10 wiedergegeben. Jede Kurve laBt sich als halber 
(III, s) oder ganzer (II) Bogen einer Hypozykloide auffassen [Grundkreisradius (®, + ®,)/p, 
Rollkreisradius @,/p]; fir p =0 liegen Bogen von gewohnlichen Zykloiden vor. Bei der 


is Ill, sa: w= —, p= 0; LUNs) 0 — 7 — 105 
1 Sa Ed 
EET sib Oy ering ID 36° Il, bi 0 = 5 
7 7 
Pie 36.8 
IL,ga 
II = ae i z 
fe (ise = ; Ay Se 
aS Ct) EP rie 7c 3 
8 x? 
1 ogy 
TIT qyes— =, p = 0 (markierte Punkte im Ab- 
i @ stand At= 0,2). 
-O4 ~G2 0 02 OF 06 08 40 a ib! 


Abb. 10. Der ausgeknickte Stab bei den Lagerungen II und III fiir die rechts angegebenen Belastungen. 


Lagerung III darf die elastische Linie auch fiir ein quasitangentiales Moment, das die x-Rich- 
tung vorschreibt, aus (4.27) oder (4.28) berechnet werden, mu8 dann aber, weil durch die Rand- | 
bedingung (4.6) das Verhaltnis A;/A, = tg a festgelegt ist, noch um den Winkel « gedreht 
werden [fiir die Kurve ITT, q (a) ist z. B. A;/A, = ov]. 

In den Knickfallen I und V, s—s erhalten wir fiir 


eine (4.13) entsprechende verallgemeinerte Bezichung, in der B und C mit Hilfe der Randbe- 


dingungen durch A auszudriicken sind. Fir A = 100/p wird die elastische Linie dargestellt 
durch 


100 sin ®,/2 100 © sin®,/2 _ @, 
i= (OO, C8 Ob = (, = et aD 1 1 2 
Pe ( 2 1 1 Sin ©, /2 cos M, Pe @, cos 5 Dies B,/2 cos —*) » 


100 
Pp 


of 


COE : sin @,/2 . 
(4) TeiP Sse be @, sin OP, t — @, — ee sin QM, i) (4.29) 


ie _ 100 1 cos 2 wt — cos w 100 sin2wt 
p=: = 3a 2"(q —*) |. i (2+ — “ah | 


sin w sin w 
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wobei Anfangs- und Endpunkt zu t = — 1/2,t =1/2 gehéren (zugehérige Knickkurven in 
Abb. 11). Diein Abb. 12 wiedergegebenen Kurven des Falles V, s—s wurden aus 


0 
Gt p (% cos ®t + Cy eT hos Dt 4 oF (% cos 2 + ®, IER, ae cos » 


in ®,/2 2 2 15in @,/2 2 
100 OD, sin ® 1/2 
y = —(4wsin—-t — @, sin @,t ®, - 5 sin @, t 
> 3 2 nage tere ap) at), (4.30) 
100 1 cos 2 wt — cos w 100 /sin 2 wt 
p02 4 = —— |2w(— — 1) 4 y= Beets | 
8 w? 4 sin w 8 w? \ sin w 
berechnet. 
1 
2 
2 
Nia 
7 
Zz 9 
3 
=f 
—2 
Abb. 11. Der ausgeknickte Stab bei der Lagerung I 3 
| fiir die Belastungen I,a: w = 4,4944, p = 0; Abb. 12. Der ausgeknickte Stab im Knickfall V, s—s fiir die Belastungen 
I, b: w = 2,885, p = 20,0; I,e: w= 1,413, V, a: w = 2,456, p = 0; V,b: w = 1,683, p = 5,0; V,c: w = 1,027, 
_ p = 34,0 (markierte Punkte im Abstand 4t = 0,1). p = 8,0 (markierte Punkte im Abstand At = 0,1). 


(Hingegangen am 20. August 1954.) 


Anschrift des Verfassers: Dipl.-Phys. Max Beck, Zollikerberg bei Ziirich (Schweiz), 
Forchstrafe 110. 
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Uber die Koppelung der Dehnungs- und Torsionsschwingungen 
von umlaufenden Scheiben. 


Von K. Zoller. 


1. Einleitung. Bei einer umlaufenden Scheibe sind die von den Dampfkraften nicht unmittel- 
bar erregbaren Dehnungsschwingungen iiber die Corioliskraft mit den Torsionsschwingungen ver- 
koppelt, und da zudem die Eigenfrequenzen beider Schwingungsarten von der gleichen GréBen- 
ordnung sind, muf die dynamische Festigkeitsuntersuchung der Scheibe beide Reihen von Eigen- 
frequenzen ermitteln!. Abschatzungen haben gezeigt, daf bei den heute zulassigen Drehgeschwin- 
digkeiten von Dampfturbinenscheiben diese Untersuchung einfacher an der ruhenden Scheibe 
durchgefiihrt werden darf, ein Problem, das man heute vollstandig beherrscht. Anders verhalt 
es sich bei hheren Drehgeschwindigkeiten, bei denen die Eigenfrequenzen durch die auftretenden 
Flieh- und Corioliskrafte stark abgeandert werden. Diese Frage hat zuerst K. Yamada? zu beant- 
worten versucht, doch kénnen seine Ergebnisse nicht befriedigen, da seine Ausgangsgleichungen 
einen wesentlichen Fehler enthalten, wie R. Grammel bemerkt hat*. In der vorliegenden Arbeit 
soll daher dieses Problem erneut behandelt werden. Dabei beschranken wir uns auf die Bestim- 
mung der Grundfrequenzen und die Abschatzung der zweiten Eigenfrequenzen im Musterfalle 
einer Scheibe mit hyperbolischem Profil, freiem AuBen- und starr eingespanntem Innenrand. 
Ausgangspunkt der Untersuchung ist das Hamiltonsche Prinzip, aus dem wir zuerst als Neben- 
ergebnis die exakten Differentialgleichungen fiir kleine stehende Schwingungen gewinnen. Die 
Eigenfrequenzen selber werden mit dem Ritzschen Verfahren berechnet. 


2. Das Hamiltonsche Prinzip. Wir beziehen (Abb. 1) die Lage der Scheibenteilchen auf ein 
mit der mittleren Scheibendrehgeschwindigkeit ~ = konst. umlaufendes Polarkoordinatensystem 
in der Scheibenmittelebene. Insbesondere seien 
r, p die Polarkoordinaten der Scheibenteilchen 
in ihrer stationaren Gleichgewichtslage P, auf der 
mit ~ =konst. umlaufenden, also vorgespannten 
Scheibe vom Profil y(r). Bei kleinen, dreh- 
symmetrischen Schwingungen um diese Gleich- 
gewichtslage fiihren die Teilchen eine zusatzliche 
Radialverschiebung (r,t) und eine Tangential- 
verschiebung um den Winkel (r,t) aus. Beide 
werden als klein von erster Ordnung angesehen. 
Wir bezeichnen noch partielle Ableitungen nach 
der Zeit ¢ durch Punkte und Ableitungen nach 
dem Radius r durch Striche, z. B. ist 


: do a 
= a” Raton () 


Aus den Polarkoordinaten r+ 0, p+y der aus- 

Abb. 1. Koordinaten und Geschwindigkeitskomponenten, gelenkten Teilchen berechnen sich ihre Geschwin- 
digkeitskomponenten in radialer und tangentialer 

Richtung zu @ und (r+o) y, so daB das Geschwindigkeitsquadrat bis auf GréBen, die klein 


von dritter Ordnung sind, 
v? = 0? + r? wy? + ((3. 0.)) (2) 
ist. An duBeren Kraften wirken auf die Massenelemente dm = — y r dr dp in radialer Richtung 
§ 


die Flichkraft (r + @) w? dm und eine Komponente 2 w (r + 0) dm der Corioliskraft, deren 
tangentiale Komponente — 2 w @ dm ist. 


' R. Grammel, Ing.-Arch. 6 (1935), S. 256 und C. B. Biezeno u. R. Grammel, Technische D ik, Bd. 2 
S.29ff. 2. Aufl. Berlin-Géttingen-Heidelberg 1953. ‘ot 

2 K. Yamada, Proc. Sec. Japan Nat. Congr. Appl. Mech. 1952, S. 343. 

3 R. Grammel, Ing.-Arch. 22 (1954), S. 145. 
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Um das Hamiltonsche Prinzip anzuwenden, berechnen wir zuerst die kinetische Energie der 
Scheibe (bezogen auf das umlaufende Koordinatensystem) 


Ta [Feta any [yr +9) dr + (3.0). (3) 
b 


Hinsichtlich der Formanderungsenergie der Scheibe und ihrer potentiellen Energie im Flieh- 
kraftfeld beachten wir, daB es sich bei der Schwingung um kleine Verschiebungen der Teilchen 
gegen ihre Gleichgewichtslage P) handelt. Daher sind einander entgegengesetzt gleich und kén- 
nen somit unberiicksichtigt bleiben: die Formanderungsarbeit gegen die Scheiben v orspannungen 
und die Arbeit gegen die stationaéren Fliehkrafte rj@?dm. Von der potentiellen Energie im 
Fliehkraftfeld bleibt so nur der der zusitzlichen Fliehkraft entsprechende Anteil 


D=— | Tourdmo=—2x 7 wt [ yrordr, (4) 
b 
der klein von zweiter Ordnung ist. Ebenso tritt eine Formanderungsenergie L derselben GréBen- 


ordnung infolge der bei der Verschiebung 0, y geweckten Zusatzspannungen 0,,0,, 7 auf. Diese 
gehoren zu den Verzerrungen 


=O, =e, y= ey (5) 


und berechnen sich aus dem Hookeschen Gesetz gemah 


) 
| 
jae +4); (6) 


east ue r 
t= Oya. 
wo E der Elastizitatsmodul, 4 = = die Querzahl und 
E 
° = 30+H) 0 
der Schubmodul ist. Damit wird die spezifische Formanderungsenergie 
1 E / Ceo Vo amen 
gy (Or & + Oy bp + 27 Ee) = oq mH (o'? 4 2 £ { mi) 4 epg (8) 
und somit ; 
L= yrdrdg: (0,8 + Oy & +27 &g) 
1 a FE 7 5 
memes [yr lilt teu % +6) + erty dr. (9) 
F 


SchlieBlich tritt noch die Arbeit der Corioliskraft bei den virtuellen Verschiebungen 6g, (r + @) Ow 


auf: ; 
6A = [] 2m (r +0) p 0g —2W0 (r + @) dy] dm 


a 


=2n-22 ae (i 60 — 0 dy) dr + (3. 0.)). 
b 


(10) 


Das Hamiltonsche Prinzip verlangt nun, dab fiir beliebige Zeitpunkte t,, t, die Variation 
t, ts 
6f(T—®—L)d+fdAd=0 (11) 
fA t 


wird, wobei die virtuellen Anderungen 60, Oy fur t, und t, verschwinden und, wie die Koordi- 
naten ¢, y selber, die geometrischen Randbedingungen starrer Einspannung am Innenrand r = b 
erfiillen miissen: | 
(6¢):, = (6e):, = (Op), = (Oy), = 9, (12) 
60 = oy = 0 fir rb. 
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3. Die Bewegungsgleichungen. Fiihrt man in (11) mit (3), (4), (9) und (10) die Variation aus, 
so laBt es sich durch Teilintegration nach der Zeitt und nach dem Radius r stets erreichen, 
da8 im Integrand von (11) nur die virtuellen Anderungen 60 und dy, aber keine zeitlichen oder 
értlichen Ableitungen davon auftreten. Wahlt man dann t, = t, t, = t + dt, so folgen in bekann- 
ter SchluBweise die Bewegungsgleichungen 


E , fi / i * 
yr = ee (yro +H ye) —KYe — Fol t2oyrhtutyre, | 
(13) 
yP pal elyryy +yry]—2oyre | 
und dazu am freien AuBenrand r = a die Randbedingungen 
ree 
eye San 
ere | fur 7a, (14) 
y=0 | 


Man iiberzeugt sich leicht, daB die Bewegungsgleichungen (13) mit den Grammelschen Glei- 
chungen iibereinstimmen, wobei nur die erste der Gleichungen (13) ein von der Fliehkraft her- 
riihrendes Glied w? y r 9 enthalt, das man bei den heute vorkommenden Drehgeschwindigkeiten w 
mit gutem Grund weglaBt. Die Yamadaschen Gleichungen unterscheiden sich von den Glei- 
chungen (13) dadurch, daB sie in der zweiten Gleichung ein mechanisch sinnloses ,,Fliehkraft- 
glied“ w? y r? y enthalten. 

Gehen wir nun zu drehsymmetrischen, stationdaren Schwingungen iiber, so geniigt es nach (13), 
fur sie den Ansatz 

0 Ka) :cosat,. “p= V(r) sia (15) 

zu machen. Durch Einsetzen von (15) in (13) findet man fiir die Amplitudenfunktionen R, Y 
das Gleichungssystem 


E / / vA 
area (yr R + uy RY —pyR —<R Te ae) a Ie aa 


2awyrR+ 2G [ye hy +yrP] tot yh =0 | a 
und dazu noch die Randbedingungen innen und auSen 
R(b) = V(b) = 0, 
R'(a) +“ R(a) = ¥"(a) =0. ee 


Damit sind wir auf ein Eigenwertproblem fiir die Eigenfrequenzen ~ gekommen, das wir aber 
nicht weiter verfolgen wollen, da sich die Eigenfrequenzen miheloser auf andere Weise ermitteln 
lassen. 


4. Das Hamiltonsche Prinzip bei stationiren Schwingungen. Wir gehen mit dem Ansatz (15) 
fiir stationdre Schwingungen in die Energie- und Arbeitsausdriicke (3), (4), (9), (10) ein und neh- 
men als virtuelle Anderungen 

doe =dReosat, dy =dVsinat. (18) 


Damit kommt 
a 


T=207-0%| yr (R? sin? « t + 1? YW? cos? « t) dr, 
b 


a 


P=—2n7 [ 97 RP cost a tdr, 
fo} 
b 
Te E RR | R a 
L=224-5 [or |ra(R” + 2 u 2 rz) coat +Grt Wt sint a t dr, 


b 
6A =2m-2 70 0 | 97 (WAR cos? a t + ROW sin® a 8) dr 
é 


XXIII. Band 1955 Zoller: Uber die Koppelung der Dehnungs- und Torsionsschwingungen usw. 297 


Fir die weitere Rechnung ist es nun zweckmaBig, dimensionslose GréBen zu verwenden. Ohne 
Anderung der Bezeichnung sollen jetzt alle Langen in Vielfachen des auBeren Scheibenhalbmessers a 
ausgedriickt sein, d. h. wir schreiben im folgenden durchweg einfach wieder r, R, y, b anstelle 
von r/a, R/a, y/a, b/a. AuBerdem fiihren wir die bezogene Drehgeschwindigkeit 

ya 
und anstelle der Eigenfrequenz « den ebenfalls dimensionslosen EKigenwert 


A=a//%% (21) 


ein. Hiermit wird aus (19) 


1 
7 T eer 
iS 2nGa? 2 12 | yr (RE sint a t + 12 Y? cos? 1) dr , 


b 


2 


p= a= 2 yr R? cos? atdr, 

Nae (22) 
tee [Polesats teres cos? a t + r?2 W’2 sin? g t| dr, 

j 1 
(A= ea Vie, [9 PCP OR costae + ROW sin? at) dr. 


b 


Bei der Ausfiihrung der Variation (11) 


ee et = ty — 
6f/(T—®—L)d+ f bAde=0 (23) 
Fe ty 
wollen wir einmal die Integrationsgrenzen t, = 7/24, ty = > + 2 7 )/x wahlen, miissen dann aber 


nach (18) mit Riicksicht auf (12) bY = 0 setzen. Das Ergebnis der zeitlichen Integration hebt 
sich wegen 


t, t 
ip sin? « t dt = if cos? « t dt (24) 


t, zs 
als unwesentlicher gemeinsamer Faktor weg und es bleibt die Bedingung 
1 1 
(2 + 0?) [yr RORdr +240 f[ yr? dRadr 
b b 
1 : (25) 
2 , , fa 1 , ' So 
een rae 6R’ + 4 (ROR + R’ 6R) + | ROR| dr =0 
b 


iibrig. Wahlen wir weiter t, = 0 und t, = 2z/a, so ist jetzt entsprechend OR = 0 zu setzen, und 
die Variationsbedingung geht in die Beziehung 


1 1 t 
Ae fy 8 WOW dr + 2A@ f yr? ROW dr — fy PW OY dr —0 tae) 
5 U : 
uber. 


5. Das Ritzsche Verfahren. Aus den Gleichungen (25) und (26) lassen sich die Kigenwerte A mit 
Hilfe des Ritzschen Verfahrens ausreichend genau bestimmen. Wir gehen hierzu von den k-glied- 
rigen Ritzansatzen 


R= a; x; (r) > VU: DE b; Y);(r) P) (27) 
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also ORa=—- > l0aht OV =, S00; r5 ! 
OR = SOG. OF ss, 0b ty 
aus. Die Genauigkeit der Rechnung wird erheblich gesteigert, wenn wir die Koordinatenfunk- 


tionen {;, ; nicht nur die geometrischen, sondern samtliche fiir R und Y bestehenden Rand- 
bedingungen (17) erfiillen lassen: 


(28) 


ti(b) = 9(b) = 0, \ (29) 
Hl) +4) =H) =0. J 
Diese Eigenschaften haben beispielsweise die Funktionen 
t= (r—b)"+ p (r—by"™ 
y n+n(l—) 
ieee Pi (20) [ict [ee 
und ) =(r—b)"+q (bye 
. n 
rape 1 GF G8)" 
Insbesondere wahlen wir 
; 7 1+ 4(1—d) 
mar—bta(—ih at = — Gaye a mT | 
> (r b) =F Pe (r b) mit Po (1 wee b) [4 ap bu al va b)] ? t (30) 
}, =r—b+ q, (r—})? mit n =—sa-—° 
}, = (r — 6)? + q, (r — b)* mit & =—7a—p- 


Setzt man nun (27) und (28) in die Variationsbedingungen (25) und (26) ein, so kommt 
(28 +) f yr(¥ oj 4) (Ya) dr + 24.2 fy ®(¥ by) (Y da, 2) dr 
rhea ye lBaj¥) (8 docx) + 2 (S a4) (Y dai x) +“ (Yaz) (Y da4) 
we or x) (¥) 6a; z,)] dr = 0 
rome 98 fy (3 by) (3486) de + 24.0 fy r® (a; wi) (3 db.) dr 


— fx (3 jj) (2 db.yi) dr = 0. 


Aus diesen Gleichungen erhalt man in bekannter Weise ein lineares, homogenes Gleichungssystem 
fir die Koeffizienten a; und b;, das mit den Abkiirzungen 


y 1 
A Ned, ots, By = fy uyjdr, 
1 
Cie May pode? 1 
ee LYyt T&ty + uty) +a uy dr , (31) 
i : 1 
ey By = fy 8 yi) dr 
i 


von der Form 


k 
i 


k 
a; [(A? + 2?) Ai; — C;;] + » 6-212 B;=0, 
jou 
(i =1,...h) (32) 
k k 
j= 
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ist. Die Wurzeln der gleich Null gesetzten Determinante dieses Gleichungssystems sind obere 
Naherungswerte der 2k ersten Eigenwerte 2 und damit der entsprechenden Eigenfrequenzen a 
der Scheibenschwingungen. Wie man sieht, zerfallt bei ruhender Scheibe (Q = 0) das Glei- 
chungssystem (32) in ein System, das nur die a; als Unbekannte, und in ein zweites System, das 
nur die 6; als Unbekannte enthalt. DemgemaB zerfallt die Frequenzgleichung in deren zwei, die 
in getrennten Reihen die Eigenwerte Aj, 43, . . . bzw. Frequenzen der Dehnungsschwingungen und 
die Eigenwerte Ay, Az, ... baw. Frequenzen der Torsionsschwingungen liefern. Da fiir beliebige 
Werte Q diese Frequenzen stetig von Q abhangen, kann man danach auch fiir 2 + 0 die Eigen- 
werte in zwei entsprechende Teilfolgen (A;) und (;’) aufteilen, die die Dehnungs- bzw. die torsio- 
nale Folge heifen mégen, wenn auch diese Bezeichnung hinsichtlich der Schwingungsformen der 
Koppelschwingungen eigentlich nur fiir (2Q—» 0 sinnvoll ist. 


1) Eingliedrige Ritzansiatze (k = 1): Schon recht gute Naherungswerte fiir die beiden 
niedersten Eigenwerte jj und Aj erhalten wir aus den eingliedrigen Ritzansaétzen 


R=a,},(r), WY = b, h,(r) - (33) 
Die charakteristische Gleichung lautet in diesem Falle 
Ay, Dy A* + [2? (Ay Dy — 4 Bi) — Ay, Ey, — Cy Dy)? — E,, (Q? Ayy— Cy) = 0. (34) 


Insbesondere bei ruhender Scheibe (Q = 0) ergeben sie die beiden Naherungswerte von oben 


y (Ee 
Axo ~ Fe Aro YS ie . (35) 


11 11 


2) Zweigliedrige Ritzansatze (k=2): Sehr gute Naherungswerte fiir 4), A; und gute 
Abschatzungen von oben fiir die nachsten Eigenwerte Aj, Ay werden durch zweigliedrige Ritz- 
ansatze geliefert. Die Auswertung der Determinante fiihrt auf eine charakteristische Gleichung 
von der Form 


co A® + c,A® + cgA* + ¢,42 +e, = 0, (36) 
deren Koeffizienten sich aus den HilfsgréBen 

S, = Ay, A» — Alp, 

S, = 2 S;— Ay Cyy — Agy Cy + 2 Ajy Cy» 

S3 = (2? Ay, am Cy) (2? Aggy a Cy») as (22 Aj, — 12)” = Cy Coy = Cre “P QQ? (S, — 2? S;) ? 
8S, = Dy Dy, — D3, 

S; — Dy Ex =e Dyy Ey ae Diy Ey, 2 

S, = E,; Ex, or E}, ° 


— 


S, = Ay, By — Aj, By ? 
Ss = Boy (22? Apts C,;) = Bay (2? Aj, — Cy) = 22? S, + By Cy. — By Cy, 
Sy = By Do, — Boo Dy» ; 
Sip = By Ey, — Boy Ex» , 
Su = Ay, By — Aj, By , 
Sip = By ($2? Ay, — Cy) — By (£2? Ayy — Cyp) = $2? Sy + By Cy — By Cy, (37) 
Sg = By Dy, — By Dy ; 
Su = By E,, — By E,, , 
| Sts = Aj, Boy — Ay By, 


eae By, (22 Aj, a Cy») 7za By, (22 Ay» eat C.9) == £2? Sts “I B,; Cy — B Ci ? 
ar By Day — By, Diy ? 


ne 
for) 
rw 
1 


WK 
q 


Sig = By Fy, — By, Ey , 

Sig ae Aj, Bo. — As, By , 
| Seo = B,, ($2? Ay, — Cy») — By (22 Avy a Cy) = 22? Sig + By Cyp — By Ci» 
| Sor = By Dy», — By, Dy, , 

Spo = Bry Ey, — By, E,,, 


we 


a (By Byy ae By, By)? 


— 
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gemah 
Cis ) 
¢ =—S8,S; + 8,8, + 4.2? (— S, Sy + Sy Sys + Sh5 Siz — Sy Sou) » | 
c = S, Sz —S, 8; + S3 Sy + 4 2? (— S; Sip — Sg Sy + Sty Sta + Sto Siz rf 
Bis Sys Sys a Si Si ca Sip Soo — Soo Soy + 4 £228 23) ? | 38) 
cs = S, Sg— S38; + 4 2? (— Se Sip + Sho Sta + Sie Sig — Soo Soo) » 
Cp = Sas J 
berechnen. Bei ruhender Scheibe (Q = 0) spaltet sich die Gleichung (36) in die zwei Gleichungen 
4S, +22 S,(0) + S,(0) = 0 (39) 
fiir die Naherungswerte von Aj, A: und 
At S,—A? S; + S, =0 (40) 


fiir die Naherungswerte von Aj, Ay auf. 


6. Beispiele. Um den EinfluB der Drehgeschwindigkeit auf die Eigenfrequenzen auch zahlen- 
maBig zu erkennen, wollen wir die vorstehenden Gleichungen fiir eine Scheibe von hyperboli- 
schem Profil 


konst. 
are oss (41) 
mit dem Innenhalbmesser b = = bei einer Querzahl wy = + (m = 4) auswerten. Die Frequenz- 


gleichung (34) wird hier 


44 — 2.42 (34,515 + 1,45 064 Q2)" 
+ 589,835 — 9,9906 22 =0. (42) 


Gleichung (42) mége fiir ein Drehgeschwindig- 
keitsintervall O0< Q2 <Q*  ausgewertet 
werden, wo (2* dadurch definiert ist, daB fiir 
Q—» Q* der niederste Eigenwert, als welcher 
sich Aj erweist, gegen Null strebt. Aus den 
Gleichungen (16) und (17) findet man leicht 
die Beziehung 

ORE a (GO) a0), (43) 
Die Wurzeln von (42) fiir verschiedene Dreh- 
geschwindigkeiten sind in der Tabelle 1 an- 
gegeben. 


Tabelle 1. Wurzeln der Frequenzgleichung (42). 


Abb. 2. Naherungswerte fiir die ersten vier Higenwerte bei Die Gite dieser Naherung beurteilen wir aus 
den exakten Werten/},/; beiruhenderScheibe. 
Die genaue Frequenzgleichung fiir die Deh- 
nungsschwingungen! liefert in unserem Fall den Wert Aj = 7,283, so daB der Naherungswert 
A, = 17,6837 mit 5,5°% Fehler noch als maBig gut gelten kann. Wesentlich besser sind die 
Naherungswerte fiir die torsionalen Frequenzen. Die genaue Frequenzgleichung? ergibt Aj = 2; 
dagegen ist der Naherungswert A = 3,1608 nur um 0,61% zu groB. 


hyperbolischem Scheibenprofil (41) mit b = = S/S + . 


1 A.a.O. R.Grammel (1935), S. 262, Gl. (31). 
2 A.a.O. C. B. Biezeno u. R.Grammel, S. 37, Gl. (17). 
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Zweigliedrige Ritzansatze fiihren zu den Naherungswerten von Tabelle 2 fiir die Eigenwerte 


(Abb. 2). 


Tabelle 2. Ndaherungswerte fiir die Eigenwerte bei zweigliedrigen Ritzansdtzen. 
a a ate A 


2 a a5 ae ie 

0 7,32267 (0,545%)| 16,563 (1%) | 3,1460 (0,15%) | 9,47428 (0,525%) 

0,04 71,3230 16,5634 3,1458 ATAL 

1 71,9791 16,7063 3,0108 9,3766 

3 9.2536 17,7356 2,2753 8,6852 

5 11,6873 19,4298 1,4415 71,6948 

fl 14,5028 Zio Loo 0,46658 6,6064 
ORIOL DON 14,9804 21,8782 0 6,4298 


Die in der Zeile Q = 0 der Tabelle 2 angegebenen prozentualen Fehler zeigen, da auch noch 
der zweite Eigenwert /5 in guter Naherung vorliegt, Man bemerkt, daB mit zunehmender Dreh- 
geschwindigkeit die torsionalen Eigenwerte abnehmen, dagegen die Eigenwerte der Dehnungs- 
gruppe ebenfalls anwachsen. Dabei fallt fiir einen Wert Q’ mit 1 < Q’ < 3 der Wert A mit dem 
Wert Ay zusammen (Abb. 2). : 

Der Wert 2 = 0,04 entspricht etwa heute zulassigen Drehgeschwindigkeiten von Turbinen- 
scheiben. Die zugehérigen Eigenwerte unterscheiden sich von denen bei ruhender Scheibe nur 
ganz unwesentlich, so dafi auch auf diesem Wege das iibliche Verfahren der Berechnung der 
Eigenfrequenzen gerechtfertigt ist. Nun ist allerdings zu erwarten, daB die Abweichungen zwi- 
schen beiden Lésungen bei dem Scheibenprofil (41) gerade besonders klein ausfallen!: Daher 
sollen die Verhaltnisse, allerdings nur mit eingliedrigen Ritzansatzen, auch noch fiir die Scheibe 


gleicher Dicke 


yy == konst: (44) 
untersucht werden. Man findet hier die charakteristische Gleichung 
At — 2A? (15,33 343 + 1,46 443 2?) + 110,77 307 — 4,18 260 2? = 0 (45) 


mit den Wurzeln nach Tabelle 3. 


Tabelle 3. Wurzeln der Gl. (45). 


vw yt 


Q A ay 
0 53,1463 2,0451 
0,04 5,1469 2,0449 
1 93,4817 1,8834 
3 17,4642 1,1457 
5 10,1897 0,2445 
5,14629 10,4037 0 


Die Eigenwerte zeigen gleiches Verhalten wie beim Profil (41), insbesondere wird nochmals be- 
-statigt, daB man bei den heute vorkommenden Drehgeschwindigkeiten von Turbinenscheiben 
die Eigenfrequenzen der ruhenden Scheibe nehmen darf. 


(Eingegangen am 20. September 1954.) 
Anschrift des Verfassers: Dozent Dr. K. Zoller, Stuttgart-Sillenbuch, Kirchheimerstr. 21. 


1 A.a.O. R.Grammel (1935), S. 260. 


262 Weidenhammer: Rheolineare Drehschwingungen in Kolbenmotoren. Ingenieur-Archiv 


Rheolineare Drehschwingungen in Kolbenmotoren. 


Von F. Weidenhammer. 


1. Zielsetzung'. Bei der Berechnung der Drehschwingungen? in den Kurbelwellen der Kolben- 
motore geht man von einem System aus, bei welchem die Kurbelwelle ersatzweise als eine nicht- 
gekrépfte, mit einzelnen Scheiben besetzte Welle angesehen wird. Die hin- und hergehenden 
Triebwerkteile bilden sich in diesem Ersatzbild als periodisch veranderliche Drehmassen ab, und 
die zugehérigen Schwingungsgleichungen sind demnach inhomogene Differentialgleichungs- 
systeme mit periodischen Koeffizienten und periodischen Stérungsfunktionen. Die Stérglieder 
dieser Gleichungen riihren im wesentlichen von den Momenten der Gaskrafte her, enthalten aber 
auch einen Bestandteil, der allein durch die Schwankungen der Drehmassen bedingt ist. In 
linearer Naherung fiir die Verdrehwinkel beschreiben diese Gleichungen erzwungene rheolineare 
Drehschwingungen. 

Zur Berechnung der resonanzkritischen Drehzahlen aus diesen Schwingungsgleichungen sucht 
man iblicherweise die periodischen Lésungen der verkiirzten Gleichungen und die zugehdérigen 
Eigenfrequenzen auf. Da namlich hier bei Ubereinstimmung der Drehgeschwindigkeit mit einem 
dieser Frequenzwerte unbegrenzt mit der Zeit anwachsende spezielle Losungen der vollstandigen 
Differentialgleichungen vorhanden sind, liegt genau dann ein Resonanzfall vor, und die so be- 
rechneten Eigenfrequenzen sind unmittelbar die Resonanzfrequenzen. Man beurteilt also das 
Resonanzverhalten ausschlieBlich nach dem Verhalten der speziellen Lésungen der unverkiirzten 
Differentialgleichungen. 

Die vollstandige Lésung eines Differentialgleichungssystems setzt sich jedoch linear aus einem 
Fundamentalsystem von Lésungen der verkiirzten Gleichungen und einer speziellen Lésung der 
vollstandigen Gleichungen zusammen. Die itibliche Drehschwingungsberechnung befaBt sich 
daher offenbar nur mit der Diskussion der speziellen Lésungen der unverkirzten Gleichungen, 
die im weiteren ihrer mechanischen Bedeutung entsprechend als erzwungene Schwingungen 
bezeichnet werden sollen. Im Gegensatz hierzu sollen die Lésungen der verkiirzten Gleichungen 
als ,,freie“ Schwingungen bezeichnet werden, obgleich diese Bezeichnung dem mechanischen Sach- 
verhalt insofern nicht gerecht wird, als auch ohne auBere Gaskrafte schon von den periodisch 
veranderlichen Drehmassen herriihrende, erregende Momente vorhanden sind und es also freie 
Schwingungen im eigentlichen Sinne nicht gibt. Im weiteren sollen ausschlieBlich die Lésungen 
des verkiirzten Differentialgleichungssystemes berechnet werden. Hierzu besteht insofern Veran- 
lassung, als bekannt ist, daB sich derartige ,,freie’ rheolineare Schwingungen unter Umstanden 
ebenso wie die erzwungenen Schwingungen mit der Zeit unbegrenzt aufschaukeln kénnen. Auf 
die Bedeutung dieser ,,freien“* Schwingungen fiir die Drehschwingungen in Kolbenmotoren hat 
bereits G. R. Goldsbrough* hingewiesen. Er hat insbesondere an einem Labormodell fiir eine Ein- 
zylindermaschine die ,,freien“* Schwingungen als gefahrlich anwachsend nachweisen kénnen, 
wobei in seinen MeBergebnissen die ,,freien“* und die erzwungenen Schwingungen eindeutig zu 
unterscheiden sind. Allerdings ergeben sich bei seinen Betrachtungen praktisch keine anderen 
Drehzahlen als kritisch, als die tibliche Beurteilung des Resonanzverhaltens iiber die erzwungenen 
Schwingungen liefert. Jedoch untersucht G. R. Goldsbrough nicht den allgemeinen Fall einer 
Mehrzylindermaschine mit mehreren Drehschwingungsfreiheitsgraden. Heute bereitet es indes 
keine Schwierigkeiten, die Untersuchungen auf den allgemeinen Fall auszudehnen, da in der 
Zwischenzeit Differentialgleichungssysteme mit periodischen Koeffizienten in anderem Zu- 
sammenhang mehrfach untersucht worden sind. So fand E. Mettler> einen mit der Zeit unbe- 
grenzt anwachsenden Typ freier Schwingungen, dessen Frequenzen bemerkenswerterweise nicht 
mit denen der erzwungenen Schwingungen zusammenfallen, sondern sich nach Art der Schwe- 


* Der Verfasser ist Herrn Prof. Dr. E. Mettler fiir die Anregung zu dieser Arbeit zu herzlichem Dank 
verpflichtet. 


3 o et cae Darstellung in C. B. Biezeno u. R. Grammel, Technische Dynamik, Bd. 2, S. 347, 
. util. e 
° Vel. z. B. K. Klotter, Technische Schwingungslehre, 2. Aufl., S. 350, Berlin 1951. 


* G. R. Goldsbrough, Proc. Roy. Soc. London A 109 (1925), S.99 und A 113 (1927), S. 259. 
> E. Mettler, Ing.-Arch. 17 (1949), S. 431. 


XXIII. Band 1955 Weidenhammer: Rheolineare Drehschwingungen in Kolbenmotoren. 263 


bungsfrequenzen aus den Eigenfrequenzen zusammensetzen. Ziel der vorliegenden Arbeit ist es, 
an einem Beispiel zu zeigen, da durch die ,,freien‘‘ rheolinearen Schwingungen unter Um- 
standen zusatzliche, bisher nicht beachtete kritische Drehzahlbereiche entstehen kénnen. 


Den nachfolgenden Betrachtungen haftet eine gewisse Unsicherheit dadurch an, daf man 
in gewissem Umfange dampfende Krifte in die Rechnung einbeziehen muB, da ,,freie Schwin- 
gungen unter Umstanden ganz weggedimpft werden kénnen. Es wird sich jedoch zeigen, dah 
selbst bei vorsichtiger Einschatzung der fiir die Dampfung bekanntgewordenen Zahlenwerte 
noch gefahrlich anwachsende Schwingungen zu erwarten sind. Allerdings werden sich dann nicht 
die unbegrenzten Schwingungsamplituden der linearen Theorie einstellen, sondern es wird vermut- 


lich nur eine gewisse Laufunruhe entstehen, die im einzelnen nur mit einer nichtlinearen Theorie 
beschrieben werden kénnte. 


2. Diskussion der Schwingungsgleichungen. a) Einleitung. Bezeichnet man mit y, den 
gesamten Drehwinkel jeder Kurbel, so setzt sich dieser aus dem mittleren Drehwinkel wt 
(w Drehgeschwindigkeit der Welle, t Zeit), dem Kurbelversetzungswinkel y, und dem Verdreh- 
winkel #,(t) zusammen: 


Y, =wt+y, + %,. (1) 


Die Bewegungsgleichungen fiir Drehschwingungen um die gleichférmige Drehung der Kurbelwelle 
lauten dann! 

oa ly ae, Ch RS Rd 
Yilpr) Ox 4 may (w + 01)? + BO, + cn—1, 2 (Ox — Ox-1) — Ch, n41 (Ox41 — Ox) = Mi, 


02 een (2) 


Hierin bezeichnen wie immer im folgenden: Y,(y;,) die Drehmassen, c;,;,,1 die Verdrehsteifig- 
keiten zweiter Art mit c,,,,,; = 0,f die Dampfungsmomentkonstante, M;, die AuBeren Momente 
der Gaskrafte einschlieBlich der stationdéren Mittelwerte und Punkte Ableitungen nach der Zeit. 


Die Drehmassen Y’, sind periodische Funktionen der Drehwinkel y,, und nur fiir das Schwung- 
rad ist einfach das Massentragheitsmoment ¥Y;, = J = konst. einzusetzen. In die Gleichgewichts- 
aussage (2) sind dimpfende Momente durch den geschwindigkeitsproportionalen Ansatz [ 0}, auf- 
genommen, obgleich dieser rechnerisch einfache Ansatz den tatsachlichen Verhaltnissen nicht voll 
gerecht werden kann. Doch findet man nur fiir diesen Ansatz einige Zahlenangaben in der 
Literatur, so da man auf diese Naherung angewiesen ist. — 

Beschrankt man sich auf kleine Verdrehwinkel 9, und kleine Verdrehgeschwindigkeiten 0j, 
so kann man (2) linearisieren und erhalt mit 


Yr= wt + x, (3) 
die Schwingungsgleichungen in der von F. Kluge® angegebenen ausfihrlichen Form: 
eee kA : 2 ee, : 
Vpn) 0, + —*w & 4 Ze Eh, + BOL — Ce, % Oe—-1 + (Ch—1, e+ Ch, h+1) Ox 
dy, 2 dp 


— Ch, +1 Pe 41= Mi(pe)— > (k-=.0; eeeeeces, 1h): (4) 
Zur Berechnung erzwungener Schwingungen kann man diese Gleichungen durch Fortlassen der 
linksseitigen zweiten und dritten Summanden jeder Gleichung (4) weiter vereinfachen, wie 
R. Grammel und W. T. Koiter® gefunden haben. Da hier jedoch die Lésungen der verkiirzten Glei- 
chungen aufgesucht werden sollen, liegt ein anderes Problem vor und man kann die Schlubweise 
von R. Grammel und W. T. Koiter zunichst nicht tibernehmen. Es muf daher vorerst ohne Ver- 
nachlassigungen weitergerechnet werden, wobei sich dann aber die Grammelsche Naherung als 
auch hier berechtigt erweisen wird. ~ 

Aus den Gleichungen (4) liest man die eingangs erwahnte Tatsache ab, da auch in dem nur 
theoretischen Falle einer rotierenden Kurbelwelle ohne Antrieb durch die Gaskrafte (M, = 0) 
hereits erzwungene £chwingungen auftreten, was durch die periodischen Schwankungen der Dreh- 
massen bedingt ist. 


1 C. B. Biezeno u. R. Grammel, Technische Dynamik, 2. Aufl., Bd. 2, S. 356, Gl. (11). 


2 F. Kluge, Ing.-Arch. 2 (1931), S. 137. 
3 R. en “7, angew. Math. Mech. 15 (1935), S. 47 und W. T. Koiter, Proc. Kon. Nederl. Akad. Weten- 


schappen, Amsterdam, Ser. B. 54 (1951), S. 464. 
18 
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wo dY, 


Zur Berechnung der Lésungen der durch Fortlassen der Zwangsmomente M,.(px) =F cig aa 


(k =0,1,...,n) verkiirzten Gleichungen (4) ist es zweckmaBig, die Gleichungen durch Ein- 
fiihrung neuer Koordinaten zu vereinfachen. 

Im weiteren werde vorausgesetzt, daB8 die Drehmasse des Schwungrades Ut = J sehr gros 
im Vergleich mit den anderen Drehmassen ist, so das das Schwungrad keine Drehschwingungen 
ausfiihrt (9,, = 0) und die zugehdrige Schwingungsgleichung mit dem Index n in (4) fortbleiben 
kann. 

b) Beseitigung der Geschwindigkeitsglieder. An Stelle der Verdrehwinkel %;, sollen 
andere Koordinaten d, so eingefiihrt werden, daB die Summanden mit den Verdrehgeschwindig- 
keiten fortfallen. Zu dem Zweck werden, wie in der Theorie der linearen Differentialgleichungen 
ublich, die Transformationen 


t/t t 

eee Dia eles (5) 
Si di, e 0 aa di, F(X) 

F(t) 

in die Gleichungen (4) eingefiihrt. Dann erhalt man fiir die neuen Koordinaten d; in denjenigen 
Fallen einfachere Differentialgleichungen, in denen die Exponentialfaktoren ganz aus den Glei- 
chungen verschwinden. Dies ist streng aufer fiir f = 0 dann noch der Fall, wenn samtliche 
Drehmassen Y,(p;) (k = 0, 1,..., — 1) gleich sind. In sehr guter Naherung kann man jedoch 
die Transformationen auch dann verwenden, wenn die Y,(y,) zwar auf Grund der Kurbelver- 
setzungen nicht identisch sind, aber doch gleiche Mittelwerte haben, wie dies bei den Motoren 
iiblicher Bauart der Fall ist. Zerlegt man namlich mit R. Grammel die periodischen Drehmassen 


0 
e 


in ihren harmonischen Mittelwert und eine periodische Schwankung um diesen Mittelwert von | 


der GroéBenordnung ¢: 


1 1 
ern 1 AN 6 
7 = all + eeu) (6 
so nehmen die in (5) vorkommenden Exponentialfaktoren die Werte 
t 
Bice Be Bt 
SR a Ce eC 
e 0 =e + 0(Be) 


an. Da man sich bei der Berechnung erzwungener Schwingungen auf eine Theorie beschranken 
kann, in welcher nur lineare Potenzen des Parameters ¢ beibehalten werden, liegt es nahe, auch 
bei der Berechnung der ,,freien“’ Schwingungen in der gleichen Naherung zu bleiben. Wegen der 
Kleinheit der Dampfungskonstanten f kénnen dann in guter Naherung auch Glieder der GréBen- 
ordnung fe vernachlassigt werden. Um diese Voraussetzung kenntlich zu machen, werde im 
weiteren B = Be geschrieben. Dann sind alle Exponentialfaktoren gleich, sofern simtliche Dreh- 
ae die gleichen Mittelwerte besitzen, und es gelten an Stelle von (4) die vereinfachten Glei- 
chungen 


d+ dy, er —( a y| —— Gi pees ae (Ch k+ Chk ) =i Diy eG) (7) 
2 Uy 2, Ue, Y |, ee —], ok+1 yy k, k+ 1 l%, a —- 
(k= 071,25, 2 == 1) 


fiir d,(t). Offenbar kann man in (7) der Kleinheit von 6 wegen die Summanden (¢ 8/2 Y,)? strei- 


chen. Auch die Glieder 
2 


yw \2 
(Gy =fanow 


enthalten nur héhere als erste Potenzen des Parameters ¢ und bleiben daher ebenfalls unberiick- 
sichtigt. Man kann also in einer Naherung, die nur lineare Potenzen des Schwankungsparameters 
beibehalt, wie in der Theorie der erzwungenen Schwingungen von vereinfachten Gleichungen aus- 


gehen. Zur Berechnung der Lésungen von (7) ist eine noch weitergehende Vereinfachung durch 
Ubergang auf Hauptkoordinaten zweckmaBig. 


c) Hinfihrung von Hauptkoordinaten. Zur Abkiirzung der formalen Rechnung bei | 
der Einfithrung von Hauptschwingungskoordinaten wird das vereinfachte Differentialgleichungs- | 


system (7) in Matrizenform geschrieben: 


te 03 (8) 
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Hierin kommen der Lésungsvektor ¢ und die symmetrische Matrix 9 vor: 


[ dy eS ete | 
; a, Pete hee, 
= > )) f — oe Co1 ah Sig = Oy | | . (9), (10) 
‘ 7 \l ee es Meee | } 


Alle Elemente der Matrix Q{ seien in Fourierreihen mit dem Argument wt entwickelt und ins- 
besondere die jeweiligen konstanten Anteile, die nullten Fourierkoeffizienten, abgespalten, so 
daf die Zerlegung 

MW = AW, + A, (11) 


entsteht, wobei 2, nur konstante Elemente enthalt und Q{, nur periodisch von der Zeit abhingige 
Elemente besitzt. Fir die Transformation auf Hauptkoordinaten ist es zweckmabBig, (8) als die 
Eulerschen Gleichungen des kinetischen Potentials 


ee ives 
b=— zt tt sr Uh +U)e (12) 


aufzufassen. In (12) bezeichnet ein Akzent wie iiblich den Ubergang zur gespiegelten Matrix. 
Es liegt nahe, in (12) eine Hauptachsentransformation zur Einfiihrung von Hauptschwingungs- 
koordinaten durchzufiihren, um das Gleichungssystem (8) hinsichtlich der Funktionen d, mit 
konstanten Koeffizienten zu entkoppeln. Fir konstante Drehmassen wire das Gleichungs- 
system dann vollig entkoppelt und im Falle veranderlicher Drehmassen jedenfalls soweit wie 
moglich vereinfacht und daher einer weiteren Auswertung leichter zuganglich. Zur Durchfiihrung 
der Hauptachsentransformation hat man die Eigenschwingungsformen ¢ von 


r+ %r=0 (8,0) 
aufzusuchen und macht zu dem Zweck den Ansatz fiir Eigenschwingungen mit der Eigenkreis- 
frequenz x: 

Coe 1: (13) 


Dann sind die Eigenkreisfrequenzen aus der Frequenzgleichung 
det |, — x? ©| = 0 
zu bestimmen, worin & die Einheitsmatrix bezeichnet. Wenn man die Frequenzen bestimmt hat, 


lassen sich mit dem Ansatz (13) die zugehérigen Schwingungsformen aus (8,0) berechnen. Die so 
gefundenen n Losungsvektoren 7, seien normiert, so daB gilt 


ae (ie Pei VSR IK 50 oy Wo 
2S ier ey ET ne ey et 
Dann kann man mit diesen Eigenvektoren die Transformationsmatrix 
| 
k= 1 lo —— Sas | (Cin) (14) 
| 
mit den konstanten Elementen a;, zusammenstellen und durch die orthogonale Transformation 
r=Xy (15) 


eer 


zu den Hauptkoordinaten 


| (16) 


wi | 


iibergehen. Fiihrt man die Transformation (15) in das kinetische Potential (12) ein, und beachtet 
x’ X =, so erhalt man einfacher 


1 ? ° 1 ff 1 / 
L=— Z)'y+zVByt+zy Gy, (17) 


j= 


1 L. Collatz, Eigenwertaufgaben, S. 246, Leipzig 1949. 
18* 


266 Weidenhammer: Rheolineare Drehschwingungen in Kolbenmotoren. Ingenieur-Archiv 


worin 


WAR UX=( (18) 


0 Met 
die Diagonalmatrix der Eigenkreisfrequenzquadrate und 
& = (gux-(t)) = x SA 4 (19) 


eine symmetrische Matrix mit periodisch von der Zeit abhangigen Elementen Snw(t) bezeichnet. 
Da in (17) die ersten beiden quadratischen Formen rein quadratisch geworden sind, ist das zu- 
gehérige System Eulerscher Differentialgleichungen nur noch in den Gliedern mit zeitabhangigen 
Koeffizienten gekoppelt: 


n—1 


VE + Hk Vb + >) Brn (t) Yu —4() (k (0h Ilse ou n—l). (20) 
k’=0 


Da die Funktionen g;,,, allgemein periodische Funktionen sind, ist (20) ein entsprechend allgemeines 
Hillsches Differentialgleichungssystem fiir die Hauptschwingungskoordinaten. Uber die Lésungen 
kann man allgemeingiiltig zunachst nur sagen, da sie nach dem Floquetschen Theorem? periodisch 
zweiter Art sind und daB sich die allgemeine Lésung, wenn nicht Ausartungen eintreten, aus 
Partikularlésungen der Form 
n—1 == yp = 
Vee (C 8! 24d, (MiG je | 892", 9 (on) ene 0 er 
1=0 
mit den 2 n Integrationskonstanten C, zusammensetzen. In den Partikularl6sungen kommen 
n Paare charakteristischer Exponenten -+ g; und periodische Funktionen @,)(t) vor’, so daB 
mit der Zeit unbegrenzt anwachsende Lésungen méglich sind, wenn auch nur einer der charakte- 
ristischen Exponenten einen nichtverschwindenden Imaginarteil besitzt. Zur Beurteilung des 
Verhaltens der ,,freien“ Schwingungen muf man daher diese charakteristischen Exponenten 
kennen. Ihre Berechnung ist zwar grundsatzlich z.B. mit unendlichen Determinanten méglich, 
doch ist diese Rechnung im allgemeinen weitlaufig. Im weiteren wird daher nur ein typisches 
Beispiel durchgerechnet, bei welchem man auf in anderem Zusammenhang durchgefiihrte Be- 
rechnungen des charakteristischen Exponenten zurickgreifen kann. 


3. Die Vierzylinder-Zweitakt-Maschine. Als Beispiel werde eine Vierzylinder-Maschine mit 
unter sich gleichen Getriebeteilen und mit unendlich groBem Schwungrad betrachtet. Dann 
erfolgen die Drehschwingungen mit vier Freiheitsgraden. Die Torsionssteifigkeiten dieser homo- 
genen Maschine werden mit c,,4,=¢ (k = 0,1, 2) zwischen den Kurbeln allgemein und mit 
C3,4 =¢/2 zwischen letzter Kurbel und Schwungrad speziell in die Rechnung aufgenommen. Die 
Veranderlichkeit der Drehmassen werde vereinfacht fiir sehr lange Schubstangen berechnet. In 
dieser sogenannten Frahmschen Naherung* erhalt man allgemein 


W=R+5 (1—cos2an), 


wenn R und R’ reduzierte Tragheitsmomente bezeichnen, die sich aus den Massen der Kurbel, 

der Pleuelstange und des Kolbens zusammensetzen. Zur Abkiirzung werden noch 

__ RR _ Ri 
RR+R Oo 


@=R4+2, ¢ a (22), (23) 


gesetzt, so da® ¢ den relativen Anteil der hin- und hergehenden Massen an der ganzen Drehmasse 
kennzeichnet. Auf Grund einer bei den Vierzylinder-Zweitakt-Maschinen iiblichen Bauart sind 


i Vgl. 7165 1B, de Horn, Gewohnliche Differentialgleichungen, Sammlung Schubert, Bd. L, § 55, Leipzig 1905. 
2 Diese Aussage ist in einer Veroffentlichung von J. O. Fleckenstein, Comment. Math. Helv. 15 (1943), 
S. 367 enthalten, kann jedoch dort infolge eines Schreibfehlers nicht herausgelesen werden. Auf S. 373, Zeile 10 
von oben muf es heifen: 

"_, (L+ x”) tg? ad. 
sk, ! Ds 

! ier! 0 te? xz 3s — x2 
woraus das paarweise Auftreten der charakteristischen Exponente h bei all i Val - 
funktionen in der dortigen Gl. (1) folgt. : ne eee ee ae ee 


3 C. B. Biezeno u. R. Grammel, Technische Dynamik, 2. Aufl., Bd. 2, S. 294, Gl. (4). 


=0, 
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die Kurbeln um die Winkel 4a =2/2, X. = 32/2, y; = ma gegeneinander versetzt und alle Mittel- 
werte 0, =@ gleich, so dab 
il 1 


1 
es 2 Bake 
Y, wy, 6 | + € cos 2 wt) + O(c?) , v= w= @ | 


1 — cos 2 wt) + O(e?) (24) 


wird, wenn nur lineare Potenzen des Parameters ¢ beriicksichtigt werden. Um die Hauptkoor- 
dinaten einzufiihren, hat man zunichst die Matrix 9{ von (8) aufzuspalten, was in der Naherung (24) 


auf 
( ol) 0 a 
es Naa et ea 
ere ainthGe ie Boe a (25) 
Ct eae } 
und 
1 0 0 0 \ 
Q= = =I 0 
W, = ecos2wt eke ie RO 0 | (26) 
0 6" 0 37 } 


fiihrt. Die Eigenwerte der Matrix Q{, und damit die Eigenfrequenzen lassen sich leicht und allge- 
mein mit den Frequenzfunktionen von R. Grammel angeben. 
Nach kurzer Rechnung sind damit die Elemente der Diagonalmatrix durch 


k | 0 | il | 2 3 
C, | 0,084| 0,81 | 2,15 | 3,46 


Cc ° 
“ie = 6 be mit den Zahlenwerten 


(27) 


gegeben und man kann nun aus (8,0) die Eigenvektoren berechnen und die Transformations- 
matrix (14) zusammenstellen. Damit ist auch die symmetrische Matrix ® von (19) zu berechnen, 
wobei Einzelheiten dieser Rechnung kein Interesse bieten: 
pti 0,536 —1,980 0,000 
—0,658 —0,737 —0,865 


c 
6 = 6 0,587 | 
Die Differentialgleichungen in den Hauptkoordinaten (20) haben hier also die spezielle Form 
d? ,, { *% \? : c ; 
ot Ga + ( 7 Vee cos sD 6 8H vs Ol)e (be 0.1 23) (29) 


wenn durch 
Set 


eine dimensionslose Zeitzahlung eingefiihrt wird und die g,, hier die Zahlwerte der in (28) vor- 
kommenden Matrixelemente bezeichnen. Die Gleichungen (29) entsprechen formal genau den 
ausfiihrlich untersuchten Stabilitatsgleichungen von E. Mettler!, so daB man alle Ergebnisse in 
dem hier benétigten Umfang dem Schrifttum entnehmen kann. Insbesondere kann man die 
Werte der charakteristischen Exponenten zahlenmaBig angeben und kann damit tber das 
Wachstum der Lésungen von (29) entscheiden. Aus diesen Ergebnissen entnimmt man, dal 
Drehgeschwindigkeiten in der Nahe der Werte 


_ Rr 
Oe 2p 


) (30) 


zu anwachsenden Lésungen y;, der Gleichungen (29) fiihren. Bei diesen Werten w, hat man also 
zu prifen, ob die Verdrehwinkel 


1 FE, Mettler, Ing.-Arch. 17 (1949), S. 426. 
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endlich bleiben oder nicht. Zu dem Zweck ist die Transformation (15) 
3 
di, = > Mi Yi 
i=0 


in (5) mit den Funktionen (21) einzusetzen, wobei zu beachten ist, daB in (29) die Zeitzahlung 
s =2wt verwendet wurde: 


Wy.) pea 3 3 —=1 0,20 
o, = Vane é 20° SS ai (s C, A I ie @;(t) + Cyrie ©, sd) A 
Y(t) i=0 i=0 


Falls nun auch nur ein charakteristischer Exponent 9; einen Imaginarteil von einem solchen 
Betrage hat, daB die Ungleichung 


= Si 0, 20t 


Be 

[Im 011 — Foo 0 (31) 
erfiillt ist, treten unbegrenzt anwachsende Schwingungen ein. Derartige ,, Resonanzen“ der freien 
Schwingungen sind gemaB (30) bei den kritischen Drehgeschwindigkeiten 


Wo — = (k — 0, iN 2, 3) (p Zz 0 ganz) (30,1) 


und zusatzlich bei den nach Art der Schwebungsfrequenzen gebildeten 


Hy + %] \(k — 0; 1, 2, 3) 
2 p \(=09, 1, 2, 3) 


méglich. Allerdings haben nur die zu p=1 gehérenden Werte Bedeutung, da fiir gréBere Ord- 
nungszahlen p der |I m g)| schnell kleiner wird und in der verwendeten Naherung 0 (¢?) tiberhaupt 
fiir alle Werte p > 1 verschwindet. 

Mit den Werten (30,1) braucht man sich nicht weiter zu befassen, da die gleichen Werte in 
bezug auf erzwungene Schwingungen kritisch sind und daher bekannt und ausfihrlich untersucht 
sind. Hingegen sind die Drehgeschwindigkeiten (30,2) in der Theorie der erzwungenen Schwin- 
gungen ganz ohne Bedeutung und daher bisher nicht beachtet worden. 


(k +1, p>0 ganz) (30,2) 


Im weiteren werden daher nur kritische Drehgeschwindigkeiten in der Nahe von @)= : (x, + x1) 


daraufhin untersucht, ob diese Werte auch bei Vorhandensein der Dampfung noch zu gefahrlich 
anwachsenden Schwingungen fiihren. Diese Frage kann durch Auswertung der Ungleichung (31) 
dahingehend beantwortet werden, da man diejenigen Werte der kritischen Drehgeschwindig- 
keiten We angibt, welche selbst gerade noch zu anwachsenden Schwingungen fiihren und nur 
kritische Werte einschlieBen. Zu dem Zweck kann man auf in anderem Zusammenhang von 
G. Grammel+ gefundene Ergebnisse zuriickgreifen, welche eine allgemeine formelmafige Aus- 
wertung in dem hier bendtigten Umfang gestatten. Mit den Konstanten g;, von (28), den Mittel- 
werten © aus (22), dem Parameter ¢ von (23) und den Zahlen ¢;, von (27) erhalt man fiir kritische 
Drehgeschwindigkeiten in der Nahe von 


die Grenzwerte 


Oc = Lte|/ = i ae p? 
it =<) AVG VG) eG ea 


die einen ganzen Bereich von resonanzkritischen Drehgeschwindigkeiten von der Breite 


2 = 
Skl ie 
ew ws See Oe 
: | (Ven + Ve)? Vee, OC? «5 (32) 
einschlieBen. Da die Eigenfrequenzen mit wachsenden Ordnungszahlen k und I schnell zunehmen, 


ist der Bereich, der zu den kleinstméglichen Indizeswerten k —0, 1 =1 gehort, gemaB (32) der- 


jenige Bereich, der die gréBte Frequenzspanne umfaBt und ist damit wohl der gefahrlichste. Er 
soll daher auch zahlenmaBig berechnet werden. 


1 G. Grammel, Ing.-Arch. 20 (1952), S. 170. 
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Zu dem Zweck mu8 man die Dampfungszahl B von (2) kennen. H. Holzer} setzt B =B, Ow 
und gibt fiir den Zahlenfaktor B, den Wert von etwa 0,04 an. Messungen von J. L. Mansa? 
haben diesen Wert gut bestatigt und jedenfalls keine gréBeren Zahlenwerte ergeben. Mit diesen 
Angaben lassen sich die Grenzen des Bereiches kritischer Drehzahlen berechnen und in einem 
Diagramm darstellen (Abb. 1), in dem als Abszisse w/w, und als Ordinate der Parameter ¢ von (23) 
als MaB fiir die GréBe der Drehmassenschwankung aufgetragen ist. Je starker der Anteil der hin- 
und hergehenden Massen, d.h. je groBer ¢ ist, um so gréBer ist der Bereich der kritischen Dreh- 
geschwindigkeiten. Aus Abb.1 ist weiter zu entnehmen, daB fiir den angegebenen Daimpfungs- 
beiwert 6, =0,04 kaum ein Einflu8 der Dampfung im Sinne einer Verminderung der Gefahrlich- 
keit der Schwingungen zu erkennen ist. Selbst fiir einen angenommenen fiinffachen Wert 6, = 0,20 
bleibt immer noch ein wesentlicher Bereich von Drehzahlen kritisch. Allerdings mu dann der 
Parameter ¢ mindestens den ,,Schwellwert“ ¢, 0,23 annehmen, um zu gefahrlichen Schwingun- 
gen zu fithren. Der Schwankungseinflu8® der Drehmassen muB also einen gewissen Mindestwert 

é gegentiber dem zeitlich konstanten Mittelwert erreichen. 
Aus Abb.2 ist die gegenseitige Lage der ein- 


O4¥ 
zelnen kritischen Frequenzbereiche héherer Ordnung 
fir k=0, 1=1 und p=1, 2, 3, 4 zu erkennen, wobei 
a diese Bereiche allerdings in der zugrundeliegenden 
Naherung mit nur ersten Potenzen von ¢ nicht zu be- 
ee rechnen sind. Sie diirften keine praktische Bedeutung 
a haben, da sie nur auf sehr viel langsamer anwachsende 
g 
0,04 
@ 
@ oO 
4 Ranere Gh Go: 40 af fiat Soe p=t 
Abb. 1. Bereich kritischer Drehgeschwindigkeiten fiir Abb. 2. Lage der Bereiche héherer Ordnung (p > 1) 
Oy = (v0 +). DampfungseinfluB: ——— : B, = 0; fiir = aie > p=1, 2, 3,4. 
—-—-— : B, = 0,04 nach H. Holzer; —— —: B, = 0,2. 


Schwingungen fiihren und daher von der Dampfung vermutlich ganz ausgeléscht werden. Uber- 
dies haben sie eine verschwindend kleine Frequenzbreite. 


4. Zusammenfassung. Die Berechnung der Drehschwingungen in den Kurbelwellen der 
Kolbenmotore fiihrt in linearer Naherung fiir die Verdrehwinkel und Verdrehgeschwindigkeiten 
auf lineare Differentialgleichungen mit periodischen Koeffizienten, wenn man den Ein- 
flu8 der hin- und hergehenden Triebwerkteile mitberiicksichtigt. Dieses Differentialglei- 
chungssystem ist inhomogen mit periodischen Stérungsfunktionen auf Grund der auferen 
Momente der Gaskrafte und als Folge des periodisch mit dem Verdrehwinkel schwankenden 
Einflusses der Triebwerkmassen. Die Lésung dieses Gleichungssystems setzt sich zusammen 
aus den erzwungenen Schwingungen als den speziellen Lisungen der vollstandigen Glei- 
chungen und den Lésungen der verkiirzten Gleichungen, kurz ,,freie“ Schwingungen genannt. 
Bei der Betrachtung dieser freien Schwingungen an dem Beispiel einer Vierzylinder-Zweitakt- 
Maschine ergeben sich zusatzliche kritische Drehgeschwindigkeiten wy = (xx +%1)/2 (ene 3) 
auch wenn man den Dampfungseinflu8 mit den Zahlenwerten von H. Holzer vereinfacht beriick- 
sichtigt. Ob diese neu berechneten kritischen Drehzahlen tatsachlich zu einer gefahrlichen Lauf- 
unruhe im Betrieb der Maschinen fiihren, kann wohl nur experimentell beantwortet werden, da 
man rechnerisch eine nichtlineare Theorie mit EinschluB der Dampfung heranziehen mufte, was 
derzeit kaum Aussicht auf Erfolg bietet. 


Aus dem Institut fiir Mechanische Schwingungstechnik der Technischen Hochschule Karlsruhe. 
(Eingegangen am 24, September 1954.) 
Anschrift des Verfassers: Priv.-Doz. Dr. F. Weidenhammer, Karlsruhe i. B., Hertzstr. 16, Bau 40. 


1 Vgl. F. Séchting, Berechnung mechanischer Schwingungen, S. 226. Wien 1951. 
2 J. L. Mansa, Die Berechnung der Dampfung von Drehschwingungen einer Flugmotorkurbelwelle, Diss. 


Karlsruhe 1932. 
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Kreisférmiger Stollen unter Temperaturbeanspruchung. 
Von K. Hirschfeld. 


1. Einleitung. Durch Messungen und Beobachtungen hat man festgestellt, daB im allgemeinen 
in einer Tiefe von 10 bis 20 m unter der Erdoberflache das Jahresmittel herrscht, und von da ab nach 
dem Erdinnern zu die Temperatur ansteigt. Im Durchschnitt nimmt sie auf 100 m Tiefe um 3° C 
zu, so da man die geothermische Tiefenstufe, innerhalb der die Vertikalstrecke eine Temperatur- 
zunahme von 1° aufweist, mit etwa 33 m annehmen kann. Dennoch ist es nicht méglich, allgemein- 
giiltige Zahlen anzugeben, da die geothermische Tiefenstufe vom geologischen Aufbau der Erdrinde 
abhingt. In der mittleren Schwabischen Alb hat man 11 m, in Kanada und Siidafrika dagegen 
125 m gemessen. Der Grund ist darin zu suchen, daB in alten Schollen, die sich seit langer Zeit 
orogenetisch und magmatisch in Ruhe befinden, die Temperatur sehr langsam zunimmt, wahrend 
junge vulkanische Gebiete mit verhaltnismaBig hoch emporreichendem Magma einen gréBeren 
Temperaturanstieg aufweisen. Durch den Zerfall der radioaktiven Elemente wird der Warmevorrat 
standig erginzt, so dafs man aus dem Gehalt der Gesteine an Uran, Thorium und Kalium auf die 
entwickelte Warmemenge schlieBen kann. Danach ist 1 cm? Gestein in der Lage, bei 


Granit Balsalt Peridotit Sedimenten 
40 10 5 10—30 


Kalorien im Jahr zu entwickeln. Die Unsicherheit in der Kenntnis der Gesteinszusammensetzung 
laBt aber uber den Verlauf der Temperaturkurve nur Vermutungen zu. 


Fur die Warmeberechnung von Stollen, die die vorher betrachteten Tiefen bei weitem nicht 
erreichen, ist man gezwungen, sich auf MeBergebnisse zu stiitzen, die an ahnlichen Bauwerken fest- 
gestellt worden sind. Eine der bekanntesten Naherungen ist die von Niethammer in Fantolis Schrif- 
ten wiedergegebene Gleichung?, die sich auf Beobachtungen am Sempione-Stollen in den Alpen 
stutzt. 


d = 11,8 — 0,00457 Hy + 0,025 He . 


Darin bedeuten Hy, die Hohe tiber dem Meeresspiegel in m und Hg der Abstand zwischen Stollen- 
achse und Gelandeoberflache in m. Dem Ausdruck liegt auBerdem die Annahme zu Grunde, daB 
die Temperatur des Gesteins in einer gréBeren Tiefe als 15 m konstant bleibt. Nimmt man z. B. an, 
daB Hy = 500m und He = 100m betragt, so wird # = 12°C, oder fiir Hy —1000m und 
He= 200m wird 7 = 12,2- C. 

Nach Hiitte I, 27. Auflage S. 593 ergibt sich die Warmettbergangszahl fiir Wasser in geraden Rohr- 
leitungen naherungsweise zu 


o = 1755 (1 + 0,015 8,,) 


w0,87 keal 
013 | m2 nec 2 


wobei fiir @,,, der Wert 
US, = 0,9 Ow = 0,1 Or 
empfohlen wird. 
In diesen Ausdriicken sind Jy die Wassertemperatur, 3p die Temperatur der Rohrwandung, 
w die Geschwindigkeit des Wassers in m/s und d der Durchmesser des Rohres in m. 

Wahlen wir als mittlere Werte dy = 7° C, Op = 12° C und damit %,, = 7,5° C, und setzen wir 
weiter w = 2 m/s und den lichten Rohrdurchmesser d = 3,0 m, so errechnet sich die Warmeiiber- 
gangszahl zu 

a ~ 3100 |_*e 
7; Far (| : 

Dieser hohe Wert der Warmeiibergangszahl liegt fast ausnahmslos in der Wassergeschwindigkeit 
begriindet. Der Stollendurchmesser hat demgegeniiber, wie aus Tabelle 1 deutlich hervorgeht, nur 
geringen Einfluf. 


1 Vgl. auch W. Grébner, L’Energia Elettrica (1939) Band 16 S. 595. 
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Tabelle 1. Warmeiibergangszahlen a9 fiir 0, = 7,5 °C. 


SS a 1 2 4 6 | 8 10 
i | 1952 3568 | 6522 9280 11 919 14 474 
24 | 1784 3261 5960 8480 10 892 1S) epee 
3 1693 3094 5654 8045 10 333 12 547 
5 | 1584. 2895 5290 1528 9 669 11 741 


Zur Bestimmung der Wassertemperaturen ist man meist auf Messungen an Anlagen, die sich in 
der Nahe des geplanten Stollens befinden, angewiesen. Fantoli setzt auf Grund finfjahriger Beob- 
achtungen im Alpengebiet Liro a Prestone die Temperatur des abflieBenden Wassers im Stollen 
als Funktion in ° C ein und bekommt 

60) = 7 +.4,5 sin i. =| =1—4,5 cos =, 
worin t die Zeit in Monaten bedeutet und t = 0 mit dem 15. Januar zusammenfiallt. 

In Abb. 1 ist die Funktion ((t) fiir den Zeitraum eines Jahres aufgetragen. Verschiebt man die 
Abszissenachse so, daB ein gleich groBer positiver und negativer Schwingungsausschlag entsteht, 
so deckt sich die Kurve mit einer Kosinusschwingung und berechtigt daher zu dem in der nach- 
stehendenRechnung gewahltenLésungsansatz. 

Das den Stollen durchflieBende Wasser wird 
dem Fels eine gewisse Anzahl Kalorien ent- 


4c 


senkrecht zur Bildebene 
| mit der Zeit veranderlich 


of ? of 
i er eee BOMOn AO 70 S77 7 Monae s ) V8 
75. Jan. 15. Jult 15, Jan. 
Abb. 1. Wassertemperaturen #(t) iiber ein Jahr aufgetragen. Abb. 2. Aufbau der Temperatur-Belastung. 


ziehen und mit zu Tal férdern. Demzufolge ergibt sich in der Umgebung des Stollens eine Volumen- 
verminderung der Felsmasse, was eine VergréBerung des Stollenradius nach sich zieht. Da der Fels 
aber in Wirklichkeit schichtig und unstetig ist, treten die thermischen Ausdehnungen ungleich- 
maBig auf. Darauf Riicksicht zu nehmen, wiirde aber groBe Schwierigkeiten in die Rechnung 
hineintragen, so daB entgegen der Wirklichkeit hier naherungsweise véllige Isotropie und Homo- 
genitat vorausgesetzt wird und daher die Gesetze der Elastizitatstheorie zugrundegelegt werden 
kénnen. Diese Annahme kann man um so mehr verantworten, als dadurch keine nennenswerten 
Fehler entstehen. 


2. Ermittlung der Temperaturverteilung. Um klarere Verhaltnisse fiir die Rechnung zu gewinnen, 
wird das Problem so aufgespalten, dafs erstens der Temperaturverlauf aus dem Unterschied zwischen 
der konstanten Gebirgs- und mittleren Wassertemperatur und zweitens der Anteil der periodisch 
schwankenden Wassertemperatur untersucht wird (Abb. 2). Beide Ergebnisse sind dann zum Schlufs 
zu uberlagern. 

Der Druckstollen kann als eine kreiszylindrische Durchbohrung des Gesteins mit dem Halb- 
messer 7; angesehen werden. Ist #(x, y, z, t,) die Felstemperatur im Punkt x, y, 2 zur Zeit ¢ unter der 
Voraussetzung, daB die eine Achse mit der des Stollens zusammenfallt, so geniigt diese Funktion 3 
der Warmegleichung 


AU IS NE tee 
Ops (ase Oy? | Oz 
Fiir den gesamten AuSenraum des Stollens ist die Temperaturleitfahigkeit 


A 
C=, 
Cy 
wobei c die spezifische Warme, y das Raumgewicht und A die Warmeleitfahigkeit des Gesteins 
bedeuten. Infolge der Symmetrie des Systems ist die Funktion unabhangig von x und in der y, z- 
Ebene unabhingig von y. Es kann daher die allgemeine Differentialgleichung der Warmeleitung 
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fiir den Fall achsensymmetrischer Belastung 


oo (#8 | 1 a8 ay, 
Ot a -a\Ortm aasaracr, 

angesetzt werden, die sich mit 9 = r/r; in der dimensionslosen Form 
ot rm (aoe e@ 00 


darstellt. Dabei erstreckt sich der Giiltigkeitsbereich der Funktion (0, t) auf das Feld r;< @ <<-- 00 


und 0< t<-+oo. Der EinfluB der periodischen Schwankung der Wassertemperatur auf das Tem- | 


peraturfeld wird bei einer Kreisfrequenz w dieser Schwankung erfaBt durch den Lésungsansatz 


Ho,t) =B Go (0 yes cos (wt — 9g) + Go (0 aa sin (w t —»)| , (2) 


i 
a 


woraus sich das Temperaturgefalle 


VP Le (oP) ema—o rene P)nio—a] 


bestimmen laBt. In diesen Ausdriicken bezeichnet G,(z) die Graysche Zylinderfunktion, die mit der 


Hankelschen Zylinderfunktion 
H(z) = I,(z) + iN,(2) 
durch die Gleichung 


G,(2) = = HY) 


verkniipft ist. In der Rechnung sind die Funktionen Gp,(r), Go9(r), G(r) und C,0(r) verwandt worden. 
Dabei weist der erste Index auf die Ordnung, der zweite auf den Real- oder Imaginarteil hin. Somit 
bedeutet 


Golr Vi) = Goulr) + i Geol) 
die Graysche Funktion nullter Ordnung mit dem komplexen Argument r Vis 
G,(r Ji) = G,,(r) + 1 G(r) 


die Graysche Funktion erster Ordnung vom gleichen Argument. Mit G* (r) und G*,(r) sind abkiirzend 
die negativen Ableitungen der Funktionen G,,(r) bzw. G,,(r) bezeichnet, die sich als Linearkombi- 
nationen der Grayschen Zylinderfunktionen dritter Art G,,(r) und G,,(r) darstellen lassen 


u(r) — Gy2(r dGo,(r “ Gi(r) + Giolr) dG),(r 
oe! 


und aus Tafeln | entnommen werden kénnen. 


Aus der Warmeiibergangsbedingung 


oo 

ie = (8 — Ow) (4) 
bzw. 

od t , 

oe ke (4’) 


in der } die Temperatur des Gebirges und #y die des Wassers im Stollen darstellt, folgt mit 9 = 0; 
und Ow = Ow max COS Wt 


TO ox rio io ro\ , rio 
- B Vy Ch by yl —B>C, to ip 22) + = Dimas 608 o| cos (w t — 9) 


\ 


rho rio ia io) ria aati 
ue E Bye Ge & / = ie Bop ie y22) — = Pymax sin o| sin(wt—qg)=0, (5) 


1 F. Télke, Besselsche und Hankelsche Zylinderfunktionen nullter bis dritter Ordnung, Stuttgart 1936. 
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ae = 
| Yen ol 22) te ( V2 
tgp = —— = (6) 
Ti w@ i@ tO i 
~ Yea 22) 98a 


ri 
— —— dymax sin p 


A a 0 
te Gh & V2) Ar a6 AG V ze 


errechnen. Damit ist das gesuchte Temperaturfeld 


. loa (o Vt EN eaerca gl 24 ants oom fem 


Pyne [aVe*)  ~/e2 on (af). 
Go2 \ Qi z ar cs / s Gh E \: a) 
gefunden. 


An dem Beispiele eines Stollendurchmessers von 3,0 m soll nun das entstehende Temperaturfeld 
erlautert werden. Hierbei mégen betragen 


woraus sich dann 


und 


die Warmeleitzahl A = 3,0 [keal/m h °C] 

die spezifische Warme ¢=0,20 [kcal/kg °C] 

die Warmeiibergangszahl « = 3100 [keal/m?h °C] 

ai K . fi am 22 — 200 ia 4 = =| 
ie Kreisfrequenz w ie eT 0.0073 LO0e* [bass 


Die HilfsgréBen ergeben sich zu 


(ses 

yee =—0,53899, “*—1550, 4 —6,452- 10-4, 

a ri 0 

26 aE ae a 
Ce\y---| = 0,7881, Gy,\) 2) = 0.0584,  Gt{V ==} = 1,6630, ¢*,\/— =] —0,3399, 
a a BI a Me a 

y = 2,4459 , sing = 0,64093 , cos pg = — 0,76760. 

Die weitere Rechnung wird tabellarisch durchgefiihrt und liefert im Endergebnis die bezogenen 


Temperaturen von Tabelle 2. 


Tabelle 2. Temperaturen Wt, a)/Pymax - 


{ 1,0 1,2 1,4 1,6 1,8 2,0 2,5 3,0 5,0 10,0 30,0 


+0,8658 |+0,7722 | +0,6996 | +0,6295 | +0,5653 | +0,5066 | +0,3808 | +0,2795 | +0,0529 | —0,0113 |—-0,0° 127 
+0,5000 |+0,4886 | +0,4754 | +0,4560 | +0,4334 | +0,4093 | +0,3453 | +0,2819 | +0,0925 | —0,0088 |—0,0° 260 
+0,0003 |+0,0739 | +0,1237 | +0,1601 | +0,1855 | +0,2033 | +0,2174 | +.0,2211 | +0,1070 | —0,0039 —0,0° 324 
—0,4995 |—0,3605 | —0,2610 | —0,1786 | —0,1121 | —0,0589 | +0,0312 | +0,0797 | +0,0930 | +0,0019 —0,0° 303 


—0,8656 —0,6983 | —0,5759 —0,4695 | —0,3797 | —0,3043 | —0,1637 | —0,0707 | +9,0540 | +0,0040 ES 199 


+0,9995 |+0,8491 | +0,7364 | +0,6345 | +0,5456 | +0,4682 | +0,3142 | +0,2022 | —0,0006 | —0,0107 emew 
| 
| 
| 
—0,9995 |—0,8491 | —0,7364 / —0,6345 —0,5456 | —0,4682 | —0,3142 | —0,2022 | +0,0006 | +0, 0107 | 


Da fiir die ganzmonatigen Zeiten t = 0 bis t = 5 und t = 6 bis t = 12 Monate cos (wt — ¢) und 
sin (w t — ¢) nur ihr Vorzeichen andern, so gilt 


) wee o eee nae 
== ) 0 = pee fea re (Sema faa Oy max /; —7 : 


In Abb. 3 sind diese Werte so aufgetragen, daB iiber der Abszisse @ jede Kurve den Temperatur- 
verlauf im Gestein wahrend eines bestimmten Monats anzeigt; die Abb. 4 ist dagegen so angelegt, 
daB die Abszisse die Zeitangabe tragt, und die Kurven je einem besonderen Kreisschnitt @ an- 
gehéren. Die Kurve fiir 9 = 1 ist mit der jahreszeitlichen Temperaturschwingung identisch. 
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3. Ermittlung der Temperaturspannungen!. Das vorliegende Problem stellt, wie auch das der 
Temperaturermittlung, eine unendliche Scheibe mit kreiszylindrischer Offnung dar. Die Scheibe 
unterliegt keinen Belastungen, sondern ist lediglich entsprechend der Temperaturverteilung einem 


Abb. 3. Temperaturverteilung in der Umgebung eines kreisférmigen Stollens (¢ als Parameter). 


o/ Vy max \| 
10 O~ Ze O } 
98 = ww Verlauf der Wassertemperatur | | 
06 = WN 
O4 oS | 


pe et NS 
VI. 7 


130 —$—~_ rinaes 
4 59 GYM OK. SBS Sen 5X. 59751. ON. 6H. 6, BY 150) 7X. 
— L - | A 
5. 


— 
77. 72 Monate j 


Abb. 4. Temperaturverteilung in der Umgebung eines kreisférmigen Stollens (9 als Parameter). 


achsensymmetrischen Spannungszustand ausgesetzt. Denkt man sich durch je zwei benachbarte 
Radial- und Zylinderschnitte ein Element aus der Scheibe herausgelést (Abb. 5), so wirken darauf | 
eine Radialspannung o, und eine Tangential- 
spannung o;. Schubspannungen sind infolge 
der achsensymmetrischen W armedehnung nicht 
miéglich. Somit sind die Normalspannungeno, | 
und o, gleichzeitig die Hauptspannungen. 

Die auf das Sektorelement wirkenden Krafte 
Gh ATAY verlaufen samtlich radial, da sich die beiden 
Tangentialkrafte auch zu einer radial gerich- 
teten Resultierenden zusammensetzen lassen. 


pray» 0 arn ag 


ip 


o,hAr 0;hAr 


4p Demzufolge lautet die auf die Radialrichtung 
bezogene Gleichgewichtsbedingung 
Abb. 5. Schnittelement einer Scheibe. a = h Ar Ag —o,h Ar Ap = 0 | 
oder nach Division durch h Ar Ap die wohlbekannte Gleichung | 
gene: 
hs dr 


1 Die Anregung zu dieser Rechnung verdanke ich Herrn F. Télke. 
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oder in dimensionsloser Form mit a oe 
Tj 
d(o o; 
O; = “om e (9) 


Sind ¢, und ¢, die Dehnungen der Scheibe in radialer und tangentialer Richtung und ¢, die 
Temperaturdehnungen, so gelten mit dem Elastizitatsmodul E und der Querdehnung y nach dem 
erweiterten Hookeschen Gesetz die Beziehungen 


1 1 
Er — £9 =F (Or — MO) » €1— €9 = = (0; —p0,). (10) 


Driickt man weiter ¢, und ¢, durch die Radialverschiebung u aus, so folgt in bekannter Weise - 


du u 
= oes gy) 
oder nach Einfiihrung von (11) in (10) 
du ih 1 
a 88 1 (Or — BM) s ~ =e +7 (—po,). (12) 


Mit (9) und (12) liegt ein vollstandiges System simultaner Differentialgleichungen vor, aus denen 
sich die drei Unbekannten u, o, und o, gewinnen lassen. Wird die zweite der Gleichungen (12) mit r 
multipliziert und differentiiert, so folgt 


d d 1 do d(ro, 
cerca ar cesta (13) 
Andererseits liefert die erste Gleichung (12) in Verbindung mit (9) 
du 1 d(r o,) 
eel aie od : (14) 
Wird (14) von (13) subtrahiert, so ergibt sich 
0 ae i (% br Gi—a,] 
oder in dimensionsloser Form 
de Lac, do o 
Osea ; f——*)., 1S 
aap alesaae fl (15) 
Multipliziert man diese Gleichung mit E und setzt o, aus (9) ein, so wird 
@(go,) , 1dgo) eo, , (Ee) __ l 
doe @ de go dg ee) 
Die homogene Lésung ist mit 
Cc 
(ooh = Coto (17) 


bekannt; die Stérungsfunktion dieser inhomogenen Differentialgleichung wird zunachst analytisch 
in der Form 
Beg = Eat =nP mit n= Hos (18) 


angegeben, wobei E den mittleren Elastizitétsmodul des Gesteins und a, den thermischen Aus- 
dehnungskoeffizienten bezeichnen. Fiir Beton betragt z.B. E ~ 2-10° [kg/cm?] und a, = 10~°, so 
daB sich dafiir der Wert n = 2 errechnet. Die Temperaturfunktion # ist dabei durch (2) dargestellt. 
Fiihrt man ihren Differentialquotienten (3) in (16) mit (18) ein, so ergibt sich 


be oa 
Beulineo) see // ra | ( \/ Ha) = 
ee ae mee Ge Ea, B oe G* le _] 08 (wt— 9) 


+Gh, E te) sin (wt — o) (19) 
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Die allgemeine Lésung dieser Differentialgleichung hat, wie man durch Einsetzen nachpriift, die 
Form 


: Sire | ao 
eo, =Ce 4 Heed ela ae sin Creeper 22) cos (or— 9). (20) 
i W 


6, 
— 030}- Abb. 6. Radialspannungen ———’——_infolge Temperatur im Gebirge um 
EX Oy max 
einen kreisférmigen Druckstollen (¢ als Parameter.) 


Pils ; Ss 
20 Abb. 7. Tangentialspannungen —_______ infolge Temperatur im Gebirge um einen 


i ; E ag oy max 
kreisférmigen Druckstollen (t als Parameter), 


Im folgenden sind nun aus den Randbedingungen die Integrationskonstanten zu bestimmen. 
Fiir den Innenrand g@ = 1 mufs co, = 0 sein. Damit geht (20) iiber in 


Ea,B rho n 
0=C,+CG+ ae ct (3) sin (wt — y) — G*, (/2) cos (os 9) : 
i 


i 
a 
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Aus der Bedingung, daB fiir @ = oo der Wert o, noch endlich sein mu, folet 


ay 


(21) 
Damit ist aber auch gleichzeitig die zweite Konstante bekannt 
Eds B\ ax re é Ro 
Cc, = — ize ler, (22) sin (wt —- p) — G* (2aks (wt — ») : 
a 


Diese in die durch g dividierte Gleichung (20) eingesetzt, ergibt die Radialspannung 


(22) 


oo Eds B 


| ees ne he CEN ie 
0, yee er GC <5 + v5 G* \o hs (wt—¢) 
a. 


XI. 
—$p— t 
ees 72 Monote 
0= 50 
- O51 
-70- 
=H 
-20 
6, 
Abb. 8. Tangentialspannungen 


Eas oy max 


infolge Temperatur im Gebirge um einen kreisférmigen Druckstollen (g als Parameter). 


Bildet man von (23) die Ableitung nach 9, so ergibt sich gemaB (9) die Tangentialspannung zu 


Tio acl ze ve ( 2) 
Oo, = Ea,B == Ci, ff a Ces OQ + — 22 Q a 
ate 


a 2 


Vee Miele, ro) 1 
— Ea, B — == G*, - ie > Gu Q 1%; C @ 
ane 


a 


In den Ausdriicken fiir o, und o, ist der Wert (7) von B einzusetzen. 
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Auch hier empfiehlt es sich, den Spannungsverlauf an einem Beispiele zu erértern. Unter den 


gleichen Annahmen wie bei der vorhergehenden Zahlenrechnung erhalt man 


Den weiteren Rechnungsgang liefert Tabelle 3 


1 = 0,973291 Owmax bzw. B/V/ P 


Erste Zeile o,/E.05 Op max > 


== 
ri Ww 


Tabelle 3. Radial- und Tangentialspannungen. 


zweite Zeile 0,/E 49 Oy max + 


JL S000 Twines 


1,0 1,2 1,4 1,6 1,8 2,0 2,5 3,0 5,0 10,0 30,. 

0 0,141 0,209 0,240 0,251 0,252 0,230 0,198 0,0934 0,0153 0,00 
1,857 ils y2lt 0,980 0,733 0,547 0,398 0,156 0,00704|—0,174 |—0,0485 |—0,00% 

0 0,125 0,189 OF2725 0,239 0,244 02235 0,213 0,121 0,0278 0,004 
1,610 1,206 0,948 0,755 0,607 0,487 0,270 0,123 |—0,126 |—0,0729 |—0,00 

0 0,0753 0,119 0,146 0,162 0,171 0,177 0,171 0,116 0,0329 0,00K 
0,931 0,769 0,661 0,576 0,504 0,441 0,312 0,206 |—0,0419 |—0,0774 |—0,00k 

0 0,00562) 0,0172 0,0299 0,0420 0,0525 0,0719 0,0827 0,0797 0,0291 0,00 
0,00320 0,125 0,198 0,242 0,267 0,277 0,270 0,234 0,0509 |—0,0614 |—0,00I 

0 —0,0656 |—0,0894 | —0,0940 |—0,0893 |—0,0803 |—0,0528 |—0,0276 |+0,0223 |+0,0175 | +0,00/ 
0,926 —0,552 |—0,319 |—0,157 |—0,042 |+0,126 |-+0,156 +0,199 |+0,131 |—0,0289 — 0,001 
0 —0,119 |—0,172 |—0,193 |—0,197 |—0,192 |—0,163 —0,131 |—0,0410 | +0,00125 +0,00) 
—1,606 —1,082 |—0,750 |—0,514 |—0,340 |—0,209 +0,000289) +0,111 |+0,177 +0,0113 |—0,00! 
0 —0,141 |—0,209 |—0,240 |—0,251 |—0,252 |—0,230 —0,198 |—0,0934 |—0,0153 —0,004 
—1],857 —1,321 |—0,980 |—0,733 |—0,547 |—0,398 |—0,156 —0,00704| +0,174 |-+0,0485 -+0,06) 


Die Eintragung wurde bis t = 6 Monate vorgenommen. Fir die eingeklammerten Monate ist zu 
beachten, daf sich die Vorzeichen der Spannungen umdrehen. Aus den Kurven von Abb. 6 ist zu 
ersehen, da die Radialspannungen o,/E a9 ymax vom Stollenrand beginnend sehr schnell zu- | 
nehmen und bei etwa 9 = 1,8 das Maximum des Maximums erreichen, um dann gegen null asympto- 
tisch abzuklingen. Zur Bestimmung der Tangentialspannung dient die Abb.7. Die GréBtwerte 
treten bei 9 = 1 auf und streben sehr schnell der Bezugsgeraden zu. Sehr aufschluBreich ist auch die 
Darstellung in Abb. 8, in der die Tangentialspannungen fiir verschiedene Kreisschnitte 9 als Funk- | 
tion von t aufgetragen sind. Die gestrichelte Kurve entspricht darin dem Verlauf der Wasser- | 


temperatur. 


(Eingegangen am 24. September 1954.) 


Anschrift des Verfassers: Prof. Dr.-Ing. K. Hirschfeld, Aachen, Technische Hochschule. 
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Uber die Berechnung der an rollenden Fahrzeugen 
wirkenden Haftreibungen. 


Von B. Stiickler. 


1. Einleitung. Die in der Dissertation des Verfassers} (im folgenden mit I bezeichnet) auf- 
gestellten Differentialgleichungen fiir die Bewegung eines idealisierten Krzftwagens gelten nur fiir 
den Fall, da8 alle Rader auf der Fahrebene rollen, ohne dabei zu gleiten. Das ist jedoch nur méglich, 
wenn hinreichend grofe Haftreibungskrafte auftreten. Dieses ist aber leider oft nicht der Fall, und 
dann wird es besonders gefahrlich, wenn das Fahrzeug auch noch ins Schleudern gerat. An die sich 
hieraus ergebenden Probleme ist man bisher vorwiegend experimentell herangegangen. Theore- 
tische Betrachtungen hielten sich in mehr oder weniger engem Rahmen und dienten mehr dazu, sich 
die beobachteten Vorginge verstandlich zu machen, als eine brauchbare Theorie aufzubauen. Zur 
Verkleinerung dieser Liicke méchte die vorliegende Arbeit beitragen. Und zwar wird auf theo- 
retischem Wege ermittelt, in welcher Weise die an den in I betrachteten Fahrzeugen wirkenden 
Haftreibungen von beliebigen eingeprigten Kraften und den Massenbeschleunigungen abhangen. 

An sich bestehen keine grundsatzlichen Schwierigkeiten, um diese Aufgabe zu lésen. Man 
kénnte beispielsweise elementar mit Hilfe des Schwerpunktsatzes und des allgemeinen Momenten- 
satzes vorgehen. Diese wiirden aber bei einer Anwendung auf das ganze Fahrzeug nur drei Gleichun- 
gen ergeben, in denen die gesuchten Haftreibungen vorkamen. Demgegeniiber treten an jedem Rad 
i. a. zwei Reibungskomponenten auf, so daB bei vier Radern acht Krafte zu bestimmen sind. 
Infolgedessen mifSten die genannten Satze auch noch auf Teile des Fahrzeuges angewandt werden. 
Dadurch kamen zwangslaufig innere Reaktionskrafte in die Rechnung hinein, die aber wieder 
eliminiert werden miiBten. Im ganzen gesehen ein jedenfalls sehr mihseliges und umstandliches 
Verfahren. 

Besser ware es schon, die Lagrangeschen Gleichungen zu benutzen ?. Aber auch dann ware ein 
umfangreiches Gleichungssystem zu lésen, und auBerdem wirde man die Ergebnisse in einer sehr 
uniibersichtlichen Form erhalten. 

Diese Nachteile fihrten zur Entwicklung des in Ziff. 2 angegebenen Verfahrens. Durch eine 
naheliegende Erganzung der sog. Lagrange-Eulerschen Gleichungen ? lassen sich die gesuchten 
Reaktionskrafte ohne Eliminations-ProzeB bestimmen, und zwar in einer bedeutend tbersicht- 
licheren Form als bei Benutzung der Lagrangeschen Gleichungen. Daher wird sich dieses Verfahren 
vor allem dann mit Vorteil anwenden lassen, wenn es sich um Systeme mit einer gréferen Zahl von 
Freiheitsgraden und Bedingungen — holonomen oder nichtholonomen — handelt. 


2. Erweiterung des Lagrange-Eulerschen Gleichungssystems. Dicse Gleichungen gehen aus dem 
Lagrangeschen Prinzip 


6A, = [ dK, or < [ b dr dm (1) 
hervor. Darin bedeutet 6.4, die virtuelle Arbeit der eingepragten Krafte, dK, die Summe der an 
einem Massenelement dm wirkenden eingepragten Krafte, b die Beschleunigung und Or eine vir- 
tuelle Verschiebung eines dm. Nun gehoren die betrachteten Fahrzeuge zu den Systemen mit end- 
lich vielen Freiheitsgraden. Fiir diese ist jeder Vektor r von einem festen Bezugspunkt zu einem 
dm des Systems eine Funktion der Lagrangeschen Koordinaten 

T= T (G2 Gore +o In)» 


so daf fiir die virtuellen Verschiebungen 
or 
6r = aq, 6q, (2) 
gilt. Bestehen nun kinematische Bedingungen von der Form 


> 54% =9 (irate a it Tt) (3) 


1 B. Stiickler, Ing.-Arch. 20 (1952), S. 337. 
2 Siehe z. B. G. Hamel, Theoretische Mechanik, S. 464, Berlin, Gottingen, Heidelberg 1949. 


3 Siehe G. Hamel, a. a. O. S. 480. 
19 
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so daB auch 
» b,, 04, = 0 (ip PAs 7) (4) 
gilt, so werden neue GeschwindigkeitsgréBen w, und virtuelle Verschiebungen ds, eingefiihrt durch 
= 0 fir w=1,2,...,m, 
= »b,,4 (5) 
Wy Pb heal fiir _ Soe 
= Osiris. at, 
= »'b 6 
bs, She | eee Om (6) 
wobei aber vorlaufig noch nicht beachtet werden darf, daB diese GréBen fiir yw = 1,..., m gleich 
Null sind. Es wird also gewissermafen zunachst der allgemeinere Fall betrachtet, daB das System 


noch nicht den Bedingungen (3) unterliegt. Die w,, sollen voneinander unabhangig sein, so dal es 
eindeutige Auflosungsformeln 


7b ee (7) 
ue 

6q, = Bee OS , (8) 
ub 


gibt. Damit bekommt man fiir dr nach (2) 


n 


ms es = or 
Tac a) Chine a Pa Bai. (9) 


ae is 


Mit diesen Beziehungen liefert 1 das Lagrange-Prinzip (1) 


al coon oT 
» K.— (a 0On (i | | ee Po 0) 
O3T 


pei 


mit 


T = T(q,,@,) als kinetische Energie, 


Sas oT oT 0 
K,= | aKa ( N= Bye 


Se acs = oq, 


= 

l 10 
A (10) 
7 

J 


ob 
= Qv 
Dube im Ds 24, (Be? Uy A= Dye eon) 3 
Daraus ergeben sich bei Bestehen der Bedingungen ds, = 0(u—=1,...,m) die Lagrange-Euler- 
schen Gleichungen 


ak OAE oT 
Kp bo, — (ae) + Dy Baer Pee (=m tse enon). (11) 
O37 


Um nun auch noch die Reaktionskrafte zu ermitteln, die bei Bedingungen w, = 0 auftreten, 
wird dieses Gleichungssystem in folgender Weise erweitert. Jede Bedingung (6) wird mit einem 
zunachst unbekannten Faktor multipliziert: 


A, 98, = 90 A ica ore trecdiE ye (12) 
Dann kann an Stelle von (10) geschrieben werden 


n 
q 


dmods oT 
= K, a Ay — (a OW (er a Dah Poor) 0s, UF (13) 
Q3T 


a 


Jetzt wird behauptet, daf alle ds,, als frei angesehen werden diirfen, so daB® neben (11) noch gilt 


donk oT 
=— K+ 5 poe — (Ge) + Dy Bucs Pee, Tse). a) 
Q,T 


Der Beweis ist hier besonders einfach; denn (14) kann einfach als Definition der J , angesehen werden. 
Damit ist jedoch noch nichts gewonnen. Denn aus dieser zunichst rein mathematischen Definition 


1 Siehe G. Hamel, a. a. O. S. 473 ff. 
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der A - geht in keiner Weise hervor, ob diese GréBen auch eine physikalische Bedeutung haben. Zur 
Klarung dieser Frage wird auf das Newtonsche Grundgesetz 


dK, + dKp =bdm (15) 


zuriickgegriffen, worin dKp die Summe der am dm wirkenden Reaktionskrafte bedeutet. Multi- 
pliziert man dieses skalar mit einem beliebigen Vektor, der zum Unterschied zu den virtuellen Ver- 
schiebungen mit dr bezeichnet werde, und summiert iiber das ganze System, so tritt an Stelle von (1) 


[dK,dr + [ dKgdr = [ bdrdm. (16) 
Wird nun speziell in Analogie zu (9) 


dr= 5) 4,ds,+ > a, OS, (17) 
w=1 w=m +1 
gewahlt, so liefert (16) offensichtlich als Erginzung von (10) 


m 


iss 9 donors oT 
aa K, a. it dKr a, — (i Bey (350) Si De Pye) DS, 
Q.t 


} 
| 
e=1 
( (18) 
Cal oT 
Ke. G oalerd, ee Pew.) ds, = 9. | 
0.7% 


Hierin sind aber im Gegensatz zu (10) alle Verschiebungen frei; folglich mu jede eckige Klammer 
Null ergeben. Damit ist zunichst bewiesen, daB fiir die durch (14) rein formal eingefiihrten A, gilt 

A= [eKee, (1,1). (19) 
Hierin treten an sich simtliche am System wirkende Reaktionskrafte auf. Man braucht aber jetzt 
nur noch an die Hamelsche Definition der Reaktionskrafte 1 zu denken, um die Frage nach der 
physikalischen Bedeutung der /,, in erschépfender Weise beantworten zu kénnen. Nach dieser ist 
eine Kraft dann eine Reaktionskraft, wenn ihr einziges Merkmal ? — das soll heiBen, das einzige, 
wodurch man auf das Vorhandensein der betreffenden Kraft schlieBen kann — eine Bewegungsein- 
schrankung ist. Nun ist das einzige Merkmal der Reaktionsgréfe 2,, gema ihrer Kinfiihrung (12) 
0s, = 9, was gleichbedeutend mit w, = 0 ist. Infolgedessen kann sich J,, nur aus solchen Reak- 
tionskraften zusammensetzen, die w, = 0 halten. Ferner ]a8t sich aus dem Vorzeichen von /,, 
etwas iiber die Richtung der betreffenden Reaktionskrafte aussagen; denn /,, Ds,, ist die Arbeit der 
Reaktionskrafte bei einer gedachten Verschiebung DS. Diese muB also bei DS, > 0 positiv sein, wenn 
1, >9, dagegen negativ, wenn J, < 0 ist. 


1 


u=m+1 


3. Anwendung auf rollende Fahrzeuge. In dieser Arbeit beschrankt sich die Anwendung auf die 
in I betrachteten Fahrzeuge, fiir die 1 = b’/] = 0 ist 3, wahrend der allgemeine Fall J ~ 0 in einer 
spateren Arbeit behandelt werden soll. 

Als erstes sind geeignete Verschiebungen zu wahlen. Da es naheliegt, die an jedem Rade wir- 
kende Haftreibung in eine Komponente tangential und in eine normal zur Bahn des Rades zu zer- 
legen, wird man zweckmaBigerweise auch Verschiebungen in diesen Richtungen wahlen. In Abb. 1 
bedeuten Ds, baw. Ds, Verschiebungen von B bzw. B, in Richtung der Vorder- bzw. Hinterradachse 
und Ds, bis Ds, Verschiebungen der Bertthrungspunkte der Rader tangential zu deren Bahnen. Dabei 
sollen die Pfeile angeben, in welchen Richtungen die Verschiebungen positiv gezahlt werden. Als 
virtuelle Verschiebungen werden die schon in I benutzten GréBen ds und dy gewahlt, und es wird 


Os, = Os , ds, = 0x (20) 


gesetzt. Entsprechendes gelte fiir die GeschwindigkeitsgréBen w,,, fiir die ja nach (5) und (6) die 
gleichen Koeffizienten mabgebend sein sollen wie fiir die ds,,. Diese b,,, sind nun zu bestimmen. In 
Abb. 2 sei vp die Geschwindigkeit von B bei einer beliebigen Bewegung des Fahrzeuges auf der 
Fahrebene. Dieser Vektor ist in verschiedene Komponenten zerlegt: 

a) in @, und w,, also in Richtung der Vorderradachse und senkrecht dazu, 

b) in vj und v,, das soll hei®en, in Richtung der Hinterradachse und senkrecht dazu, 

c) in x und y. 

1 Siehe G. Hamel, a. a. O. S. 517. 

2 Zum Merkmal-Begriff s. G. Hamel, a. a. O. S. 510 u. S. 43. 

31, 5S. 338. 2 

Gye 
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Fiir diese Komponenten gilt 
w,=—sind x+cosdy, w, =cosPx+sindy, (21) 
vj =— sind, x + cost, y, v, = cos 9,4 + sin By. (22) 
Andrerseits ist v; = W) + 1d, , da die Relativbewegung von B gegeniiber B, nur eine Drehung 
um B, mit der Winkelgeschwindigkeit by sein kann. Daraus folgt 
W, = — sin 9, x + cos 0, y—19,. (23) 
SchlieBlich ergibt sich aus I (9) 
Wy = —cd—1@, Os i Os aig: | (24) 
Wy, =O, +cO—Tq, Wg =v, + 64, — Qs - | 
Wird nun fiir die benutzten wahren Koordinaten die Reihenfolge 


h=%*%, B=Y¥. BHA UWH—x%> b—P> WH Per WM Is —Ps (25) 


Abb. 1. Die gewihlten Verschiebungen. Abb, 2. 


gewahlt, so fiihren die Beziehungen (21) bis (24) zu der Koeffizienten-Matrix 


—sin? cos? 0 0 0 0 0 0 
—sin?, cost, —I 0 0 0 0 0 
cos}? sing —c —ce —r 0 0 0 
cos? sin } Cc c Q =F 0 0 
b — 
(us) cos#, sind, —a 0 0 0 —r, O |’ ee) 
cos?, sin #, b 0 0 0 0 —r, 
cos sin } 0 0 0 0 0 0 
{ 0 0 0 1 0 0 0 0) 
deren Reziproke lautet 
(— sin & 0 0 0 0 0 cos? 0 } 
cos 0 0 0 0 0 sin?’ 0 
il 1 1 
7 008 x i 0 0 0 0 7 sin x 0 
0 0 0 0 0 0 0 1 
ee Sie ls J oe Q 
Coa ey Thea Pes ae rate a sin) — Fl, ay 
1 
= 00s x —< Ear 0 0 + (1 +-+sinz) c 
1 5 i: 
— 7, (70082 + sin) ai 0 0 ee 0 7, (0082 — 7 sin x) 0 
1 
ab. 7 b 1 1 
=a; eos 7, — sin z)— 0 0 0 — z, (008 + sin 2) 0 | 
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Hiermit werden die £9, errechnet. Da diese in (14) mit w, zu multiplizieren sind, von denen aber 
beim reinen Rollen nur w, = v und w, = ¥ ~ 0 sind, kommen zur Ermittlung der ReaktionsgréBen 
nach (14) nur die Buo7 und Cos in Frage, und zwar fiir u = 1,...,m—=—6. Die Rechnung ergibt, 
dafi nur folgende Koeffizienten von Null verschieden sind: 


Buy = 2 CO8iX.3 Boz = — = ’ Biss 

Bi = = ? Boog = = COS X » Piss $= 1. 

Bis Pos = 2 sin ¥ , (tee | : (28) 
Biaz ( = - sin x , Boer = 2 sin x » 

Beas J 


Als nachstes ist die kinetische Energie T, des Fahrzeuges zu bestimmen und daraus die Impuls- 
komponenten P,,. Bei der Bildung von T, darf zwar noch nicht von den Bedingungen WO, =9 
Gebrauch gemacht werden, weil man sonst auch P|, = 0 bekame, wohl aber diirfen quadratische 
Glieder mit zwei verschwindenden w,, fortgelassen werden. So tritt jetzt an Stelle von I (25) 


2 Tp = mpv? + Ji + 2mp (sv, + 1,0,) 0, + C, (v2? + $2). 
Dazu kommen vom Vorderradsystem die Glieder 
2mav? +248? + C(p? +92), 
wobei hier 2A das Haupttragheitsmoment des ganzen Vorderradsystems um die Lenkachse ist und 


2 mp dessen Masse, deren Schwerpunkt auf der Lenkachse angenommen wird. Da nach Abb. 2 
Vv, = — sin ¥ W, + cos ¥ Ww, ist, so bekommt man fiir das ganze Fahrzeug 


2 T) = mp (— sin yw, + cos yw)? Cent yet 0A) 8: + 29, (— mrs siny w, + mpl, | 


ee eee (29) 
+ mps cosy, +24Aa,) +2mgut+240t+ C(P++G(e +9). J 
a Catt cae 
Hieraus folgt fiir P, = Ba T(q,.@,,): 
OT, OT, OD, OT) Oh | OT, My, OTy Hy | OT, Os rr 
oe CO : 08, Ou P= Op 80, a OP, Oy 891 Gory, ° Ops By ” eo 


Die hier ttbergangene Rechnung liefert unter Benutzung der Abkiirzungen! I (21), I (28), I (31) und 


1 fia OF ls 
Ms =p(J+24 ote re C,)—me > Gp) 
die Werte 
je 2 ed Be eras w | 
i= (mr cos X + > sin 4 \ Wy i 1 Ws » 
: DeAd 
P, = — (mp, cos x + ms sin ¥) W, — —— Ws» 
Py=— 7 [Lf sing)er ce), 
C Coe 
P, 4 ie eu) (2 1. TE sin 7) -- za 9: (32) 
Py=— {Ft (0084 — sing) on 
ry 
by: 
Pe=— Fh (0084 +p sin x) on 
ry 
ae 2A’, 
P= (rm, + mea sin2 z+ masin® 2) oy + pe 
J 


2 
1 In I (21) muB es 4’ = A + = C heifen. 
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Da leicht einzusehen ist, daB fiir 1 =1,...,7 (@T/ds,,) =0 ist, so liefert (14) nach Einsetzen von 
(28) zunachst 


= Le L fooe y Py 4 P, desing (Py Bet] oy 2 (Pee ee 
[ 


~~ Oy l 

dy By Paap rt Cee Pea Fea ie) ei Ore (33) 
: > 

A, = — K, + P,, 

Ag=—K,+ P,, 

As = — Kye 

Ag =— K, ar P, 


Interessant ist die Bedeutung der eckigen Klammern in A, und/,. Fir diese erhalt man nach kurzer 
Rechnung 


° I; . - 
cosy P, +'P, + (P,+ Pi P,) siny =m —- sinyv, (34) 


* 2 . . 
P, + cosy Py— (Ps + Po) sing = —m (1-7) cos x + sing sin yYv, (35) 


worin /J,« und s* die Koordinaten des Schwerpunktes des ganzen Fahrzeuges sind (Abb. 3). Wie sich 
aber leicht nachweisen laBt, ist (34) bzw. (35) nichts anderes als die Komponente des Impulses mv,» 
in W,- bzw. w,-Richtung. 

Nun ist noch zu klaren, aus welchen Reaktionskraften sich die A,, zusammensetzen. Da nach 
Ziff. 2/,,0s, die Arbeit der Reaktionskrafte bei einer Verschiebung Ds, ist und /,, nur solche Reak- 
tionskrafte enthalt, die 
w, =9 halten, ergibt sich 
folgendes. 

Von allen am Fahrzeug 
wirkendenReaktionskraften 
kénnen nur die von der 
Fahrebene auf die Vorder- 
rader ausgetibten Reibungs- 
krafte ein seitliches Gleiten 
dieser Rader verhindern. 
Also tragen nur diese Krafte 
zuw, = 0 bei. Da aber von 
diesen bei einer Verschie- 
bung Ds, nur die Kompo- 
nenten senkrecht zu den 
Radern Arbeit leisten, ist 2, 
gleich deren Summe, die mit 
Ry bezeichnet wird. Ent- 
sprechend ist A, die Summe 
der Komponenten der an 
den Hinterradern angreifen- 
den Reibungskrafte senk- 
recht zu den Radern, fiir 
die Ry gesetzt wird. Dar- 
aus ersieht man, dak die 
Theorie der Systeme starrer 
Korper nichts tiber die Gré- 
Be der genannten Kompo- 
nenten im Ejinzelnen aus- 
sagen kann. Dagegen sind 


Abb. 3. Bezeichnungen und positive Richtungen der Komponenten der Haftreibungen pe Une : 
2 ee tad Carnie ® s die tibrigen /,, nur Einzel- 


krafte, denn wenn z.B. nur 
Ds; ~ 0 gewahlt wird, was einer Drehung des linken Vorderrades um dessen Achse entspricht, so 
leistet nur die an diesem Rade wirkende Reibung Arbeit, und zwar nur deren Komponente R,; in 
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Radrichtung (Abb. 3). Entsprechendes gilt fiir 2, bis J, so daB man schlieBlich aus (32) und (33) 
fiir die genannten Reaktionskrafte erhalt 


— 


A = Ry K, + (mry 0082.4 + Stsin2 2) + mS sin ¥ v? 


2 A“ ae 2 
-- 88 XK + (m +m + Se — mr sin 2 7 — my sin? 7) sae 
1 


A, = Ruy = — K, — (mp, cos y + mg sin y) 0 + > (1p }e08 2 + sing] sin y v? 


2 A’ 
TO ea ala 
(36) 
A, = R,. =— K,— = 7 a Sa ea | 
A,=R,,=— Ko at ae sin x) +ei], 
A;= Ri = —K,— Gee at hate 
J 


C, d , 
A= eniesc = (cos x + 7 Tsing). 
Zu diesem System gehéren fiir 4 = 7,8 die Bewegungsgleichungen I (78) und I (79). 


4, Die Komponenten K_,, der wichtigsten eingeprigten Krifte. a) Schwerkraft. Diese liefert 
nur dann einen Beitrag zu den K,,, wenn der durch den Schwerpunkt des Fahrzeuges gehende resul- 
tierende Schwerkraftvektor Komponenten senkrecht zur Normale der Fahrebene hat. Diese seien 
mit G, und G, bezeichnet (Abb. 3). Die zugehérigen K,, ergeben sich aus der Arbeit dieser Krafte 
gemaB 

=8 
DA, = K psa oa; 1 0s, (37) 
p=m+1 
Zunachst wird nur Ds, ~ 0 gewahlt, was einer Drehung des ganzen Fahrzeuges um M, entspricht, 
und zwar um den Winkel 


ye -~ bs, == os cos y Ds, . (38) 


Dann 1aB8t sich an Hand der Abb. 3 ee nachweisen, da dabei von der Schwerkraft die Arbeit 

DA, = [G, lw + G (Q sin y — s*)] Dd, (39) 

geleistet wird, woraus sich nach (38) die erste der in (40) zusammengefaBten Beziehungen ergibt. 

K, bekommt man, indem nur Ds, # 0 gewahlt wird, was nur durch Drehung des ganzen Fahrzeuges 

um B méglich ist. Und schlieBlich leistet die Schwerkraft noch Arbeit bei einer Verschiebung ds,, 
also bei einer Drehung um M. Die Ergebnisse lauten 

K, = (6.7 —G, 7) cosy + G,siny, | 

( 


ee Q-F)+G5 = 


(40) 
) (6. — 6,7) sin y —G, cosy. | 
b) Luftwiderstand. Aus I, Abb. 4, ist ersichtlich, daB in (40) nur P; cos statt G,, 
P;sin « statt G,, 0 statt 1,* und e statt s* zu setzen ist, um zu erhalten 


Kn Pr cos (sin x — 7087), | 


k= Le [sina + 5 cos «) ( (41) 


K, = — P; cos « (cos + sina). | 
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c) Momente um die Radachsen, wie Antriebsmomente, Brems- 
momente, Momente der Rollreibung. Fiir diese gilt zundchst 


pA, = M10 + M,, Do + Mir de, + Mar 093 - (42) 


Hieraus bekommt man durch Umrechnung auf die Ds,, mit Hilfe der Matrix (27) 


il : 
kK, = - (M,,— M,,) + _ (6M,,— aM] CE Oi (M,, + Mri) siny ; 
is 1 

; 1 

eae E (M,,— My) + = (6M — aM], 
1 
Mhr 
gee ed, Kees, Ree K, 2 22s (43) 
* if r ry ry 


1 1 F 
K,= : (Mi els M,1) ais ae 08 ae je Mer Mut) at rl (6 M),,—aMyj,i) s1n X > 


Cc 


K,=“ (M,,—M,)). 


5. Uber den Einflu8 der Schwerpunktslage auf die wichtigsten Massenbeschleunigungsglieder von 
Ry und Ry. Es sind dies die Glieder mit ¢ und v? in (36), da normalerweise ¥ und ¥ klein sind und 
in der reinen Geradeausfahrt und der reinen Kreisfahrt ganz verschwinden. Die Glieder mit v? 
lassen sich verhaltnismaBig leicht tibersehen; sie riihren von der Normalbeschleunigung des Fahr- 
zeuges her, denn nach Abb. 3 ist (1/l) sin y v? = v?/o, das ist aber die Normalbeschleunigung der 
Lenkachse. Ihre vektorielle Summe muB nach dem Schwerpunktssatz in der reinen Kreisfahrt mit 
konstanter Geschwindigkeit gleich m b,« sein, wenn alle K,, = 0 sind. Komplizierter ist der Aufbau 
der Koeffizienten von ?. Fiir die in ihnen auftretenden Konstanten gilt nach I (28), I (31), I (21) + 
und (31) 


m s* C, b—a ) 
Ua Nessie fica ilar: Sr peak! 
yh 2A Pogo Cee eo a ane 
nm = Pp (J ot. — MF {- 2 r ( aS a ) r) ( (44) 
1 Destine I. | 
my =p (I+ 2AF2CS + A+) — me FZ, | 


In der Praxis sind aber gewohnlich die Tragheitsmomente der Rader klein gegeniiber dem 
Gesamttragheitsmoment des Fahrzeuges 7. Beachtet man ferner, da nach dem Steinerschen Satz 


J+i2A+2mpglP = Jy« + m (1,3 + s*?) (45) 
gilt, so bekommt man naherungsweise fiir den Fall, daB s*/l1< 1 ist, 
Tes ie Jx Is Is 
Mz — —™m Lae : ms © 5 m ii Fl (46) 


Um ein ungefahres Bild der Verhaltnisse zu erhalten, wird folgendes Zahlenbeispiel benutzt 2: 
m = 130 kgm-*sec?.. J,, = 230 komsec?, 1 — 2.6, ae 


Diese Werte ergeben fiir die Koeffizienten in den genannten Gliedern die Diagramme der Abb. 4 
und 5. Zu diesen ist folgendes zu bemerken: In Abb. 4 bzw. 5 sind die Koeffizienten von v2 bzw. 
v in Abhangigkeit von x dargestellt. Als Kurvenparameter ist das Verhaltnis 1,./I gewahlt. Wo 
nicht anderes angegeben, ist s* = 0. Von dieser Schwerpunktskoordinate hangt (Ry), tiberhaupt 
nicht und (R,), nur in geringem Mabe ab, da in der Praxis s*/l1< 1. Beachtlicher ist jedoch der Ein- 
fluB von s* auf die Glieder mit 6. Es ist nun besonders im Hinblick auf die Gefabr des Schleuderns 
wichtig, daf} die fiir das Rollen erforderlichen seitlichen Haftreibungskrafte Ry und Ry moglichst 
klein bleiben. An sich kénnten auftretende Massenbeschleunigungen durch eingepragte Krafte 
wenigstens teilweise kompensiert werden. Z. B. konnen bekanntlich bei einer reinen Kreisfahrt mit 
konstanter Geschwindigkeit Ry und Ry = 0 sein, wenn die Fahrbahn eine entsprechende Neigung 
aufweist. Auch innere eingeprigte Krafte kénnen nach Ziff. 4, c) Ry und Ry verkleinern. Aber i.a 


1 Siehe FuBnote von S. 283. 
2 Vel. Zahlenbeispiel I, S. 348. 
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wird es sehr schwierig sein, durch eingepragte Krafte den jeweils erforderlichen Ausgleich zu 
schaffen. Man denke beispielsweise nur daran, da® fiir den Stadtverkehr die Kurven gewohnlich 
nicht iiberhéht werden kénnen. Daher sollte man auch darauf bedacht sein, die Massenbeschleuni- 
gungsglieder so klein wie méglich zu halten oder zumindest eine unginstige Verteilung auf Vorder- 
und Hinterrader zu vermeiden. So ersieht man aus Abb. 5, daB® im Zahlenbeispiel (Ry); nahezu 


(hy), 
kgmi“Isek? 
30 


20 


10 


Abb, 4. 


Abb. 5. 


verschwindet, wenn s* = 0 und I,./] = 1/2. Der letzte Wert ist zwar fiir (Ry) ; ziemlich ungiinstig, 
es ist aber folgendes zu bedenken: Beim Bremsen wird infolge der Zunahme des von der Fahrbahn 
auf die Vorderrader ausgeiibten Normaldruckes der Héchstwert der Haftreibung an den Vorder- 
radern vergroBert, dagegen an den Hinterraddern verkleinert. Und von diesem Gesichtspunkt aus 
betrachtet, kénnte sich |,«/] = 1/2 durchaus als giinstiger Wert erweisen. 


(Eingegangen am 4. Oktober 1954.) 
Anschrift des Verfassers: Dr.-Ing. B. Stiickler, Berlin-Steglitz, Mariendorfer Str. 32b. 
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Zur Beulung des flachen kreiszylindrischen Schalenstreifens 
bei beliebiger orthogonaler Anisotropie. 


Von F. Driickler. 


1. Aufgabenstellung. Stabilitatsuntersuchungen diinnwandiger anisotroper Rechteckplatten, 
bei denen die elastischen Hauptachsen nicht in die Richtung der Langs- und Querrander fallen, sind 
bereits zahlreich vorhanden. Fiir fast alle praktisch vorkommenden Falle der Belastung und Lage- 
rung haben Green und Hearmon! die kritischen Lasten berechnet und ihre Ergebnisse mit denen 
anderer Autoren verglichen. Fiir den wichtigen Fall der Schubbelastung eines Plattenstreifens 
hat Miiller-Magyari? ebenfalls die kritische Schubspannung berechnet, wobei er einen anderen Weg 


als Green und Hearmon benutzt hat. Alle bisherigen Untersuchungen beziehen sich jedoch nur auf | 


den ,,ebenen“ Plattenstreifen. Die vorliegende Arbeit? befaBt sich dariiber hinaus mit der Berech- 
nung der Beulspannungen ,,gekriimmter“ Plattenstreifen unter Druck- und Schubbelastung bei 
Neigung der elastischen Hauptachsen. Unter einem gekriimmten Plattenstreifen, fir den in dieser 
Arbeit der vielleicht etwas deutlichere Ausdruck ,,flacher kreiszylindrischer Schalenstreifen“ 
benutzt sei, ist ein durch zwei Achsenschnitte aus einer unendlich langen diinnen Kreiszylinder- 
schale herausgeschnittener Streifen zu verstehen, dessen Breite klein gegenttber dem Krimmungs- 
radius ist (Abb. 1). 

Berechnungen fiir isotrope zylindrische Schalen haben unter anderem Leggett 4, Kromm ® 
und Marquerre ® vorgenommen. Wahrend sich Leggett und Marquerre von vornherein mit der flachen 
Zylinderschale befassen, hat Kromm mit Hilfe des vollstandigen Ansatzes die Giltigkeitsgrenze der 
Ergebnisse bei flacher Schale fiir die verschiedenen Lastfalle nachgewiesen. Man erhalt bei der 
flachen Schale ein System von zwei partiellen Differentialgleichungen, im allgemeinen Fall ein 
System von drei partiellen Differentialgleichungen. Dazu kommt noch, da sich bei der flachen 
Zylinderschale die sogenannten Lastglieder in dem System der Differentialgleichungen in wesent- 
lich giinstigerer Form ergeben, was bei der Auswertung der Beuldeterminante von Bedeutung ist. 
Fur die im folgenden zu behandelnden geneigten Hauptachsen wird gegeniiber dem isotropen Fall 
das Elastizitatsgesetz und damit das Endergebnis allerdings wesentlich komplizierter. Die errech- 
neten Formeln gestatten auch die Ermittlung der Beulspannungen bei versteiften Zylinderschalen 
aus Stahl und anderen isotropen Baustoffen, bei denen die Anisotropie durch orthogenale Steifen- 
kreuze, die nicht in Richtung der Langs- und Querrander fallen, zustandekommt. Die Schwer- 
achsen der Versteifungen miissen allerdings in der Mittelflache liegen. 


2. Bezeichnungen. Im folgenden sollen 
bedeuten 

b die Plattenbreite, 

1 die Halbwellenlange der Verschiebungs- 
gréBen in x-Richtung, die beim un- 
endlich langen Plattenstreifen peri- 
odisch verlaufen miissen, 

t die Plattendicke, 

u,v, w die Verschiebungen eines Punktes 
der Mittelflache, Atha ee ; te aaled 

a By; 7 Hie Verzerrungen Ace Mittelflache, ae acher zylindrischer peo bei beliebiger orthogonaler 

Ex3 €y3y die Verzerrungen im Abstand z von der Mittelflache, 

Ny = ot die gleichmaBig verteilte Druckbelastung je Langeneinheit, 

T, = tt die gleichmaBig verteilte Schubbelastung je Langeneinheit, 

Nx, Ny, Nxyy Mx, My, Mxy die Schnittreaktionen am Plattenelement, 

Ci Sj,  Elastizitatskonstanten. 


1 A, E. Green und R. F. S. Hearmon, Phil. Mag., (7) 36 (1945) S. 650. 

® F. Miiller-Magyari, Osterr. Ing.-Arch., 4 (1950), S. 22. 

3 Gekiirzte Fassung einer von der Fakultat fiir Bauwesen der Technischen Hochschule Hannover genehmig- 
ten Dissertation, Referent: Prof. Dr.-Ing. habil. A. Pfliiger, Korreferent: Prof. Dr.-Ing. W. Zerna. : 

* D. M. A. Leggett, Proc. Roy. Soc. Lond. (A) 162 (1937) S. 62. 

° A. Kromm, Luftfahrtforschung 15 (1938), S. 517. 

° Marquerre, Neuere Festigkeitsprobleme des Ingenieurs S. 229. Berlin 1950. 
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3. Beziehungen zwischen Schnittkraften und Verschiebungen. Die Schnittkrafte der flachen 


Zylinderschale sind } 
Ne on a, Neve) Oy :02 Ne Ns Uz, | 1 
M,=fo.zdz, M,=fo,zdz, Miy=Myz=[rzdz. | ©) 


Fir den Fall der geneigten elastischen Hauptachsen muf das verallgemeinerte Hookesche Gesetz ? 
angesetzt werden: 


On = Cyy Ex TF Cyp Ey + gy » 
Oy = Cyy Ex 1 Cyn Ey TH Cog | (2) 
T = yg Ex + yg Ey + C337 


oder aufgelést nach den Verzerrungen 


Ex = $4, Ox + Sy9 Cy + S43 T , 


Ey = Sig Ox “1 Seg Dy 1 Sog T > | (3) 


(eee Roe yr a OG) Say 21 = 


Dabei gilt im Fall orthogonaler Anisotropie fiir die c;, und s;, folgendes: Betrachtet man zunachst 
den Sonderfall eines Koordinatensystems <, ¥, dessen Achsen mit den elastischen Hauptachsen 
zusammenfallen, so bestehen fiir die zugeh6rigen ¢;, und 5,;, die Beziehungen 


ee Me erent ef gs Pe 
ji Lp? @ 1 — vz v7’ 33 2 12 1l— vz v5 1 — vz v5” 
(4) 
5 = gt ~ 1 m 1 x ve vy 
cel jae 2D = See BEN oer) SS aie, ae 
Be Es G Ee re 


Ist nun x, y ein rechtwinkliges Koordinatensystem, das nicht mit den elastischen Hauptachsen 
zusammenfallt, und « der Winkel zwischen x und x, dann driicken sich die c;, durch die ¢,, folgen- 
dermafen aus: 


Cy, = €,, costa + ©, sint a + (4 €,, + 2 €,,) cos® « sin? a , 

C33 = Cop Co8* & + C,, sint x + (4 Cyg + 2 C,,) cos? & - sin? a , 

C33 = Cg3 cost x + yg sint a + (6,1 + Cog — 2 Cg — 2 Cyy) cos? a sin? x , 
C19 = G4, cost a + G, sin? a + (C1 + Co — 4 Coq) cos? « sin? a , 


C13 = (C4 — 2 S33 — 19) cos? w sin « + (— Cy + 2 Cag + Gy) cosa sin? a , 
Cog = (Cop + 2 Cag + C19) cos? x sin a + (C4, — 2 C33 — C49) cos a sin? a . 
Entsprechendes gilt fiir die s,;,: 
$41 = &1 c08t & + Sy sin* « + (S35 + 2 8,5) cos? a sin? a , 
Sop = Soy C084 & + 54, sin*t ~ + (833 + 2 81.) cos? a sin? a , 
= 55, cost a + 85, sint ~ + (4 84, + 4 “5 — 2 833 — 8 8,5) cos « sin? «, 
Sip = Sq co84 & + yy sint % + (811 + Sop — 853) cos” & sin? a , 
Sig = (2 81 — S53 — 2 349) cos? w sin & + (— 2 S49 + S33 + 2 819) cos a sin? a , 


S93 = (— 2 Sop — Sg3 + 2 8,9) cos? x sina + (2 °11 — Sg3 — 2 849) cos x sin? x . 


wH 
oc) 
© 
| 
| 


| 


Der Zusammenhang zwischen Verzerrungen und Verschiebungen ergibt sich wie bei der ebenen 
Platte zu 


= ee ay ie 
Ey = Ey — ZW, Ey = Ey — SW 9 y= y—2zw 9 (7) 


worin die iiberstrichenen Groen die Verzerrungen der Mittelflache bedeuten. Fiir die flache 


Zylinderschale gilt 
Ex =U , éyaut—, ee ee ee (8) 


1 Siehe FuBnoten 4 und 5 von Seite 288. 
2 W. Geckeler, Elastizitatstheorie anisotroper Kérper, Handbuch der Physik, Band 6, S. 104; Berlin 1928; 


A. E. H. Love und A. Timpe: Lehrbuch der Elastizitat, Deutsche Ubersetzung, Leipzig und Berlin 1907; siehe 
auch Miiller-Magyari, FuBnote 2 von Seite 288 
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Im Vergleich zur Platte kommt fiir é, ein Zusatzglied hinzu, das durch die VergréBerung des Kreises 
mit dem Halbmesser r infolge w entsteht. Man erhalt also 


&, =u —zw’, ey =v + —— aw, y=u' + v' —2zw". (9) 
Setzt man (2) und (9) in (1) ein, so erhalt man das Elastizitatsgesetz fiir die Schnittkrafte 


N,=t lew w + C5 (o" af 7) + C3 (u* + “) > 


L 


Ny = les w + Co9 (» aie = + C93 (u’ + “) 2 


L 


Nxy = his wu’ co” ae 2 + C33 (u* + v’) ’ 


(10) 


3 


t ” a Te 
M,=—7lai¥ + ¢,w" + 2¢,w''], 

pe yr ee le 
My =— 75 le12¥ “+ Cy W"* + 2 Cy, w'*] , 

e i eo le 
Miy = — 75 le3¥ + ¢,,w +2 cs, W *] - 


4. Aufstellung der Beulbedingungen. Die Liésung des vorliegenden Stablilitatsproblems erfolgt 
mit Hilfe der Energiemethode !. Dazu werden die Formanderungsarbeiten der inneren und auBeren 
Krafte benétigt. Die Forminderungsarbeit der inneren Krafte setzt sich zusammen aus der Form- 
anderungsarbeit der Zusatzspannungen, die in der Schale infolge der Verschiebungen u, v und w 
entstehen und aus der Formanderungsarbeit, welche die Grundspannungen bei den zusatzlichen 
Verzerrungen der Schale leisten. Fiir die Formanderungsarbeit der Zusatzspannungen bezogen auf 
die Flacheneinheit gilt 


1 
ty = 7 | (nee + oy ey +07) de. (11) 
Mit (2) und (7) erhalt man daraus die gesamte Formanderungsarbeit der Zusatzspannungen zu 
+4/2 b 21 
Ai Sucy [ox Ex — (Cy ex + Cyey + Gygy) ZW" + Oy Ey — (Ca Ex + Cop &y + lpg) ZW" 
ype WY 0) a 
+ © Y — (CygExz + Cog by + Cygy) 220 *| dx dy dz. (12) 


Die Formanderungsarbeit der Krafte des Grundzustandes ist beim Vorhandensein von Normal- und 
Schubkraften gegeben zu 

a;. = —tt(1 + é,) (1+ €,) siny, Gig — 0 Vee (13) 
Hierbei miissen sin y und é, bis auf quadratische Glieder berechnet werden. Sie ergeben sich genau 
wie bei der Platte zu 
ut+o’t ww Sey 1 12 
(Cama ae am) 


Damit wird die Formanderungsarbeit der Krafte des Grundzustandes 


Sin4. 
b 21 b 21 1 
Aig + Aig——rt{ [(w +0! tu! wi) dedy tar] | (w + yw) de dy. (15) 
00 OF 0 


Die Formanderungsarbeit der a4uBeren Krafte lautet fiir den Fall des unendlich langen Schalen- 
streifens 


21 b b 
A,=t tf (Uyp— Uy) dx +7 tf (Vs=21— Vx=o) dy —o tf (Ux 2)— Ux=o) dy . (16) 

0 0 0 
Mit Beriicksichtigung von (8) und (9) und Integration iiber der Schalenstirke ergibt sich nun die 


gesamte Formanderungsarbeit, die gleich der negativen ersten und zweiten Variation des Potentials 
des Grundzustandes ist. 


* Eine zusammenfassende Darstellung der Energiemethode fiir Stabilitatsprobleme findet sich z. B. bei 
A. Pfliiger, Stabilitatsprobleme der Elastostatik, Berlin 1950. 
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Wegen Gleichgewicht des Grundzustandes wird die erste Variation gleich Null, was aus (15) und 
(16) zu ersehen ist, indem sich die linearen Anteile herausheben. Fiir die zweite Variation, die mit 


a5 6°29 bezeichnet sei, gilt dann 


l 


b 2 
Pg tt 
0 0 


o.u + dy [r + al +t (v +u°) + - ee? + Cog w? + Acsgw'*? + 2¢,.0 w" 
+ 4, Ww" w’* + 4c, ww’) + 27 w’ w'— ow] dx dy. (17) 


Fir das zu lésende Variationsproblem 6(622,) = 0 ist hier das Verfahren nach Galerkin benutzt 
worden. Fiihrt man in (17) die Variation 6(62,) durch, so erhalt man fiir die mit u und v behafte- 
ten Glieder die folgenden Ausdriicke: 
b 21 b 21a. 
R, Se: Ou’ eee ee) S* bu dx dy, 


00 


Be est ost | eens = bu dx dy. 
0 0 0 00 


0 


b 21 b 21a | 
a tJ foydv'dedy =f 0,80 dx —e] | 2 bude dy, | 
00 0 0 
b 21 | 
R, = —| [ 0, dw dx dy, ' (18) 
00 | 
b 21 
Beet ff x de" dedy = tr» | dye | bv dx dy, | 
0 
| 
| 


Die unter den Doppelintegralen verbleibenden Ausdriicke sind fiir verschwindende du, dv die 
Eulerschen Differentialgleichungen. Man kann sie durch Einfiihrung einer Spannungsfunktion 


identisch befriedigen: @” =6,, Crag), Nye ae (19) 
Die mit w behafteten Glieder ergeben Ausdriicke entsprechend (18). Bei Anwendung des Galerkin- 


Verfahrens hat man also nur fiir ® und w, und nicht fiir u, v, w, Ansatzfunktionen zu wahlen. 
Hinsichtlich der Randbedingungen gilt dabei folgendes: In x-Richtung soll die Schale unendlich 
lang sein (Abb. 1). Es wird daher in dieser Richtung periodischer Verlauf der u, v, w und ihrer 
Ableitungen verlangt. An den Langsrandern wird momentenfreie Auflagerung, Verschwinden der 
Tangential- und Radialverschiebung und das Verschwinden der Normalspannung gefordert, d. h. 
berry — 05 y= b: u=0, w= 05° Ny=0, Vive 0%. (20) 
Zur Naherungslésung nach Galerkin werden folgende Ansatze gemacht: 


i pee i Bee b 
P=Uye, w=WMyje*, P=z (21) 
mit Uy) = A, sin, Wy) => B,, sin. (22) 


Die Ansatze befriedigen nur angenahert die geforderten Randbedingungen. Sie sind trotzdem 
brauchbar, da die Ausdriicke fiir M, und u zeigen, dafs nur die Glieder an den Randern nicht ver- 
schwinden, die aus der Neigung der elastischen Hauptachsen herriihren. Nach Variation unter dem 
Integralzeichen in (17) — wie bereits fiir u und v gezeigt — erhalt man mit den Ansatzen (21) die 
folgenden Gleichungen, in denen nun auch die Randausdriicke, welche infolge der nur angenaherten 
Befriedigung der Randbedingungen nicht verschwinden, enthalten sind: 

b 


; eoce a 4 m7 27 Tee mes 3 U° 9 uw. U 
[ Juv + So GaP U— (251. + 833) 73 B U + 2 85 F518 Fs S135 UB 
0 capes 
a t Tu. ; a 
+ 5p BW 5 [aU dy — 545 5 ip U Aven = 10", 


*_ pt 4 7 pe mt Sy ET ery | 6 
hela W + C59 We"? — 2 (2 33 + C12) Ga B W—4 e375 18 + 4 95 5 iB 


6 b 


==). 


Uo ep W205 ip W|OWdy +27, 53 ip WOW" 


—) b? 


| 
| 
| 
; 0 2 Pa(23)) 


0 
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Wie noch weiter unten gezeigt wird, ist jedoch der Einflu8 der Randausdriicke auf das Endergebnis 
von untergeordneter Bedeutung. Mit den Ansatzen (22) und ihren Ableitungen erhalt man unter 
Beriicksichtigung der Integrationsformeln 


b 2n 
b — —_—_. fiir n-- m ungerade 
(a ee n? — m? : : 


0 : 0 fir n+ m gerade , 


ia funny Srl (24) 


b b 
J sin sin ™* J cos "= cos =" = fii nat 
A le Wie 7 == The 


b 
C06 ee cose | = 
b baie 


i 2 fiir n-+_m ungerade , 
0 fiir n-+- m gerade 


folgendes Gleichungssystem fiir die Bestimmung der Freiwerte A und B: 


> th 3) n-m =e 
Andy +i— S) Am 75 em + Bab =0, (25a) 


m2 


/ a6 : / 
—A,e + B, (d,—o') —i— >| Bn =a (fm + v') = 9, (25b) 
mit den Abkiirzungen 


1 
By = Sy, 24 + So 4 + (2 So + Sg) 1” B? , Cm = Sy B® + > S138 (n? + m?) , 
2 2 62 j b2 
be aps CO Se ON a (26) 


d, = - [¢11 B* + Coy m4 + 2 (2eg3 + Cy) 2? 87], fn = 5 [2 cy3 B® + C93 B (n? + m?*)]. 


Hierbei durchlaufen n und m gerade und ungerade Zahlenfolgen, n-- m ungerade. Das Gleichungs- 
system ist homogen, es hat nur dann eine von Null verschiedene Lésung, wenn die Koeffizienten- 
determinante verschwindet. Die Entwicklung dieser sogenannten Beuldeterminante liefert die 
gesuchten kritischen Lasten. Aus dem Gleichungssystem (25a) kann man fiir n = 1, 2,...,% die 
A durch B ausdriicken und in das Gleichungssystem (25b) einsetzen. Man erhalt so eine Koeffi- 
zientendeterminante i-ter Ordnung. 


5. Auswertung fiir reine Schubkraft. Streicht man in (25a) und (25b) die Glieder mit den Kon- 
stanten C5, 5,3, und C55, $3, die nur bei geneigten Hauptachsen auftreten und setzt fiir eine Neigung 
von a = 0 fiir s;, ¢;, die Werte (4) fiir Isotropie ein, so entkoppeln sich die beiden Systeme (25a) 
und (25b). Man erhalt ein Gleichungssystem mit den Konstanten B,, B,,, dessen Determinante 
dieselbe ist, die Kromm! als Ergebnis seiner Berechnung erhalt. Bei Green und Hearmon? findet man 
fiir die unendlich lange Platte unter Schubbelastung bei geneigten Hauptachsen als endgiiltige Glei- 
chung dieselbe wie hier in (25b), wenn man die Kriimmungsglieder streicht. Numerische Rechnun- 
gen sind allerdings nur fiir « = 0 und « = 90° zu finden unter Verwendung der zweiten Naherung. 
(Die ersten neun Elemente der unendlichen Determinante werden betrachtet.) Angaben iiber die 
Differenz der ersten und zweiten Naherung sind nicht vorhanden. Fir denselben Fall gibt Miiller- 
Magyari® Werte an, die auf anderem Wege gefunden worden sind. Es ist deshalb dieser Lastfall 
naher untersucht worden. Man erhalt aus (25b) fiir den ebenen Fall: 

Erste Naherung: 


3 \2 1 
[2 ¢13 8? + Sey5 + T*]? = (73) B? nb aI “22 BE + 2 (eyo + 2 ex) lenb* a on 55 + 8 (ey. +2 ex) 
mit t= 7 ee (27) 


Zweite Naherung: 


1 
B ab 1 Cop B t 2 (Cy, + 2 =) ip nb ae cons + 8 (ey, 4+ 2 ex) 


81 
x en BP + Cp pe + 18 (cy. + 2 ex) a (5a) B? (2 e436? + 13 ey. + 7) (28) 


16 \2 81 
aa (53) Pa leab™ T Cy9 Be + 18 (cy. + 2 =) [2 e158? + 5 ey, + 7*]?. 


1 Siehe FuBnote 5 von Seite 288. 2 Siehe FuBnote 1 von Seite 288. 3 Siehe FuBnote 2 von Seite 288. 


XXIII. Band 1955 Driickler: Zur Beulung des flachen kreiszylindrischen Schalenstreifens usw. 293 


Numerische Rechnungen sind wie bei M iiller-Magyari fiir Sperrholz (Buche, drei Lagen mit 
E, = 120 000 kg/cm? , E,, = 60 000 kg/cm? , », = vy ~0, G=10 000 kg/cm’, geltend fiir eine 
Gesamtwandstarke t = 0,8 — 1,5 mm) durchgefiihrt worden. Man erhailt 
fiir «= 0°, erste Naherung: 7*/40 000 = 10,05, 
fiir « = 90°, erste Naherung: 1*/40 000 = 7,11, ; 
fiir « = 45°, erste Naherung: 177/40 000 — — 8,52, 75/40 000 = 12,83 , 

zweite Naherung: 17/40 000 = — 7,67, 13/40 000 = 12,59. 


Im folgenden sind diese Werte, mit denen von Miiller-Magyari und Seydel+ verglichen. Als 
BezugsgréBe fiir die Kurve 


nach Miiller-Magyari gilt 


we 4 if 
x = 5 Ven cb ? 
wobei 
ao TU, 12 b2 
LS —$—$ 


2 24 3 
me" t Vey ¢8, 


ist; 


die Werte der vorliegenden 
Arbeit sind in einfacher 
Weise umgeformt worden: 


t*/40 000 = x 2 Die er- 


sten Naherungen liefern fiir 
« = 0unda = 90° gute Er- 
gebnisse. Da fiir « = 45° ge- 
naueVergleichswerte fehlen, 
sind hierfiir die erste und 
zweite Naherung berechnet 
worden; sie differieren um 


a . 
rund 10 Yo Es bleibt zu Abb. 2. Kritische Schubspannungen eines unendlich langen Plattenstreifens bei beliebiger ortho- 
erwarten, daB der Unter- gonaler Anisotropie nach F. Miiller-Magyari. Beispiel: Sperrholz, E (langs) = 120 000 kg/cm, 


schied von der zweiten zur —— Werte EET ee ened Sees Seydel, ote des Verfassers. 
dritten Naherung ziemlich 
bedeutungslos ist. Auf die numerische Auswertung der dritten Naherung ist deshalb verzichtet 
worden, zumal der Aufwand dafiir jenseits des Ertraglichen liegen diirfte. 

Der geringe EinfluB der Randausdriicke (23) auf die Ergebnisse sei fiir die erste Naherung fiir 
& = 45° nachgewiesen. Ohne die Randausdriicke verandern sich die Abkiirzungen c,, und f,, in 


(26) zu 
Cn (P* 5,, pom eq) 5 


ie 0 Cn oe BT Cs) a) (9) 
Damit erhalt man fiir den ebenen Fall die erste Naherung 
[2 (¢13, B® + ¢o3) + 7*] [2 (crs 6? + 4 3) + 7°] 
= (75) ig en ie toa + (4 e353 + 2 es) ln B? + Co ze + 4 (4es5 + 2 es)}. (30) 


Die numerische Rechnung fiir « = 45° ergibt mit 
t{/40 000 = — 8,58 , 73/40 000 = 12,89 


gegentber den oben angegebenen genauen Werten nur eine bedeutungslose Abweichung. 

Gibt man sich mit den am Sonderfall der Platte gewonnenen Erkenntnissen zufrieden und 
nimmt die Differenz zwischen der ersten und zweiten Naherung in Hinblick auf die ungenaue Er- 
fiillung der theoretisch vorausgesetzten Eigenschaften anisotroper Baustoffe auch beim Schalen- 
streifen als noch tragbar hin, so erhalt man in diesem Fall fiir die kritische Beulspannung durch reine 
Schubbelastung aus der Auflésung der Systeme (25a) und (25b) fiir n =1, 2; m=1,2; nm 


ungerade 


1 FE, Seydel, Ing.-Arch. 4 (1933) S. 169. 
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E C15 B? + Sey, + T* 
| ow B4 (B? S93 + 2,9 513) |? 
T 16\2 
[S11 + Sao B4 + (2 Sy +593)B?] [16531 + Sop f* + 4 (2 $12 + $33) B?] — (=) B? (B? S93 + 2,9 S13)” | 


2 1 
= (53) len pre “22 Be + 2 (2 ¢33 + Cp) 
w B® [16 sy + See B* + (2 Sie 4 $33) 487] 
16\? . 
[S11 + See 8* + (2 Sig + $33) B?] [16 sy, + Soo B* + 4 (2 Sig + S23) B?] — (==) B? (B? S93 + 2,9 S13) | 


16 
oe 


(31) 


a. 


x leu B? + Opp a + B (2 e535 + C49) 


| 


w B?-[s11 + $2284 + (2 Sie + $33) B?] \p? 
| A 16\? 3 
[sir + So B* + (2512 + $33) B?] [16 Su, + Soo B4* + 4 (2 Sip + S33) B?] — (5) B? (B? Seg + 2,95 S13)” 
: a 12 b? 12 b4 
Be eae ae = Fae 


Man sieht, daB fiir den flachen kreiszylindrischen Schalenstreifen praktisch nur die erste Naherung 


zur numerischen Berechnung in Frage kommt. 
Numerische Rechnungen sind fiir Sperrholz (mit denselben Elastizitatskonstanten wie beim 


ebenen Fall) bei einer Hauptachsenneigung von « = 0 und « = 90° durchgefiihrt worden. Die 
kritischen Schubspannungen 1*/40000 und die dazugehérigen 6?- Werte sind fiir verschiedene Kriim- 


mungen in Abhangigkeit von yo aufgetragen. 


x*/40000 ie 


60 


70 
60 
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Abb. 3. Kritische Schubspannungen eines flachen zylindrischen Schalen- 
streifens bei beliebiger orthogonaler Anisotropie fiir ~« = 0 unda = 90°, 


ee 
in Abhangigkeit von |/w . 


E (langs) = 120 000 kg/em, E (quer) = 60 000 kg/em, 
G (langs-quer) = 10 000 kg/cm. 


== 
0 1s RENT SE 6 ean 6 mee OO) 
Abb. 4. Die zu den kritischen Schubspannungen eines flachen zylindrischen Schalenstreifens bei beliebiger orthogonaler Anisotropie fiir 


a= 0 und « = 90° gehdrenden bezogenen Halbwellenlingen ? in Abhangigkeit von |’ . 


5. Zusammenfassung. Es wird die Stabilitatsgrenze eines flachen kreiszylindrischen Schalen- 
streifens unter Druck- und Schubkraften bei beliebiger orthogonaler Anisotropie untersucht. An 
den Langsrandern ist er angenahert momentenfrei auf biegeweichen Langstragern gelagert. Die 
Ermittlung der kritischen Last geschieht mit Hilfe der Energiemethode unter Benutzung des 
Galerkinschen Verfahrens, wobei fiir die Verschiebungen Fourier-Reihen angesetzt werden. Man 
erhalt als Ergebnis der Rechnung eine unendliche Beuldeterminante. Diese wird fiir den ebenen 


Fall (Plattenstreifen) fiir Neigung der elastischen Hauptachsen ausgewertet und mit Ergebnissen 
aus bereits bekannten Untersuchungen verglichen. 


(Eingegangen am 6. Oktober 1954.) 
Anschrift des Verfassers: Dr.-Ing. Friedrich Driickler, Bremen, Gropelinger HeerstraBe 399. 
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W andnahe Geschwindigkeitsverteilung turbulenter 
Grenzschichtstroémungen mit Druckanstieg. 


Von W. Szablewski. 


1. Einleitung. Vernachlassigt man in der Bewegungsgleichung einer turbulenten Grenzschicht- 
strémung fiir eine wandnahe Zone die Triagheitsglieder, so erhailt man durch Integration die 
Gleichung 
dy 


Ty ay (1) 


D 
x 
(t) Wandschubspannung, sonstige Bezeichnungen siehe Abschn. 2). Aus dieser Gleichung ist un- 
mittelbar ersichtlich, daf} bei einer Grenzschichtstr6mung gegen wachsenden Druckanstieg, der 
schlieBlich die Ablésung (t) = 0) erzwingt, der Druckgradient immer weiter zur Wand hin bestim- 
mend in das Gleichgewicht der Krafte eingreift und nicht mehr gegeniiber den Reibungskraften 
vernachlassigbar ist. Fir das bekannte logarithmische Geschwindigkeitsgesetz, das aus der An- 
nahme konstanter Schubspannung (t ~ T)) gewonnen wird, bedeutet das, da®B seine Existenzschicht 
mit wachsendem Druckanstieg immer weiter eingeengt und schlieBlich ausgeléscht wird [fiir die 
Ablésungsstelle kann nach (1) von einem logarithmischen Gesetz offenbar nicht mehr die Rede sein]. 
In einer kiirzlich erschienenen Arbeit 1 habe ich auf der Grundlage des Prandtlschen Schubspan- 
nungsansatzes ein dem Druckgradienten Rechnung tragendes Wandgesetz als Erweiterung des 
logarithmischen Gesetzes abgeleitet, was voraussetzte, daB eine Schicht mit logarithmischen Gesetz 
zumindest noch rudimentiar vorhanden ist. Es blieben also Bereiche starkeren Druckanstiegs und 
insbesondere die Ablésung ausgeschlossen. In der vorliegenden Abhandlung lassen wir die genannte 
Einschrankung fallen und beschreiben die wandnahe Geschwindigkeitsverteilung turbulenter 
Grenzschichtstr6mungen mit Druckanstieg in voller Allgemeinheit gemaf (1). Wir erfassen also 
insbesondere auch die Ablésung. 


2. Approximative Berechnung eines elliptischen Integrals. Im folgenden bedeuten: x die Bogen- 
lange der Profilkontur, y den Wandabstand, 6 die Grenzschichtdicke, u, v die Geschwindigkeits- 
komponenten in x, y-Richtung, @ die Dichte, u die Zahigkeit, y = y/o die kinematische Zahigkeit, 
p den statischen Druck, t die Schubspannung. 

In Gleichung (1) konstituiert sich die Schubspannung aus einem laminaren und turbulenten 
Anteil 

T = Tham > Trurb 


mit Tiam = [4 0u/dy und — bei Zugrundelegung des Prandtlschen Ansatzes — Tru,» = 0 [ (du/dy)? 
mit dem Mischungsweg I. Fiir | setzen wir in Wandnahe unter Beriicksichtigung der laminaren 
Unterschicht 


wo 6, die Dicke der laminaren Unterschicht ist. 


Fuhren wir die folgenden dimensionslosen GroBen ein: 


we ye webs ato By ee CPN 
ree ga age age Nd Ae 
wo u, die Geschwindigkeit der Potentialstrémung am Rande der Grenzschicht ist, so transformiert 
sich (1) in 
Dey (UP Cpe a A 3 
Res dn | (3) Fale: g- (3) 


Mit (2) 1aBt sich (3) integrieren und ergibt, wenn wir noch 


A A A : rt) 
a) (4— 9 no) 5 (Nn — No) mit no = a 


1 W. Szablewski, Ing.-Arch. 22 (1954), S. 268. 
20 


296 Szablewski: Wandnahe Geschwindigkeitsverteilung usw. Ingenieur-Archiv 


beachten, 


n—No ; | i 
1 oe le + 4 x0? Res (7 — 1%)” (« =a n) 2 (4 no 1 


tome i 
ite yom — - —— a = — d 
WV) 2 «2 Res / (7 — No)? (7 — No) 
+ Re; (, =e S dy "No (4) 


fiir 7 — <1. Das letztere Glied gibt dabei die Geschwindigkeit am Rande der laminaren Unter- 
schicht an. 


Fiir nicht zu kleine Wandabstinde erhalt man, da Re; normalerweise eine sehr grobe 


Zahl ist, 
1—Ne ee j) 
1 Ve > m) 2 (7) — %) | 
FS - —— — d 
er a) rane e (1 — No) | 
fiir (5) 
7 poeenen SSS KH SL | 
i | 
2x Re, / cy — 3 %o | 


Hier 1a8t sich die Quadratur mittels elementarer Funktionen durchfiihren und ergibt 
2 Aa ee es ara eH ie | 
P= (61 — $10) 3 1") Sioa) | 
no) A A weer | 
io Tne 0 V(o 9 m) iy (1 — Mo) Vo “9 Mo (6) 
2 In — ———_— — — 
: i j | 


wobei die Integrationskonstante jedoch wegen der in (5) bei 7 — % = 0 auftretenden Singularitat 
nicht bestimmbar ist. Die in (5) angegebene Begrenzung des Wandabstandes nach unten ist physi- 
kalisch mit der Vernachlassigung der laminaren Reibung gegeniiber der turbulenten Reibung 
gleichbedeutend: man erhalt (6) auch sofort aus (3), wenn man dort die laminare Reibung vernach- 
lassigt. Die obige Gleichung gibt also die funktionale Gestalt an, die fiir 7— 7) ~ 7 an Stelle des 


logarithmischen Gesetzes zu treten hatte! Fir y)<7 und — > M0< c, erhalt man die bereits von 


uns friiher abgeleitete Funktion !; die zweite Bedingung ist jedoch im Bereich der Ablésung offen- 
bar nicht erfillbar. 


Zur Berechnung der Integrationskonstanten ist die Berechnung des vollstandigen Integrals (4) 
nicht zu umgehen. Wir fiihren diese im folgenden approximativ durch, wobei sich der Fehler als 
praktisch belanglos ergeben wird. 


Fir die Berechnung erweist es sich zunachst als zweckmafig, das Integral (4) auf die Variablen 


V hy rae 
Be eee 2 Ver 1011 


Vs i y 
fs 5 “le 


zu transformieren®, (Auf die physikalische Bedeutung dieser Transformation kommen wir in 


1 W. Szablewski, a. a. O., Gleichung (6): 
ue ae ae 

2, A a=\\. Patel ee RCE ee: , 

Pi (W/o > 7] V ai} =| Mes In — hares — — + konst. 


a = 
\/« cry n+Ve 


care Trin led 
® Hierbei bedeutet V, = |/—° + — = 6 > Elk : 
ierbei bedeutet V,, ‘ =F odes vel. Gl. (24) 


| 
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Abschn. 3 naher zu sprechen.) Man erhalt dann bei voriibergehender Einfiihrung der Abkiirzungen! 


oe ae = Ve uf V,. 89 
soe ee % v v 
fiir (4) die einfachere Gestalt 
11 /V14+422(0 + Fx) —1 V3 F V,6 
yaaa CAFO te 4 Moh(1 Fah), () 
0 
_ wobei der Parameter 
has — 2/2 


Ee 2) 
Re, (« Sara m) 


technisch eine positive kleine Zahl (< 1) ist (vgl. Abschn. 4). Es handelt sich also um ein elliptisches 
Integral. Die Berechnung kénnte im Einzelfall mittels der Theorie der elliptischen Integrale oder 
auch einfach durch graphische Integration erfolgen. Nachfolgend werden wir in sehr guter Approxi- 
mation eine einfache Beschreibung mittels elementarer Funktionen geben. 

a) Wir nehmen folgende identische Umformung des Integrals (7) vor: 


Ll fjle4a4e—i V6 FV,6 
a al a dx | aff) — ; * ) 
0 


v vy 


aes Jl 442220 + F x)—J14+ 42 x2 
z dx 
Dy Reh x 
0 


| 
| 
| @ 
| 


Das erste Integral laBt sich elementar berechnen und ergibt einschlieBlich der obigen Konstante 


111—/14 42 2 1 ——_—. V0 FV,,6 
a + —In (2x x +71 + 40 x) 4 2415 oo 


(9) 


x 
In der Asymptote ergibt diese Funktion eine logarithmische GesetzmaBigkeit °. 


Im zweiten Integral setzen wir approximativ 


1+ 4220 +Fx)—J1 +402 2 Sp ees 


x2 x 


(10) 


Physikalisch hei®t das, daB bei der Berechnung des Einflusses des Druckgradienten die lami- 
nare Reibung auBer Acht bleibt. Das ist also weniger als die in (5) vollzogene Vernachlassi- 


Vises 
—— 


gung der laminaren Reibung tberhaupt. Das Integral dx weist nunmehr im 


Gegensatz zu (5) keine Singularitat bei 7 — 1 = 0 bzw. x = 0 auf. Die Integration ergibt 


4° Fx 1 


2 ae 
%, = — (Yl Fe —1)4 Oe Se 1 (11) 
Insgesamt erhalten wir also nach (9) und (11) die approximative Lésung 
a ee ho en VL ae 
2 x? x " i 
eae ae 1, 41/14+Fe—1, V,6 FV,6 
+ —(/1 + Fx—1) | _ in I fs Ss ses A Oerueral | 


b) Wir fiihren nunmehr eine Fehlerabschitzung durch. Die Fehlerfunktion lautet nach (10) 


1 1| “V1 + 40 21+ Fx) dean | E* dx | ee FF sta {¢ | > 03) 
zm | 2 x eee «|| 
0 


0 ; 
1 Eine Verwechslung mit den eingangs eingefiihrten Bezeichnungen x Bogenlange, y Wandabstand dirfte 
wohl ausgeschlossen sein. 
2 Vel. W. Szablewski, Z. angew. Math. Mech. 31 (1951), S. 309. 


20* 
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Wir stellen fest: 
1. Die Fehlerfunktion (13) ist eme negative, monoton fallende Funktion. Das folgt aus der 


Bemerkung, dab a i= 
(jl ta—fa)—(y1+b —yb)<0 fir a>b>0 


ist. Der negative Maximalfehler wird demnach fiir x —> oo erreicht. Wir beschranken uns im folgen- 
den auf die Ermittlung einer Schranke fiir den asymptotischen Fehler. 

2. Die drei letzten Integrale der Fehlerfunktion (13) sind exakt auswertbar und ergeben einen 
negativen Beitrag. Das erste (elliptische) Integral 


= 
0 


x? (1 + Fx) 
5 dx , 


x 


das einen positiven Beitrag liefert, wird numerisch abgeschatzt. Dabei machen wir Gebrauch davon, 


daB yl + é fir 0 <& <1 in eine gleichmahig konvergente alternierende Reihe entwickelbar ist. 
Infolge des alternierenden Charakters liefert das jeweilig letzte Glied der Reihe eine Fehlerschranke. 
Unter Festlegung von x, durch 


Ax? x! (1 + Fx,) =1 (14) 
1aBt sich die Wurzel y1 + 4 x2 x2(1 + Fx) fur die Intervalle 0 <* <x, und x >x, in folgender 
Weise durch eine Reihe darstellen: 


1) Fir x <x, kommt 


1 (42F x3 1 /4:2F x3 \? 
Vises i224) 8 (ee ane 


4x2 F x3 
L+4 x2 x? 


yl + 4x? x2 (1+ Fx) =jl +4 x? x? 


Nach (14) ist 
eA Mee 


Wir beriicksichtigen nur das erste Glied der Reihe. Infolge des alternierenden Charakters der Reihe 
bleiben wir damit auf der sicheren Seite der Abschatzung; denn wir nehmen damit einen zu kleinen 
positiven Beitrag zur Berechnung des negativen Fehlers. Wir erhalten dann insgesamt fiir x < x, 
folgende Schranke der Fehlerfunktion 


ion —————— 41 Jlrs 1 
C7) 14+ Fx,—1)—2x1 * 
2) Fir x > x, lautet die Reihenentwicklung 
eared 1 1 1 1 2 
dex Fa/ | rar esa] ae 8 Peer TF5| a: a], 
Wir beriicksichtigen wieder nur das erste Glied und erhalten unter Beriicksichtigung von 
1 
uAeh 
LENE Woe Tee, 
nach (14) fiir x > x, die Fehlerschranke 
it == ao 
sg l2x(L—y2+Fx,)—2%m242x%m(1+4/2+Fx,)]. (16) 


3. Insgesamt haben wir so fiir den (maximalen) asymptotischen Fehler nach (15) und (16) die 
Schranke 


1 =| —————— Awl Sy Tesi oad. 
|= 4x (fore, — 1) — 2 * 
2 x2 NAV erat 
(17) 
+ 2x(1—y2+4+Fa,)—2xIn2+2xn(1+/2+4Fx,) 
gefunden. 


Diesen Ausdruck schatzen wir weiter mittels Entwicklung in alternierende Reihen ab, die wir 
jeweils im Sinne der Gewinnung einer Schranke abbrechen. ; 
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. ah 1 oe) aa. eee a 
Mit! —y14Fx, ==// +7F x), Se aare = halt +4Fs,) 
(zu grofe Werte der negativen Beitraige) und 
In (1 +/2 + F x,) 
1 


=In(1+y72)4 — 2a ig ee , : 5) 
O+12)+ Saag waa eta) 


- 1 
Fx, Jl+Fx,+1 2 Bp 


(zu kleine Werte der positiven Beitrage) erhalten wir die Schranke 


Ll oe 1+) 2 rs 1 (4+/2 
2x Pll ya + In 5 = 7 j2f*,— 7. Gtr -4- 7) (Fx,)*| . (18) 
Setzen wir hier nach (14) Co see 
x 


(zu groBer Wert), so folgt schlieBlich mit ES eel! 
3 252 
fiir den negativen Maximalfehler die Schranke 
1 eel 1 /F\2 

— 5, [0225 Y oye x 1 cae (1?) 
Diese Schranke ist natiirlich verbesserungsfahig. Beriicksichtigt man in 1) und 2) je drei Glieder 
der Reihenentwicklung, was einen betrichtlich gréBeren Rechenaufwand erfordert, so erhalt man 
in analoger Weise (wir geben nur das Resultat an) die wesentlich scharfere Schranke 


— 2 foo + HE=1E 4 o( 2). 20 


% 


Der Betrag der von F freien Konstante ist nunmehr stark herabgedriickt worden; exakt mu er 
gleich Null sein, wie aus der Integraldarstellung (13) der Fehlerfunktion folgt. Der sich nach (20) 
ergebende Fehler ist fiir technische Strémungen sehr klein (vgl. Abschn. 4, c). 


3. Wandnahe Geschwindigkeitsverteilung turbulenter Grenzschichten mit Druckanstieg. a) Be- 
nutzen wir wieder unsere Ausgangskoordinaten, so haben wir nach (12) folgendes Ergebnis erhalten: 


ib the a 1 A 
P= 38 Re, 7— MN ay + 4x? Res (4—F m)(1—m? 


Fi 
Yo ao: "lo ] D) A oat 9 2 A 2 
= n | 2 xReg ce, — > No(n —M) + |/ 1 + 4x? Reg @ ery no) (7 — No) 


+ 2/3) — Fw) a9 | Oe) 


fiir 17 — <1, d.h. Tragheitsglieder vernachlassighar und | ~ x (y — 6) - 

Das ist die Geschwindigkeitsverteilung, die an die laminare Unterschicht anschlieBt. Sie erfabt 
insbesondere die Ubergangsschicht von gleicher GréBenordnung der laminaren und turbulenten 
Reibung. 

In der vollturbulenten Schicht nimmt die Verteilung (21) fiir nicht zu kleine Wandabstande die 
einfachere Gestalt an 


1 Den Abschatzungen wird Fx, < 1 bzw. nach (14) F < 2 x vorausgesetzt. 
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2 A A ey) 
9. = ee (Va— 7h |= "9 hte Va— oy ) 


Fie, sn Aa te eh i 
¥ In| 4 7 ; 5 : 
G 2 (0 rae mw) Pe gad oF \/« ei) 


: (22) 


A 
(1 — > "No en ao 
Ya-3 
t 9 5 1+ In (4 és Res |/ = + Res (.— zy i "o 
i 
2x Re/e— 5m <n<l. 


"No 
Das ist dieselbe funktionale Form, die wir bereits kiirzlich unter einer gewissen Einschrankung 
(siehe Ziff. 1) abgeleitet hatten 1. Es tritt nur an die Stelle des Parameters c, bemerkenswerterweise 


der Parameter — Ny - Verteilung (22) stellt die Verallgemeinerung des logarith- 
mischen Gesetzes auf Grenzschichtstroémungen mit Druckanstieg 
dar: fiir schwachen Druckanstieg //2 ~ 0 erhalt man wieder das logarithmische Gesetz. Wir wollen 
die durch (22) beschriebene GesetzmaBigkeit im folgenden kurz als ,,Wandgesetz bezeichnen. 

Wir geben noch die Verteilung in der laminaren Unterschicht an: 


p = Res(4— 3-3)" fur OS tes (23) 


b) Bekanntlich hat Prandtl dem logarithmischen Gesetz mittels der sog. Schubspannungs- 
geschwindigkeit uv, = r/o die universelle Gestalt 
u 1 
a= 


Ux Y 1 Ux Do 
In “= bel L+ Inga +e") 
geben kénnen. Fiir die Ablésungsstelle v, = 0 versagt diese Formulierung offensichtlich. 

Es ist nun bemerkenswert, daB den oben von uns abgeleiteten Gesetzmafigkeiten eine ent- 
sprechende Formulierung gegeben werden kann, die nunmehr auch an der Ablésungsstelle gilt und 
die universelle Gestalt des logarithmischen Gesetzes miteinschlieBt! Dies wird durch Einfiithrung 
der verallgemeinerten ,,Schubspannungsgeschwindigkeit“ 


* x 


eG Lodpoe 
mittels der Transformation eae ae ws : es = Re, V1 — 2 No 4 (25) 


[Y. oe ay ae 
eee oe d 


erreicht. Man erhalt dann fiir den Ubergang nach (21) 


i =p ae alee Vs =) 


u Ik dl v- y 
Ve 2 x Vig yo Vn Go 
v Vv 
Le aia eee V. Veo, —§.\ 
xt 2 BOY Ve RO ae i 0 Vay: Van, 
Tepe x (fs t+ Ya tae AN 4 Za pee “| 
_ 26 
yj ee a 1 (26) 
+ Se inlee eeae Z é a as 
PEPE a F SEE TEE 5 v 2 
v Vv v Vv 
V. VO V.8 
fur +J * “0 * 
v y < y 
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ca v dp/dx —jA/2 : 
Dabei ist F — P| -i= | und hat die Bedeutung eines Formparameters. 
Yo, (OW Re (« A 2 
6 1 


a No 


Fiir die vollturbulente Schicht erhalten wir nach (22) die angekiindigte Verallgemeinerung des 
logarithmischen Gesetzes in seiner universellen Gestalt 


ee Vay 
po Vier 


4 F V : v Py HD 
yj ut Fo+% Spe (27) 


* il 
i : a V. V,6 
<<. 


Bei schwachem Druckanstieg ergibt (27) mit V, ~v, und F ~ 0 die universelle Gestalt des log- 
arithmischen Gesetzes. 


Man kann auch dem Wandgesetz (27), das den Formparameter F enthalt, durch Einfiihrung der 
Variablen 


ee View, i (dp/dx) y 
Se Fee Ve (28) 


eine universelle Formulierung geben: 


u 2 = 1 Nile ici 9§ 1 
Seen se a] a heresy | — ; 
7,— Ce = 4, V+ E—1) + pn 4 $1 + nx) (29) 


hierbei ist die Konstante 
u 1 
wobei (u/V,,)) die Geschwindigkeit am Rande der laminaren Unterschicht bedeutet. 


Fiir die laminare Unterschicht lautet die Darstellung in den Koordinaten (25) 


uw ay Vg Ogee 
ee v A F| v ee! )|; (20) 
die fur F ~ 0 in die bekannte Relation 
u 0, ¥ 
Vy v 


ibergeht. 

c) Als derzeit empirische Groen treten im Formelapparat die Koeffizienten x und V, 09/y auf. 
Der Koeffizient x erscheint bereits im logarithmischen Gesetz. Experimentelle Analyse ergab seine 
Unabhangigkeit vom Druckgradienten!; dagegen ist x als von der Reynoldsschen Modellzahl ab- 
hangig anzusehen, nimmt jedoch fir sehr groBe Werte derselben den Charakter einer universellen 
Konstante an. Auch die von mir * angegebene Erweiterung des logarithmischen Gesetzes um den 
Einflu8 des Druckgradienten ergab mit den obigen Eigenschaften des Koeffizienten x gute Uber- 
einstimmung mit dem Experiment. 


Die genannten Eigenschaften des Koeffizienten x weist auch der Koeffizient v, 09/v des log- 
arithmischen Gesetzes auf. Wenn man bedenkt, da fiir schwacheren Druckanstieg V,, 09/v ~ Vx 09/, 
so ist es naheliegend anzunehmen, da auch V,, 6)/y dieselben charakteristischen Kigenschaften be- 
sitzen wird und also fiir ein vorliegendes Stromungsmodell eine Konstante vom Betrage der Konstan- 
ten v, 0)/v sein wird. Physikalisch laBt sich diese Hypothese so deuten: Die Variable v, y/y kann 
aufgefaBt werden als eine mit der Schubspannungsgeschwindigkeit gebildete lokale Reynoldssche 


1 Vgl. hierzu die Messungen von H. Ludwieg und W. Tillmann, Ing.-Arch. 17 (1949), S. 244. 
2 W. Szablewski, Ing.-Arch. 22 (1954), S. 268. 


302 Szablewski: Wandnahe Geschwindigkeitsverteilung usw. Ingenieur-Archiv 


Zahl in der laminaren Unterschicht, deren kritischer Wert an der Turbulenzgrenze mit der Dicke 69 
der laminaren Unterschicht v, 6,/v betragt. Dabei ist gemaB der Annahme t ~ 7) in der wand- 
nahen Zone die Schubspannungsgeschwindigkeit v, ~ /t/o . Im Bereich mittleren oder starken 
Druckanstiegs, den wir hier betrachten, ist die Annahme konstanter Schubspannung nicht mehr 
haltbar: der Schubspannungsverlauf in Wandnahe wird durch die vollstandigere Gleichung 


dp 
tm + hy 

beschrieben. Damit tritt nunmehr anstelle von v, do/y die mit 
_ |/%o aed 
Ve = ees 0 dx ° 
gebildete kritische Reynoldssche Zahl V, dy/¥ . 

Die Hypothese V,, d)/y = konst. fiir ein vorliegendes Strémungsmodell hat eine wichtige Konse- 
quenz. Der Formparameter F' nimmt an der Stelle der Ablésung, wenn man die sich hier aus (24) 


ergebende Identitat 


dp/dx _—il 
eVe  % 
beachtet, nach der in (26) gegebenen Definition den Wert 
1 
F 4, = V,, dof? (31) 


an und ist durch die Modellkonstante V, 6)/y festgelegt. Dieses Ergebnis findet in Abschnitt 4 eine 
gewisse experimentelle Bestatigung. 

Wir merken weiter an, dafs nach der Hypothese V, 69/vy = konst. sich die Dicke 6) der 
laminaren Unterschicht aus der Gleichung 


6o\2 c 
63 ar 2h 62 79 = y a) bestimmt. 


4. Vergleich mit dem Experiment.* Zum Vergleich ziehen wir die von Nikuradse ! in divergenten 
Kanalen durchgefiihrten Messungen heran. Es bedeuten im folgenden: 6 die Kanalbreite, « den 
halben Offnungswinkel des ebenen Diffusors, u, die Geschwindigkeit in Kanalachse. 

a) Ermittlung der Wandschubspannung. Wie wir schon kiirzlich 2 darlegten, 
kénnen die theoretischen GesetzmaBigkeiten der Geschwindigkeitsverteilung in Wandnahe dazu 
dienen, die schwer zu messenden sehr kleinen Wandschubspannungen tT, im Druckanstiegbereich zu 
ermitteln. In der zitierten Arbeit hatten wir diese Ermittlung durchgefihrt, indem wir die Kurve 
des ,,eingeschrankten“ Wandgesetzes * jeweils passend durch die wandnahen MeBpunkte der Ge- 
schwindigkeit legten. Wir fihren diese Prozedur nochmals mittels des Wandgesetzes (22) durch und 
werden gegeniiber den damals ermittelten Werten 7) kleine, mit dem Druckanstieg wachsende, 
Korrekturen erhalten. 

Das Wandgesetz (22) lautet hier mit 


en CEie und Re = —— 


if = 
A (61) ak ; 


Seria 


folgendermaBen: 


| / 2 (32) 
Gar |i Gi 4 ‘0 


ee 
Cpe Uy ae ee 
4 
ae 4 | +m (24 Re Ye, — 4 ny) + a 


v v 


* Die praktische Durchfiihrung dieses Abschnitts lag in den Handen meines Mitarbeiters Herrn H. Moch. 
* J. Nikuradse, Untersuchungen tiber die Strémung des Wassers in konvergenten und divergenten Kana- 
len, VDI-Forschungsheft 289 (1929). 


2 W. Szablewski, Ing.-Arch. 22 (1954), S. 268. 3 Vgl. FuBnote 1 von Seite 296. 
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wobei 


F = 


(Re) (a— mm) 


Die empirischen Koeffizienten x, V,, 6)/v sind hier bei der mittleren Reynoldsschen Zahl Re ~ 81500 
den Messungen der konvergenten Renae vuneen entnehmbar, die bei der nur wenig verschiede- 
nen Zahl Re ~ 67500 durchgefiihrt 


wurden: 


A me (M0. 


M9" get 40 
= 9,51, 


Unter Zugrundelegung der von Niku- 
radse gemessenen Druckgradienten A er- 
mittelten wir zunachst durch Anpassen 
(Abb. 1) die Werte des Parameters 


Mg oar ae i 
Die dimensionslosen Wandschubspan- Messung Nikuradse 
nungen c, konnten dann unter Benut- oe ° 0=7° 
d der Definitionsgleich om cee 
zung der aus der Definitionsgleichung Sm oe 
(25) fiir V, do/v flieBenden Beziehung sae 
V, 5p Wandgeserz 24lg-y 
A (33) : ne 1 ; 
(Re/2) a eg ai 
berechnet werden. Es ergaben sich folgende Werte: 
2 
& exp Vea No o <a 
Iie —0,043 0,0375 0,0062 0,00134 (137) 
2g —0,085 0,0332 '  0,0070 0,00095 (102) 
ase —0,120 0,0298 0,0078 0,00065 ( 75) 
4° —0,140 0,0241 0,00967 0,00026 ( 42) 


(Vgl. ferner Abb. 2 ?.) 


In Klammern haben wir die erwahnten friuheren Werte zum Vergleich eingetragen. 


b) Dicke der laminaren Unterschicht. Nach (33) ergibt sich der in Abb. 3? 
angegebene Verlauf der dimensionslosen Dicke 7 = 6,/(b/2) der laminaren Unterschicht uber a, 
der das Anwachsen der laminaren Unterschicht mit gréfer werdendem Druckanstieg zeigt. 


c) Wandnahe Geschwindigkeitsverteilung. In Abb. 4 bis 7 sind die 
Ubergangsverteilung (26) sowie das Wandgesetz (27) fiir folgende sich gema8 der Ermittlung in a) 


ergebenden Parameterwerte 
Cz) anh eae) ee a ee ee 
F | 0,0050 | 0,0142 | 0,0278 | 0,0615 
dargestellt [Formel fiir F siehe (32) ]. 
Die MeBpunkte wurden gema{ der Ermittlung in a) nach der (25) nay Formel 


es Pees». ar 27 = (Re/2) Va—4n ZN 


eingetragen. 


1 W. Szablewski, Ing.-Arch. 20 (1952), S. 37. 
2 In den Abb. 2 und 3 sind noch die Werte fiir « = 0° eingetragen: c, = 0,0018 (vgl. W. Sz., Ing.-Arch. 20) 
Vy Oo/? 


d 7) = 0,0055 i8 der Formel 7, = —*———. 
und 7% gemaB de ne (Re)2)Va 
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Es fallt auf, daB die wandnachsten Punkte schlecht zum Ubergangsbogen liegen. Diese Diskre- 
panz diirfte nicht zu Lasten der Theorie gehen; man kann vermuten, daB mit den damals (1929) zur 


To 
G 
Dicke der laminaren Unterschicht 
0,0020 —Wandschubspannung 0010 . ; 
e 
001 | 
@ 
e fo} $ 
e 
0010 + i 005 J 
e 
s a EL I 
C0005 o Impulssatz (Wikuradse) 
© Wandgesetz l 
WS alr ain el 0 af eee 
0 1 ° 2 °o 3 oO 4 ° & ° if te) 2 Oo 3 oO X o 5 ° 
Abb. 2. Abb. 3. 
we 
Ihe hn 
40 1s : 5 4 40 ii 
Geschwindigkeltsverteilung Geschwindigkeitsverrellung 
30 ° Messung Mkuradse 0c=7° 30 ° Messung Nikuradse Ce Se 
Ybergang of — lbergang eA 
——— Wandgesetz ——— Wandgesefz oe 
20 20 CE = 
10 ! 70 
| - 
laminar 
| , 
| te ins 
= -_ 
0 i, a Ps 0 if 2 3 
Abb. 4. Abb. 5. 
wb 
u Vi 
Veg i IF | 
1) —— ; Geschwindigkertsverteilung 
Geschwindigkeitsverteilung es 
30 2 Messung Nikuradse oc=8° | 
libergang 30 . Messung Nikuradse a=4° 
———Wandgesetz Ybergang 
——— Wondgesetz ye 
20 aE = 
me 
70 
\y 
ee | 
0 1 y 3 | 
Abb. 6, Abb. 7. 


Verfiigung stehenden Mitteln die Messung dieser wandnachsten Punkte, welche die Messung sehr 
geringer Druckdifferenzen erforderte, nur unvollkommen durchgefihrt werden konnte !. 


1 Vgl. hierzu H. Reichardt, Z. angew. Math. Mech. 20 (1940), S. 297. 
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Abb. 8 zeigt die Gesamtheit der Messung und der theoretischen Verteilung nach dem Wand- 
gesetz (27). Dazu haben wir noch das logarithmische Gesetz und die sich gemaB (31) mit 


1 
F= 5-57: (= 9,1053) 
ub _/b 
Vo. eh 
50 —— 
Geschwindigkeitsverteilung 
40 = 
Messung Nikuradse 
° o=7° 
° Ze 
30 L Ss 3° 
a ae 
Wanagesetz 
———— log. Geserz 
20 : 
70 


30 + 
Messung Mikuradse 
o 8 6@=7° 
e Ye | 
a aaeesl ° 40 + 
Qa o : 
o ue A, 
universelles Wandgeserz < 
e 
70 a 
o 
O° oe p 
025 
"e 
fae Gpfdz) 
. ee wary 
- 0 Z Zi 


Abb, 9. 


ergebende theoretische Verteilung fiir die Ablésung eingetragen. (Um den EinfluB des Druck- 
gradienten deutlicher hervorzuheben, baben wir hier (u/V,) — (u/V,,)) aufgetragen, wobei 


u\ _ Va bo/y FV d 
(7) i z 


v 
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nach (30) die Geschwindigkeit am Rande der laminaren Unterschicht angibt. ) Der Abb. 8 kann 
man entnehmen, inwieweit dem Einflu8 des Druckgradienten durch das Wandgesetz Rechnung 
getragen wird. Zur Verteilung fiir die Ablésung ist zu bemerken, daf ihre relative Lage gut zu der 
Angabe von Nikuradse paBt, daB die Ablésung bei a ~ 4,8° einsetzt. 

SchlieBlich enthalt Abb. 9 noch das Wandgesetz in seiner universellen Form (29). Auch hier 
zeigt sich eine befriedigende Ubereinstimmung zwischen Theorie und Experiment. 


5. Zusammenfassung. Unter Vernachlassigung der Tragheitsglieder wird aus der Differential- 
gleichung turbulenter Grenzschichtstrémungen bei Zugrundelegung des Prandtlschen Mischungs- 
wegansatzes, wobei der Mischungsweg dem Wandabstand proportional gesetzt wird, die wandnahe 
Geschwindigkeitsverteilung fiir Strémungen mit Druckanstieg bis einschlieBlich der Ablosung in 
geschlossener analytischer Form approximativ berechnet. Die Fehlerabschatzung, welche auf die 
Approximation eines elliptischen Integrals hinauskommt, ergibt einen fiir technische Stromungen 
vernachlassigbaren Fehler. Dem als Verallgemeinerung des logarithmischen Gesetzes erhaltenen 
Wandgesetz kann bei Einfiihrung einer verallgemeinerten Schubspannungsgeschwindigkeit eine die 
universelle Form des logarithmischen Gesetzes enthaltende Formulierung gegeben werden. Die 
empirischen Koeffizienten des Wandgesetzes sind die des logarithmischen Gesetzes. Der Vergleich 
mit den Nikuradseschen Messungen divergenter Kanalstrémungen ergibt befriedigende Uberein- 
stimmung. 


(Eingegangen am 8. Oktober 1954.) 
Anschrift des Verfassers: Dr. W. Szablewski, Deutsche Akademie der Wissenschaften, 
Forschungsinstitut fiir Mathematik, Abt. Angewandte Mathematik, Berlin W 8, Jagerstr. 22—23. 
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